quinta-feira, 3 de maio de 2012

METODO DE NEWTON


MÉTODO DE NEWTON
!
! NOME: METODO DE NEWTON
! AUTOR: DR. CARLOS SANTANA
! DATA.: 11/09/08 18:31C
!
      program newton
      real p0,tol,n0,p
      integer i
      i=1
      n0=10000.
      print*,'entre com p0'
      read*,p0
      print*,'entre com a tolerancia'
      read*,tol
      do while(float(i).le.n0)
      p=p0-f(p0)/df(p0)
      if(abs(p-p0).lt.tol )then
         print*,'procedimento ralizado com sucesso'
         print*,'Raiz  = ', p
      go to 100
      endif
      i=i+1
      p0=p
      enddo
      print*,'O procedimento falhou depois de', n0,'interacoes'
 100  stop
      end
      !!!!!!!!!!!!!!!!!!funcao!!!!!!!!!!!!
      real function f(x)
      f=exp(x)-3.*x**2. !se preciso troque a funcao
      return
      end
      !!!!!!!!!!!!!!!derivada!!!!!!!!!!!!!!!!!!!!
      real function df(x)
      df=exp(x)-6*x  !se preciso troque a derivada
      return
      end

Nenhum comentário: