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 26
| DO i=1, n
! do j=1,10
Y(i+1) = Y(i)+dY
! alpha=0.1
! alpha=alpha*j
C Le systeme differentiel descriptif du systeme
f1 =(X2(i)/X1(i))-(alpha*(X1(i)**(1+gama)))
f3= X2(i)
f2 = -exp(-X3(i))+(v0/X1(i))
C L'algorithme d'EULER
X1(i+1) = X1(i)+f1*dY ! Equation (1) du systeme
X2(i+1) = X2(i)+f2*dY ! Equation (2) du systeme
X3(i+1)=X3(i)+f3*dY !Equation (3) du systeme
end do
! write(*,*) REAL(Y(i)),REAL(X1(i))
" alpha c'est une parametre variable".
"cest un systeme d'equation monodimensionell suivant y"
C On sauvegarde les données pour tracer la courbe X1(i) = f(Y) mais on pourrait tracer d'autre courbes,en particulier X2(i)=f(Y),,. |
Partager