Programmation équa diff en scilab

Fermé
ss_maling Messages postés 19 Date d'inscription jeudi 2 avril 2009 Statut Membre Dernière intervention 19 avril 2009 - 9 avril 2009 à 18:49
Bonjour,
Oui t'as raison
J'ai repris le script que tu trouveras ci-après.Seulement j'ai un problème pour la sortie si je veux que le résultat affiché soit une matrice 5x5. en faisant par exemple M=ss_maling1(1:5), il affiche une matrice 16x16
function M=ss_maling1(n)
h=0.01;
k=0.1;
m=zeros(1,1);
p=zeros(1,1);
u=0:0.1:1.5;
[m,p]=size(u);
for i=1:p
e=u+0.1;
a=1/h^2+u/(2*h);
b=-2/h^2-2/k;
b(2:p)=b(1);
c=1/h^2+e/(2*h);
M=diag(a)+diag(b(1:p-1),1)+diag(c(1:p-2),2);
end
endfunction