Hello everyone.
I'd like to vectorize this code:
for i = 1:N-2
for j = 1:N
Jacobian(i,j) = D2Mat(i+1,j)+(2/XCheb(i+1))*...
D1Mat(i+1,j)+(xi^2)*n*Y(i+1)^(n-1)*D0Mat(i+1,j);
end
end
where D2Mat is (NxN), XCheb is (N,1), D1Mat is (NxN), D0Mat is (NxN), Y is (N,1) and xi, n are real values.
Thanks

 採用された回答

David Hill
David Hill 2020 年 10 月 18 日
編集済み: David Hill 2020 年 10 月 18 日

0 投票

Jacobian=D2Mat(2:end-1,:)+(2./XCheb(2:end-1)).*D1Mat(2:end-1,:)+((xi^2)*n*Y(2:end-1).^(n-1)).*D0Mat(2:end-1,:);

1 件のコメント

Jesus Octavio
Jesus Octavio 2020 年 10 月 18 日
Thank you !!!
It was so easy... 🤦‍♂️

サインインしてコメントする。

その他の回答 (0 件)

Community Treasure Hunt

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

Start Hunting!

Translated by