Modelo Gráfico y Punto Fijo
eduardo19cocoTarea10 de Abril de 2020
708 Palabras (3 Páginas)126 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)) ...