Metododos Numericos Ejercicios - 11-15
Metododos Numericos Ejercicios - 11-15
Metododos Numericos Ejercicios - 11-15
5=80
a) En forma analítica
x 3.5=80
3.5
√ x3.5 =3.5√ 80
x=3.497357243
b) Con el método de la falsa posición dentro de ε a=2.5 %, los valores iniciales de 2 a 5
x 3.5=80
ITERACIÓN 1:
Evaluamos x i=2 ; xu =5 en f ( x )
f ( x i ) =(2)3.5−80=−68,6862
f ( x u ) =(5)3.5−80=199.5084
f ( x r )=(2.768317204)3.5−80=−44.70156253
f ( x i )∗f ( x r ) > 0
ITERACIÓN 2:
f ( x i ) =(2.768317204)3.5 −80=−44.70156253
f ( x u ) =(5)3.5−80=199.5084
f ( x r )=(3.176816952)3.5−80=−22.8557325
f ( x i )∗f ( x r ) > 0
ITERACIÓN 3:
f ( x i ) =(3.176816952)3.5−80=−22.8557325
f ( x u ) =(5)3.5−80=199.5084
Encontramos x r y evaluamos x r en f (x)
f ( x r )=(3.36421311)3.5 −80=−10.16192727
f ( x i )∗f ( x r ) > 0
ITERACIÓN 4:
f ( x i ) =(3.36421311)3.5−80=−10.16192727
f ( x u ) =(5)3.5−80=199.5084
f ( x r )=(3.44349351)3.5−80=−4.229973213
f ( x i )∗f ( x r ) > 0
5.12 Dada
Δx 1
n=log 2 ( )
εa
=log 2 ( )
0.05
=4.321928 ≈ 5
x=1:0.01:2;
y=-2*(x.^6)-1.5*(x.^4) +10*x+2;
>> plot(x,y)
ITERACIÓN 1:
Encuentro x r
x u−x i 1+2
x r= = =1.5
2 2
Evaluamos x i=1 ; xu =2; x r =1.5 en f ( x )
calculo f ( x i )∗f ( x r )
ITERACIÓN 2:
Encuentro x r
x u−x i 1.5+1
x r= = =1.25
2 2
Evaluamos x i=1 ; xu =1.5; x r =1.25 en f ( x )
calculo f ( x i )∗f ( x r )
ITERACIÓN 3:
Encuentro x r
x u−x i 1.25+1.5
x r= = =1.375
2 2
Evaluamos x i=1.25 ; xu =1.5 ; xr =1.375 en f ( x )
calculo f ( x i )∗f ( x r )
ITERACIÓN 4:
Encuentro x r
x u−x i 1.25+1.375
x r= = =1.3125
2 2
Evaluamos x i=1.25 ; xu =1.375 ; x r =1.375 en f ( x )
calculo f ( x i )∗f ( x r )
Δx 1
n=log 2 ( )
εa
=log 2
0.001( )
=9.9657 ≈ 10
>> x=2:0.01:5;
y= (9.8/15) *x(1-exp(-135/x));
plot(x,y)
ITERACIÓN 1:
f ( x i )∗f ( x r ) < 0
ITERACIÓN 2:
f ( x i )∗f ( x r ) < 0
5.14 Se carga una viga de la manera que se aprecia en la figura P5.14. Emplee el método de bisección
para resolver la posición dentro de la viga donde no hay momento.
ITERACIÓN 1:
x u−x i 2.5+ 3
x r= = =2.75
2 2
ITERACIÓN 2:
x u−x i 2.75+ 2.5
x r= = =2.625
2 2
ITERACIÓN 3:
x u−x i 2.5+ 2.625
x r= = =2.5625
2 2
ITERACIÓN 4:
x u−x i 2.5625+ 2.625
x r= = =2.59375
2 2
xrnuevo−X ranterior 2.59375−2.5625
ε a= = =0.01204
xnuevo
r
2.59375
ITERACIÓN 5:
x u−x i 2.59375+ 2.625
x r= = =2.609375
2 2
ITERACIÓN 6:
x u−x i 2.59375+ 2.609375
x r= = =2.6015625
2 2
ITERACIÓN 7:
x u−x i 2.6015625+2.609375
x r= = =2.60546875
2 2
ITERACIÓN 8:
x u−x i 2.60546875+2.609375
x r= = =2.607421875
2 2
ITERACIÓN 9:
x u−x i 2.607421875+2.609375
x r= = =2.608398438
2 2
xrnuevo−X ranterior 2.608398438−2.607421875
ε a= = =0.000374
xnuevo
r
2.608398438
ITERACIÓN 10:
x u−x i 2.608398438+2.609375
x r= = =2.608886719
2 2
ITERACIÓN 11:
x u−x i 2.608886719+2.608398438
x r= = =2.608642579
2 2
5.15. Por un canal trapezoidal fluye agua a una tasa de Q=20 m3 / s. La profundidad critica y para dicho
canal satisface la ecuación
Q2
0=1− B
g A 3c
y2
B=3+ y y A c =3 y +
2
Resuelva para la profundidad critica con el uso de los métodos a) grafico b) bisección y c) falsa posición.
En los incisos b) y c), haga elecciones iniciales de x i=0.5 y xu =2.5 y ejecute iteraciones hasta que el
error aproximado caiga por debajo del 1 % o el número de iteraciones supere a 10. Analice sus resultados
a) Grafico
(20)2
0=1− (3+ y )
y2 3
(9.81)(3 y + )
2
(20)2
1= ( 3+ y)
6 y + y2 3
( 9.81)( )
2
( 20 )2
1= 3
( 3+ y )∗8
( 9.81 ) ( 6 y + y 2 )
(20)2∗8
¿
9.81
(20)2 (6 y+ y 2)3
0= −
9.81 24 +8 y
x=2:0.01:5;
y=40.77471967-((6*x+(x.^2))^3)/(24+8*x);
plot(x,y)
(6 y + y 2)3
0=40.77471967−
24+8 y
b) Bisección
3
(6 y + y2 )
f ( x )=40.77471967− ; ε a=1 %=0.01 ; x i=0.5 ; x u=2.5
24+ 8 y
f ( x i )∗f ( x u )=−7012.3606<0 existe una raiz
ITERACIÓN 1:
Encuentro y r
y u − y i 0.5+ 2.5
yr = = =1.5
2 2
Evaluamos y i=0.5 ; y u =2.5; y r=1.5 en f ( x )
2 3
f ( y i) =40.77471967−
[ 6 ( 0.5 ) + ( 0.5 ) ]
=39.54871521
24 +8 ( 0.5 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 ) ]
=−177.3094423
24 +8 ( 0.5 )
2 3
f ( yr ) =40.77471967−
[ 6 ( 1.5 ) + ( 1.5 ) ]
=1.22393842
24+8 ( 0.5 )
calculo f ( y i)∗f ( y r )
ITERACIÓN 2:
Encuentro y r
y u − y i 1.5+2.5
yr = = =2
2 2
Evaluamos y i=1.5 ; y u=2.5 ; y r =2 en f ( x )
3
f ( y i) =40.77471967−
[ 6 ( 1.5 )+ ( 1.5 )2 ] =1.22393842
24 +8 ( 0.5 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 ) ]
=−177.3094423
24 +8 ( 0.5 )
2 3
f ( yr ) =40.77471967−
[ 6 ( 2 ) +( 2 ) ]
=−61.6252
24+8 ( 2 )
calculo f ( y i)∗f ( y r )
y nuevo
r − y anterior
r 2−1.5
ε a= nuevo
= =0.25>0.01
yr 2
ITERACIÓN 3:
Encuentro y r
y u − y i 1.5+2
yr = = =1.75
2 2
Evaluamos y i=1.5 ; y u=2 ; y r =1.75 en f ( x )
2 3
f ( y i) =40.77471967−
[ 6 ( 1.5 )+ ( 1.5 ) ]
=1.22393842
24 +8 ( 0.5 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 2 )+ ( 2 ) ]
=−61.6252
24 +8 ( 2 )
2 3
f ( yr ) =40.77471967−
[ 6 ( 1.75 ) + ( 1.75 ) ]
=−24.8754
24+8 ( 1.75 )
calculo f ( y i)∗f ( y r )
y nuevo
r − y anterior
r 1.75−2
ε a= nuevo
= =0.1428>0.01
yr 1.75
ITERACIÓN 4:
Encuentro y r
y u − y i 1.5+1.75
yr = = =1.625
2 2
Evaluamos y i=1.5 ; y u=1.75 ; y r =1.625 en f ( x )
2 3
f ( y i) =40.77471967−
[ 6 ( 1.5 )+ ( 1.5 ) ]
=1.22393842
24 +8 ( 0.5 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 1.75 ) + ( 1.75 ) ]
=−24.8754
24 +8 ( 1.75 )
2 3
f ( yr ) =40.77471967−
[ 6 ( 1. 625 ) + ( 1.625 ) ]
=−10.6388
24+8 ( 1. 625 )
calculo f ( y i)∗f ( y r )
y nuevo
r − y anterior
r 1. 625−1.75
ε a= nuevo
= =0.07692>0.01
yr 1.625
ITERACIÓN 5:
Encuentro y r
y u − y i 1.5+1.625
yr = = =1.5625
2 2
Evaluamos y i=1.5 ; y u=1.625 ; y r =1.5625en f ( x )
2 3
f ( y i) =40.77471967−
[ 6 ( 1.5 )+ ( 1.5 ) ]
=1.22393842
24 +8 ( 0.5 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 1.625 )+ (1. 625 ) ]
=−10.6388
24 +8 ( 1.625 )
2 3
f ( yr ) =40.77471967−
[ 6 ( 1. 5625 )+ ( 1.5625 ) ]
=−4.42797
24 +8 ( 1. 5625 )
calculo f ( y i)∗f ( y r )
y nuevo
r − y anterior
r 1.5625−1. 625
ε a= nuevo
= =0.04> 0.01
yr 1. 5625
ITERACIÓN 6:
Encuentro y r
y u − y i 1.5+1.5625
yr = = =1.53125
2 2
Evaluamos y i=1.5 ; y u=1.5625 ; y r =1.5625en f ( x )
3
f ( y i) =40.77471967−
[ 6 ( 1.5 )+ ( 1.5 )2 ] =1.22393842
24 +8 ( 0.5 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 1.5625 ) + ( 1. 5625 ) ]
=−4.42797
24+ 8 ( 1.5625 )
2 3
f ( yr ) =40.77471967−
[ 6 ( 1. 53125 )+ ( 1.53125 ) ]
=−1.5341
24+8 ( 1. 53125 )
calculo f ( y i)∗f ( y r )
y nuevo
r − y anterior
r 1.53125−1. 5625
ε a= nuevo
= =0.0204>0.01
yr 1. 53125
ITERACIÓN 7:
Encuentro y r
y u − y i 1.5+1.53125
yr = = =1.516075
2 2
Evaluamos y i=1.5 ; y u=1.53125 ; y r =1.516075en f ( x )
2 3
f ( y i) =40.77471967−
[ 6 ( 1.5 )+ ( 1.5 ) ]
=1.22393842
24 +8 ( 0.5 )
3
f ( yu ) =40.77471967−
[ 6 ( 1.53125 ) + ( 1. 53125 )2 ] =−1.5341
24+ 8 ( 1.53125 )
2 3
f ( yr ) =40.77471967−
[ 6 ( 1. 516075 )+ ( 1.516075 ) ]
=−0.1781
24 +8 ( 1. 516075 )
calculo f ( y i)∗f ( y r )
y nuevo
r − y anterior
r 1.516075−1.53125
ε a= nuevo
= =0.0100009>0.01
yr 1. 516075
ITERACIÓN 8:
Encuentro y r
y u − y i 1.5+1.516075
yr = = =1.5080375
2 2
Evaluamos y i=1.5 ; y u=1.516075 ; y r =1.5080375en f ( x )
3
f ( y i) =40.77471967−
[ 6 ( 1.5 )+ ( 1.5 )2 ] =1.22393842
24 +8 ( 0.5 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 1.516075 ) + ( 1. 516075 ) ]
=−0.1781
24+ 8 ( 1.516075 )
2 3
f ( yr ) =40.77471967−
[ 6 ( 1. 5080375 )+ ( 1.5080375 ) ]
=0.5272
24 +8 ( 1. 5080375 )
calculo f ( y i)∗f ( y r )
y nuevo
r − y anterior
r 1.5080375−1.516075
ε a= nuevo
= =0.005329<0.01
yr 1. 5080375
c) Falsa Posición
3
(6 y + y2 )
f ( x )=40.77471967− ; ε a=1 %=0.01 ; x i=0.5 ; x u=2.5
24+ 8 y
f ( x i )∗f ( x u )=−7012.3606<0 existe una raiz
f ( y i) =40.77471967−
[ 6 ( 0.5 ) + ( 0.5 ) ]
=39.54871521
24 +8 ( 0.5 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 ) ]
=−177.3094423
24 +8 ( 0.5 )
Encontramos y r y evaluamos y r en f ( y )
f ( y u ) ( y u− y i) (−177.3094423 )( 2.5−0.5 )
y r = y u− =2.5− =0.8647427025
f ( y u )−f ( y i ) −177.3094423−( 39.54871521 )
3
f ( xr ) =40.77471967−
[ 6 ( 0.8647427025 ) + ( 0.8647427025 )2 ] =34.00861398
24 +8 ( 0.8647427025 )
f ( x i )∗f ( x r ) > 0
ITERACIÓN 2:
f ( y i) =40.77471967−
[ 6 ( 0.8647427025 )+ ( 0.8647427025 ) ]
=34.00861398
24 +8 ( 0.8647427025 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 ) ]
=−177.3094423
24 +8 ( 0.5 )
Encontramos y r y evaluamos y r en f ( y )
f ( y u ) ( y u− y i) (−177.3094423 ) (2.5−0.8647427025 )
y r = y u− =2.5− =1.127913939
f ( y u )−f ( y i ) −177.3094423−( 34.00861398 )
3
f ( xr ) =40.77471967−
[ 6 ( 1.127913939 )+ ( 1.127913939 )2 ] =25.03871005
24 +8 ( 1.127913939 )
f ( x i )∗f ( x r ) > 0
y nuevo
r − y anterior
r 1.127913939−0.8647427025
ε a= nuevo
= =0.23332> 0.01
yr 1.127913939
ITERACIÓN 3:
f ( y i) =40.77471967−
[ 6 ( 1.127913939 )+ (1.127913939 ) ]
=25.03871005
24 +8 ( 1.127913939 )
3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 )2 ] =−177.3094423
24 +8 ( 0.5 )
Encontramos y r y evaluamos y r en f ( y )
f ( y u ) ( y u− y i) (−177.3094423 ) (2.5−1.127913939 )
y r = y u− =2.5− =1.297696883
f ( y u )−f ( y i ) −177.3094423−( 25.03871005 )
2 3
f ( xr ) =40.77471967−
[ 6 ( 1.297696883 )+ ( 1.297696883 ) ]
=16.07161846
24 +8 ( 1.297696883 )
f ( x i )∗f ( x r ) > 0
y nuevo
r − y anterior
r 1.297696883−1.127913939
ε a= nuevo
= =0.12083> 0.01
yr 1.297696883
ITERACIÓN 4:
f ( y i) =40.77471967−
[ 6 ( 1.297696883 )+ (1.297696883 ) ]
=16.07161846
24 +8 ( 1.297696883 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 ) ]
=−177.3094423
24 +8 ( 0.5 )
Encontramos y r y evaluamos y r en f ( y )
f ( y u ) ( y u− y i) (−177.3094423 ) (2.5−1.297696883 )
y r = y u− =2.5− =1.397618545
f ( y u )−f ( y i ) −177.3094423−( 16.07161846 )
2 3
f ( xr ) =40.77471967−
[ 6 ( 1.397618545 )+ ( 1.397618545 ) ]
=9.359925276
24 +8 ( 1.397618545 )
f ( x i )∗f ( x r ) > 0
y nuevo
r − y anterior
r 1.397618545−1.297696883
ε a= nuevo
= =0.07149> 0.01
yr 1.397618545
ITERACIÓN 5:
f ( y i) =40.77471967−
[ 6 ( 1.397618545 )+ (1.397618545 ) ]
=9.359925276
24 +8 ( 1.397618545 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 ) ]
=−177.3094423
24 +8 ( 0.5 )
Encontramos y r y evaluamos y r en f ( y )
f ( y u ) ( y u− y i ) (−177.3094423 ) (2.5−1.397618545 )
y r = y u− =2.5− =1.45289386
f ( y u )−f ( y i ) −177.3094423−( 16.07161846 )
3
f ( xr ) =40.77471967−
[ 6 ( 1.45289386 ) + ( 1.45289386 )2 ] =5.13419048
24 +8 ( 1.45289386 )
f ( x i )∗f ( x r ) > 0
y nuevo
r − y anterior
r 1.45289386−1.397618545
ε a= nuevo
= =0.038044>0.01
yr 1.45289386
ITERACIÓN 6:
Evaluamos y i=1.45289386 ; y u=2.5 en f ( x )
2 3
f ( y i) =40.77471967−
[ 6 ( 1.45289386 )+ ( 1.45289386 ) ]
=5.13419048
24 +8 ( 1.45289386 )
2 3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 ) ]
=−177.3094423
24 +8 ( 0.5 )
Encontramos y r y evaluamos y r en f ( y )
f ( y u ) ( y u− y i ) (−177.3094423 ) (2.5−1.45289386 )
y r = y u− =2.5− =1.482360727
f ( y u )−f ( y i ) −177.3094423−5.13419048
2 3
f ( xr ) =40.77471967−
[ 6 ( 1.482360727 ) + ( 1.482360727 ) ]
=2.722426367
24 +8 ( 1.482360727 )
f ( x i )∗f ( x r ) > 0
y nuevo
r − y anterior
r 1.482360727−1.45289386
ε a= nuevo
= =0.0198>0.01
yr 1.482360727
ITERACIÓN 7:
f ( y i) =40.77471967−
[ 6 ( 1.45289386 )+ ( 1.45289386 ) ]
=2.722426367
24 +8 ( 1.45289386 )
3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 )2 ] =−177.3094423
24 +8 ( 0.5 )
Encontramos y r y evaluamos y r en f ( y )
f ( y u ) ( y u− y i ) (−177.3094423 ) (2.5−1.482360727 )
y r = y u− =2.5− =1.497749381
f ( y u )−f ( y i ) −177.3094423−2.722426367
2 3
f ( xr ) =40.77471967−
[ 6 ( 1.497749381 )+ (1.497749381 ) ]
=1.417449702
24 +8 ( 1.497749381 )
f ( x i )∗f ( x r ) > 0
y nuevo
r − y anterior
r 1.497749381−1.482360727
ε a= nuevo
= =0.0102745> 0.01
yr 1.497749381
ITERACIÓN 8:
f ( y i) =40.77471967−
[ 6 ( 1.497749381) + ( 1.497749381 ) ]
=1.417449702
24+ 8 ( 1.497749381)
2 3
f ( yu ) =40.77471967−
[ 6 ( 2.5 ) + ( 2.5 ) ]
=−177.3094423
24 +8 ( 0.5 )
Encontramos y r y evaluamos y r en f ( y )
f ( y u ) ( y u− y i) (−177.3094423 ) (2.5−1.497749381 )
y r = y u− =2.5− =1.505698044
f ( y u )−f ( y i ) −177.3094423−1.417449702
2 3
f ( xr ) =40.77471967−
[ 6 ( 1.505698044 ) + ( 1.505698044 ) ]
=0.730959265
24+8 ( 1.505698044 )
f ( x i )∗f ( x r ) > 0
y nuevo
r − y anterior
r 1.505698044−1.497749381
ε a= nuevo
= =0.005279< 0.01
yr 1.505698044