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

 Réponse acceptée

David Hill
David Hill le 18 Oct 2020
Modifié(e) : David Hill le 18 Oct 2020

0 votes

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 commentaire

Jesus Octavio
Jesus Octavio le 18 Oct 2020
Thank you !!!
It was so easy... 🤦‍♂️

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

Community Treasure Hunt

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

Start Hunting!

Translated by