Bonjour,

Je dois faire une fonction qui me calcul le jacobien, et qui me donne et sortie la matrice sous cette forme bien précise :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
J=[0,1;-2000*x(1)*x(2)-1,1000-1000*x(1)^2];
Je calcul le jacbien de la manière suivante :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
syms x1 x2
h=[x2,1000*(1-x1^2)*x2-x1];
v=[x1 x2];
J=jacobian(h,v);
Le problème c'est que j'obtiens le jacobien avec des x1 et x2 au lieu de l'avoir sous forme vectorielle avec x(1) et x(2). Ca paraît bête mais j'ai essayé beaucoup de truc, notamement avec "subs" mais rien n'y fait...

Est-ce que quelqu'un saurait comment à la fois calculer le jacobien dans ma fonction, et avoir la sortie sous la bonne forme?