Jacobi method giving an error

2 次查看（过去 30 天）
Isaac Al-rai
Isaac Al-rai 2018-2-17
评论： Image Analyst 2018-2-17
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 个评论
显示 1更早的评论
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 Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by