Bonjour voila je dois intégrer une matrice qui se résume simplement a intégrer chaqu'un des termes de la matrice. Cependant voila il n'i a pas de fonction toute faite pour intégrer une matrice et il faut pour chaque intégrale appeler le terme de la matrice et dire de l'intégrer, résultat je doit me taper 18 ligne de commande répétitive ce qui m'enchante pas vraiment. Je souhaiterais réslaiser le calcul a l'aide d'une boucle For faisant intervenir le i et j correspondant a la ligne et la colonne de la matrice. Voila ce que j'ai pour l'instant :
q:=evalm(linalg[matrix](4,1,[psi1,psi2,psi3,psi4]));
> qT:=evalm(linalg[matrix](1,4,[psi1,psi2,psi3,psi4]));
> Phi:=evalm(q&*qT);
>
> S1:=22485.44;
> L:=450;
> psi1:=Z->16/L^4*Z^4-32/L^3*Z^3+16/L^2*Z^2;
> psi2:=Z->16/L^4*Z^5-40/L^2*Z^4+32/L^2*Z^3-8/L*Z^2;
> psi3:=Z->-24/L^5*Z^5+52/L^4*Z^4-34/L^3*Z^3+7/L^2*Z^2;
> psi4:=Z->4/L^4*Z^5-8/L^3*Z^4+5/L^2*Z^3-1/L*Z^2;
> Phi[1,1]:=int(psi1(z)*psi1(z),z=0..168 .6);
> Phi[2,2]:=int(psi2(z)*psi2(z),z=0..168 .6);
> Phi[3,3]:=int(psi3(z)*psi3(z),z=0..168 .6);
> Phi[4,4]:=int(psi4(z)*psi4(z),z=0..168 .6);
> Phi[1,2]:=int(psi1(z)*psi2(z),z=0..168 .6);
> Phi[1,3]:=int(psi1(z)*psi3(z),z=0..168 .6);
> Phi[1,4]:=int(psi1(z)*psi4(z),z=0..168 .6);
> Phi[2,3]:=int(psi2(z)*psi3(z),z=0..168 .6);
> Phi[2,4]:=int(psi2(z)*psi4(z),z=0..168 .6);
> Phi[3,4]:=int(psi3(z)*psi4(z),z=0..168 .6);
> Phi[4,3]:=int(psi4(z)*psi3(z),z=0..168 .6);
> Phi[2,1]:=int(psi2(z)*psi1(z),z=0..168 .6);
> Phi[3,1]:=int(psi3(z)*psi1(z),z=0..168 .6);
> Phi[4,1]:=int(psi4(z)*psi1(z),z=0..168 .6);
> Phi[3,4]:=int(psi3(z)*psi4(z),z=0..168 .6);
> Phi[3,2]:=int(psi3(z)*psi2(z),z=0..168 .6);
> Phi[4,2]:=int(psi4(z)*psi2(z),z=0..168 .6);
Ma fonction psi est fonction de z.
En gros je souhaite obtenir un truc du genre :
q:=evalm(linalg[matrix](4,1,[psi1,psi2,psi3,psi4]));
> qT:=evalm(linalg[matrix](1,4,[psi1,psi2,psi3,psi4]));
> Phi:=evalm(q&*qT);
>
> S1:=22485.44;
> L:=450;
for i j from ..........
Phi[i,j]:=int(psi[i](z)*psi[j](z),z=0..168.6);
Cordialement.
-----