bonjour a tous
je dois resoudre sous matlab 3 equations differentielles . Prenons la premeire l'equation est :x’(t)=-tx(t)+cos(at) , une equation du premier ordre non lineaire .
Sous matlab j'ai programé:(en utilisant la methode d 'euler)

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
clear all
 
n=100
a=0;%borne inferieure%
b=5 ;%borne superieure%
a1=0;
b1=-3;
 
t(1)=0.1; %condition initiale%
x(1)=0.1;%condition initiale%
u(1)=0.1; %condition initiale%
v(1)=0.1;
 
h=(b-a)/n%pas%
h1=(b1-a1)/n
 
 
for i=2:n
  t(i)=t(i-1)+h%
  x(i)=x(i-1)+h*(cos(t(i-1))-t(i-1)*x(i-1))%dérivée de x par rapport par à t%
  u(i)=u(i-1)+h1%incrementation de t par le pas%
  v(i)=v(i-1)+h1*(cos(u(i-1))- u(i-1)*v(i-1))%dérivée de x par rapport par à t%
end
 
plot (t,x,u,v)
Pouvez vous me dire deja si cette resolution est juste et m'expliquer l'interpretation graphique de resultat(point singuliers,points max, chalmps des tangeantes)

Merci , j'attens vos suggestions.