programme fortran Fermé

Question

Bonjour,
je dois résoudre une équation en utilisant la methode de newton mais en programmation fortran 90
et g besoin d'aide!!voici la fonction f(x)=x**3 + 3x - 1.
merci d'avance

juanpablo · Accepted Answer

bonsoir

je suppose que tu dois programmer le code de cette methode et non pas utiliser une librairie...


c'est une methode iterative pour trouver la racine de f(x)=0

xn+1=  xn - f(xn)  / f'(xn)

il faut donc definir une procedure pour calculer f f'

subroutine F_FP(x,f,fp)

f=x**3 + 3x - 1. 

fp=3*x**2 + 3  

end

 
et  faire une boucle avec des criteres  de controle  


Program MethodeNewton
!-----declaration 
      
!----valeur de depart
      x=1
!----boucle
      do i=1, imax
       
        call F_FP(x,f,fp)

!---test derivé fp
         
!-- iteration
        xi= x - f/fp      

!--test convergence  xi-x     

!--x accepté       
       x=xi      
    enddo

!--afficher x  si  ok
     if  ....  
        write(*,*) ' solution x=
     else
        write(*,*) ' erreur ...  
     endif
   end


voila c'est juste le principe, à toi de detailler le code

DL@ · Answer

Bonsoir,
De mémoire :
par exemple, si on veut le résultat dans une variable Y,
Y=X**3  + X*2 - 1 (** pour la puissance, et * pour la multiplication)
DL

sma · Answer

c Program MethodeNewton 
real x 

  x=1 
  do i=1, imax 
  call F_FP(x,f,fp) 
1 xi= x - f/fp 

  if (abs(xi-x).gt.10-3) then 
   x=xi 
goto 1
  else 
write(*,*) ' solution x=', x
endif 
enddo
end 
 

subroutine F_FP(x,f,fp)
real x
f=x**3 + 3x - 1 
fp=3*x**2 + 3 
return
end

Programme fortran

3 réponses

Discussions similaires