Trabajito Estadistica

Descargar como xlsx, pdf o txt
Descargar como xlsx, pdf o txt
Está en la página 1de 5

f(x)=x^3+2x^2+10x-20

𝑔(𝑥)=20/(𝑥^2+2𝑥+10)
x f(x)
-3 -59
-2 -40 x g(x)
-1 -29 1 1.53846154 pert[1,2]
0 -20 2 1.11111111 pert[1,2]
1 -7
2 16 Hay por lo menos un punto fijo
3 55
Valor inicial:Xo= 1.5
|g'(x)<1| x g'(x)
1.5 -0.4299919377 CONVERGE

ITER Xi EA
0 1.500000 ------------ tol=0.0001
1 1.311475 0.188525
2 1.394416 0.082941
3 1.357476 0.036941
4 1.373844 0.016369
5 1.366574 0.007270
6 1.369800 0.003226
7 1.368368 0.001432
8 1.369003 0.000635
9 1.368721 0.000282
10 1.368847 0.000125
11 1.368791 0.000056 4.0587503E-05
Si ER<0.5*10^(-n)
4.09057767 4 0.00005
5 0.000005
sol=1.368
𝑥_1=𝑥_0−(𝑓(𝑥_0))/(𝑓′
MÉTODO DE NEWTON RAPHSON (𝑥_0))
Xo 1.5
iter Xi f(xi) f'(xi) EA
0 1.5 2.875 22.75 ------------- tol=0.0001
1 1.373626 0.10178868 21.1550537 0.12637363
2 1.368815 0.00014159 21.0962213 0.00481155
3 1.368808 2.75083E-10 21.0961393 6.71179E-06

𝑥_2=𝑥_1−(𝑥_1−𝑥_0)/(𝑓(𝑥_1 )−𝑓(𝑥_0 ) )
MÉTODO DE LAS SECANTES 𝑓(𝑥_1) tol=0.0001

iter X0 X1 f(xo) f(x1) X2 EA


0 1 2 -7 16 1.30434783 ---------------
1 2 1.30434783 16 -1.33475795 1.3579123 0.053564479
2 1.30434783 1.3579123 -1.33475795 -0.22913573 1.36901333 0.011101021
3 1.3579123 1.36901333 -0.22913573 0.00432957 1.36880746 0.000205866
4 1.36901333 1.36880746 0.00432957 -1.3672E-05 1.36880811 6.48061E-07
solucion
1.3688
B 3 2 4 3.6
Q 12 20 13 30
H 1.19425807 2.52778031 1.01881963

𝑥_1=𝑥_0−(𝑓(𝑥_0))/
(𝑓′(𝑥_0)) 𝐹(𝐻)=3.33(𝐵−0.2𝐻) 𝐻^(3/2)−𝑄
METODO DE NEWTON RAPHSON

iter Xi f(xi) f'(xi) EA


0
1
2
3
𝐹(𝐻)=3.33(𝐵−0.2𝐻) 𝐻^(3/2)−𝑄

%EC BAZIN CON NEWTON RAPHSON


clc,clear
B=3,Q=12
f=@(H) 3.33*(B-0.2*H)*H^(1.5)-Q
df=@(H) -333/500*H^(3/2)+3*sqrt(H)*(3.33B-333/500)/2
ezplot(f)
grid on
tol=0.0001,EA=10
x0=input('Xo=')
disp(' raiz aprox EA')
while(EA>tol)
x1=x0-f(x0)/df(x0);
EA=abs(x1-x0);
fprintf("%12.8f %14.9f \n",x1,EA)
x0=x1;
endwhile

También podría gustarte