• Jueves 18 de Abril de 2024, 16:41

Autor Tema:  metodo steffensen?  (Leído 4647 veces)

peppe

  • Nuevo Miembro
  • *
  • Mensajes: 1
    • Ver Perfil
metodo steffensen?
« en: Miércoles 13 de Mayo de 2009, 17:32 »
0
Buenas tengo un problema, nos piden el emtodo de setffesen en fortrain, y el codigo me da error.

Cual seria el codigo de metodo de steffensen?

dodne esta mi error?

la funcion es: f(Xn+1)=(Xn-(f(xn))^2)/(f(xn+f(xn))-F(Xn))



       program steffensen


      real xo, fxo, fxo2, sum, fxo3, xi, xixo, fxi
      integer iter
     
      Write(*,*) 'Escriu el valor de xo:'
      Read(*,*) xo

      iter = 1
     
      do i = 1,10


fxo = cos(xo) - xo
         fxo2 = fxo * fxo
         sum = xo + fxo
         fxo3 = cos(sum) - sum


         xi = xo - fxo2/(fxo3-fxo)

         xixo = xi - xo
         xixo = ABS(xixo)
         fxi = cos(xi) - xi

         if ( xixo .lt. 0.0000005 ) then


            if ( fxi .lt. 0.000001 ) then
               solucio = xi
               Write(*,*) 'La solucio es xi igual a ',xi
            end if
         else
             iter = iter + 1
             xo = xi
             write(*,*)xo
         endif
         
      end do
     
      end