¡Descarga Búsqueda de raíces de funciones no lineales: Métodos de Newton y Bisección - Prof. Gonzále y más Ejercicios en PDF de Física solo en Docsity!
Hoja ecuaciones no lineales. Búsqueda de ceros de funciones
1.Obligatoria
program Newton implicit none integer i,n real*8 x0,x1,f,fprima,funcion,derivada,norma,s
write (,)'Valor de x' read(,)x
do i=1,
funcion=f(x0) derivada=fprima(x0) x1= x0-(funcion/derivada) !Es la sumatoria del cuadrado de la diferencia if(abs(f(x1).le.1e-10)goto 100
x0=x !x(i+1)=x(i)-(x(i)3-2.d0x(i)-5.d0)/(3.d0(x(i)**2)-x(i)) enddo
100 continue write(,)'la solucion es: ', x
pause stop end
real8 function f(x0) funcion problema implicit none integer n,i real8 x
do i=1,n f=x0*3-2.d0x0-5.d
enddo !write(,)'Funcion', f return
end
real8 function fprima(x0) !primera derivada de la funcion problema implicit none integer n,i real8 x
do i=1,n fprima=3.d0(x02)-2.d enddo !write(,*)'derivada', fprima return
End
program Biseccion implicit none integer k,n real*8 f,mult1,a,b,c,c0,mult,norma,funcion
100 continue write(,)'introduzca el valor del intervalo' read(,) a,b
mult=f(a)f(b) !write(,)'f(a)f(b)', mult !comprobaci¢n de la multiplicacion
if (mult.lt.0.d0)then c=dble(a+b)/2.d else if (mult.gt.0.d0)then goto 100 endif
!write(,)c,f(c) do k=1,
mult1=f(a)*f(c) norma=abs(f(c))
if (mult1.lt.0.d0)then b=c !goto 200
