Jacobi method giving an error

2 views (last 30 days)
Isaac Al-rai
Isaac Al-rai on 17 Feb 2018
Commented: Image Analyst on 17 Feb 2018
For a numerical methods class I am in I have to create a function that solve an arbitrarily system of linear equations using the jacobi method. I am testing out a code from my textbook
%
function X=jacobi(A,B,P,delta,max1)
N = length(B);
if true
% code
end
for k=1:max1
for j=1:N
X(j)=(B(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j);
end
err=abs(norm(X'-P));
relerr=err/(norm(X)+eps);
P=X';
if (err<delta)||(relerr<delta)
break
end
end
if true
% code
end
X=X';
end
end
Where delta is my tolerance and max1 is number of iterations.
And when I enter it (after defining my A,B, and P terms) I am giving the following error:
if true
% code
X=@jacobi(A,b,P,0.001,1);
X=@jacobi(A,b,P,0.001,1);
Error: Unbalanced or unexpected parenthesis or bracket.
end
What does this mean? Dont I need to put the function in parenthesis in order for it to run?

Accepted Answer

Image Analyst
Image Analyst on 17 Feb 2018
Get rid of the @ symbol. Just call it by it's name alone.
X = jacobi(A,b,P,0.001,1);
  3 Comments
Isaac Al-rai
Isaac Al-rai on 17 Feb 2018
Thought I might ask this too, is there a way I can make this function run using only one for loop?
Image Analyst
Image Analyst on 17 Feb 2018
Looks like you might be able to do it with conv() or conv2(), but I've not gotten into the details to figure it out. Maybe you can.

Sign in to comment.

More Answers (0)

Categories

Find more on MATLAB in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!