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

En savoir plus sur MATLAB dans Centre d'aide et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by