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