je cherche a comprendre un programme bien compliqué qui est une methode numerique de resolution d'equation aux derivé partiellles.
au passage le coup de la grille je n'ai rien compris bref)
voici mon programme si quelqu'un connait un peu ce logiciel :
function [tnew,U]=transport(t,x,u0,v)
tnew=0;
vectv=v(x)';
vectu=u0(x);
U=vectu;
for i=1:length(t)-1
vectu=[vectu(1)-(dt/dx)*vectv(1)*(vect(2)-vectu(1)),
vectu(2:$-1)-(dt/(2*dx))*vectv(2:$-1)+(vectu(3:$)-vect1:$-2),
vectu($)-(dt/dx)*vectv($)*(vectu($)-vectu($-1))]
if modulo(t(i+1),1)==0
then U=[U,vectu];
tnew=[tnew,t(i+1)];
end
end
endfunction
et le message d'erreur :
then U=[U,vectu];
!--error 34
Syntaxe de l'instruction de contrôle incorrecte.
at line 11 of function transport called by :
endfunction
at line 15 of exec file called by :
voila voila j'aimerais bien reussir a le faire tourner pour mieu comprendre ce qu'on y fait
ps: si vous ne trouver pas l'erreur mais que vous avez toute expliquation concernat ce programme je suis prenneur
merci