Salut, j'ai trouvé un système non linéaire en utilisant la méthode crank nicolson comment puis-je le résoudre en utilisant la fonction fsolve mon système dépend de temps et l'espace . (si vous avez des exemple, des cours ..)
merci pour m'aider.
for j=1:Nt+1
for i=2:Nx
Un(i,j)=Un(i,j-1)+D*(perm(xc(i+1),n)*aa(Un(i+1,j))+(perm(xc(i),n)*aa(Un(i,j))))*(Un(i+1,j)-Un(i,j))...
-D*(perm(xc(i),n)*aa(Un(i,j))+(perm(xc(i-1),n)*aa(Un(i-1,j))))*(Un(i,j)-Un(i-1,j))...
+D*(perm(xc(i+1),n)*aa(Un(i+1,j-1))+(perm(xc(i),n)*aa(Un(i,j-1))))*(Un(i+1,j-1)-Un(i,j-1))...
-D*(perm(xc(i),n)*aa(Un(i,j-1))+(perm(xc(i-1),n)*aa(Un(i-1,j-1))))*(Un(i,j-1)-Un(i-1,j-1));
end
end
Partager