Modelo Gráfico y Punto Fijo
Enviado por eduardo19coco • 10 de Abril de 2020 • Tareas • 708 Palabras (3 Páginas) • 78 Visitas
REPÚBLICA BOLIVARIANA DE VENEZUELA
MINISTERIO DEL PODER POPULAR PARA LA DEFENSA
UNIVERSIDAD NACIONAL EXPERIMENTAL POLITÉCNICA
DE LA FUERZA ARMADA NACIONAL BOLIVARIANA
UNEFA-NÚCLEO NUEVA ESPARTA
INGENIERIA CIVIL 03S - 1344 - D1
[pic 1]
Modelo
Gráfico y Punto Fijo
Profesor: Jesús Suniaga Bachiller:
Albert Maldonado.
27.261.528
Juan griego, Abril del 2018
Comandos en R.
Pto. Fijo
> g=function(x)(x*log(x)/0.1)^0.5
> curve(g,from=-10,to=10)
> abline(h=0)
[pic 2]
> x0=1.37
> tol=0.0001
> n=20
> for (i in 1:20)
+ {x=g(x0)
+ print (c(i, x0,x))
+ x=g(x0)
+ print (c(i, x0,x))
+ e=abs(x-x0)
+ if (e [1] 4.00000 7.27862 12.01983 [1] 4.00000 7.27862 12.01983 [1] 5.00000 12.01983 17.28815 [1] 5.00000 12.01983 17.28815 [1] 6.00000 17.28815 22.19720 [1] 6.00000 17.28815 22.19720 [1] 7.00000 22.19720 26.23177 [1] 7.00000 22.19720 26.23177 [1] 8.00000 26.23177 29.27430 [1] 8.00000 26.23177 29.27430 [1] 9.00000 29.27430 31.44055 [1] 9.00000 29.27430 31.44055 [1] 10.00000 31.44055 32.92569 [1] 10.00000 31.44055 32.92569 [1] 11.00000 32.92569 33.91912 [1] 11.00000 32.92569 33.91912 [1] 12.00000 33.91912 34.57315 [1] 12.00000 33.91912 34.57315 [1] 13.00000 34.57315 34.99934 [1] 13.00000 34.57315 34.99934 [1] 14.00000 34.99934 35.27523 [1] 14.00000 34.99934 35.27523 [1] 15.00000 35.27523 35.45307 [1] 15.00000 35.27523 35.45307 [1] 16.00000 35.45307 35.56740 [1] 16.00000 35.45307 35.56740 [1] 17.00000 35.56740 35.64077 [1] 17.00000 35.56740 35.64077 [1] 18.00000 35.64077 35.68781 [1] 18.00000 35.64077 35.68781 [1] 19.00000 35.68781 35.71794 [1] 19.00000 35.68781 35.71794 [1] 20.00000 35.71794 35.73723 [1] 20.00000 35.71794 35.73723 Regla Falsa > a=1; b=2; N=20; E=0.0001 > f=function(x) 0.1*x^2-x*log(x) > curve(f, from=0.1,to=10,col="red") > abline(h=0,col="blue") [pic 3] > if (f(a)*f(b) > 0) { + stop()} > xant=0 > for ( i in 1:N){ + xm=(b*f(a)-a*f(b))/(f(a)-f(b)) + err=abs(xm-xant) + cat("\n","iteracion",i) + cat("\n","a,b,xm",a,b,xm,"imagenes",f(a),f(b),f(xm), "\n", "raiz=", xm, "error =",err) + if(abs(f(xm)) ...