methode de resolution d'equations différentielles par methode de MEEM
صفحة 1 من اصل 1
methode de resolution d'equations différentielles par methode de MEEM
! methode de resolution d'equations différentielles par methode de MEEM
implicit none
!------------------------partie declaration------------
integer::i,n
parameter(n=1000)
real,dimension(n):yi,yi1,xi
real::xp,W,f,fp,f2p
!------------------------partie saisie-------------
x0=0.0
y0=1.0
!------------------------partie tgraitement-------------
xi(1)=x0
yi(1)=y0
W=(xi-xp)/n
do i=1,n
Yi(1)=Yi+W*f(xi,yi)+(w*w/2)fp(xi,yi)+(w*w*w/6)*f2p(xi,yi)
enddo
!------------------------partie affichage---------------
write(*,*)'le resultat final est',Yi1
end
!------------------------bloc fonctionnel----------------
function f(xi,yi) result(f1)
implicit none
real xi,yi,f1
f1=-xi*yi
end
function fp(xi,yi) result(fp1)
implicit none
real xi,yi,fp1
fp1=yi*(xi*xi-1)
end
function f2p(xi,yi) result(fp2)
implicit none
real xi,yi,fp2
fp2=xi*yi*(3-xi*xi)
end
implicit none
!------------------------partie declaration------------
integer::i,n
parameter(n=1000)
real,dimension(n):yi,yi1,xi
real::xp,W,f,fp,f2p
!------------------------partie saisie-------------
x0=0.0
y0=1.0
!------------------------partie tgraitement-------------
xi(1)=x0
yi(1)=y0
W=(xi-xp)/n
do i=1,n
Yi(1)=Yi+W*f(xi,yi)+(w*w/2)fp(xi,yi)+(w*w*w/6)*f2p(xi,yi)
enddo
!------------------------partie affichage---------------
write(*,*)'le resultat final est',Yi1
end
!------------------------bloc fonctionnel----------------
function f(xi,yi) result(f1)
implicit none
real xi,yi,f1
f1=-xi*yi
end
function fp(xi,yi) result(fp1)
implicit none
real xi,yi,fp1
fp1=yi*(xi*xi-1)
end
function f2p(xi,yi) result(fp2)
implicit none
real xi,yi,fp2
fp2=xi*yi*(3-xi*xi)
end
مواضيع مماثلة
» un program qui permet de calculer integral par methode de simpson
» un program qui permet de calculer les integrales par la methode de trapez
» un program qui permet de calculer les integrales par la methode de trapez
صفحة 1 من اصل 1
صلاحيات هذا المنتدى:
لاتستطيع الرد على المواضيع في هذا المنتدى