Jacobi method giving an error

1 次查看(过去 30 天)
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?

采纳的回答

Image Analyst
Image Analyst 2018-2-17
Get rid of the @ symbol. Just call it by it's name alone.
X = jacobi(A,b,P,0.001,1);
  3 个评论
Isaac Al-rai
Isaac Al-rai 2018-2-17
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 2018-2-17
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.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Argument Definitions 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by