3 RTMNR2



Compilación de programa


C
C PROGRAMA PARA CALCULAR UNA RAIZ DE UNA ECUACION TRASCENDENTE
C CON EL METODO DE NEWTON DE SEGUNDO ORDEN
C
100 FORMAT (F10.4,2I5)
101FORMAT ("1 ITERACION APROXIMACION RELATIVA VALOR APROXIMADO
1DE LA RAIZ")
102 FORMAT (1 X,45 (1H-))
103 FORMAT (2 6X, E20.10)
104 FORMAT (I 6,F20.10,E20.10)
105 FORMAT("EL METODO NO CONVERGE EN",I5,"ITERACIONES")
106 FORMAT ("UNA RAIZ DE LA ECUACION ES",E1.10)
F (X)= X**2.0-0.5
DF(X)=2.0*X
DDF(X)=2.0
OPEN (UNIT=1,FILE="DATOSRTMNR2",STATUS="OLD")
OPEN (UNIT=2,FILE="RESULTADO RTMNR2",STATUS="UNKNOWN")
READ (1,*) X,NMI,NAPROX
APROX=10.**(-NAPROX)
WRITE (*,101)
WRITE (2,101)
WRITE (*,102)
WRITE(2,102)
WRITE (*,103) X
WRITE (2,103) X
DO 1 L=1,NMI
Y=X+1.0/(0.5*DDF(X)/DF(X)-DF(X)/F(X))
REL=ABS((X-Y)/Y)
WRITE (*,104) L,REL,Y
WRITE (2,104)L,REL,Y
IF(REL.LE.APROX) GO TO 2
1 X=Y
WRITE(*,102)
WRITE(2,102)
WRITE(*,105)NMI
WRITE(2,105)NMI
GO TO 3
2 WRITE(*,102)WRITE(2,102)
WRITE(*,106)Y
WRITE(2,106)Y
3 CALL EXIT
END

Datos de inicio de programa


Valor aproximado de la raíz
X
0.6
Número Máximo deIteraciones
NMI
10
Número de cifras significativas exactas con que se desea el resultado
NAPROX
8





Resultados de programa


1 ITERACION APROXIMACION RELATIVA VALOR APROXIMADODELA RAIZ
---------------------------------------------
0.6000000238E+00
1 0.1505375952 0.7063291073E+00
20.0010997799 0.7071067691E+00
3 0.0000000000 0.7071067691E+00
---------------------------------------------
UNA RAIZ DE LA ECUACION ES 0.7071067691E+00