Page 1 sur 1

valeurs et vecteurs propres avec maxima

MessagePosté: Samedi 12 Janvier 2013, 12:37
par lassaib
Bonjour, je veux calculer les valeurs et vecteurs propres d'une matrice ( a*v = e*s*v) où s est une matrice, e et v les valeurs et vecteurs propres respectifs. Mon code tourne mais ça ne termine pas les calculs sur la dernière commande.
Y a t il un grain de sable dans mon code ? Merci de votre aide.
Code: Tout sélectionner

/* [wxMaxima: comment start ]
Methode des valeurs propres
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
ze : 1.75$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
f : exp(-ze*s)*[1,s,u,s**2,t**2,u**2,s*u,u*t**2,u**3,(t*u)**2]$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Element de matrice de recouvrement
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
s1[n,m]:= f[n]*f[m]*u*(s**2-t**2)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
s[n,m]:=integrate(integrate(integrate(s1[n,m],t,0,u),u,0,s),s,0,inf)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Energie cinetique
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
k1[n,m]:=(diff(f[n],s)*diff(f[m],s)+diff(f[n],t)*diff(f[m],t)+
diff(f[n],u)*diff(f[m],u))*u*(s**2-t**2)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
k[n,m]:=integrate(integrate(integrate(k1[n,m],t,0,u),u,0,s),s,0,inf)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
p1[n,m]:=(diff(f[n],s)*diff(f[m],u)+diff(f[n],u)*diff(f[m],s))*s*(u**2-t**2) +
(diff(f[n],t)*diff(f[m],u)+diff(f[n],u)*diff(f[m],t))*t*(s**2-u**2)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
p[n,m]:=integrate(integrate(integrate(p1[n,m],t,0,u),u,0,s),s,0,inf)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Energie potentielle
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
z : 2$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
v1[n,m]:=f[n]*f[m]*(s**2-t**2-4*z*s*u)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
v[n,m]:=integrate(integrate(integrate(v1[n,m],t,0,u),u,0,s),s,0,inf)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Hamiltonien
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
h[n,m]:=k[n,m]+p[n,m]+v[n,m]$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Creation des matrices
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
hh : genmatrix(h,10,10);
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
ss : genmatrix(s,10,10);
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Matrice final
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
iss : invert(ss)$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: input   start ] */
aa : iss*hh$
/* [wxMaxima: input   end   ] */
/* [wxMaxima: comment start ]
Valeurs et vecteurs propres de aa vals : solve([determinant(aa)=0],[e]);
   [wxMaxima: comment end   ] */
/* [wxMaxima: input   start ] */
[vals,vecs] : eigenvectors(aa);
/* [wxMaxima: input   end   ] */