Alumna: Neira Alvarado Rosario Alexandra

NEWTON RAPHSON

 CODIGO

1. f(x)=x2-2x-3 Donde: x0 =2

f´(x) = 2x-2

x1=2

f (2)
x2=x1-
f ´(2)

−(2)2−2(2)−3
x2=2
2(2)−2
x2= 3.5

( 3.5)2−2(3.5)−3
x3= 3.5−
2(3.5)−2
x3=3.05
 n x F(x) Fp(x)
1 3.5000000000 2.2500000000 5.0000000000
2 3.0500000000 0.2025000000 4.1000000000
3 3.0006097561 0.0024393962 4.0012195122
4 3.0000000929 0.0000003717 4.0000001858

3.0000000929 es una aproximación a la raíz con un error máximo de 0.0010000000

2. f(x)=sin(x)+3cos(x) Donde: x0 =1

f´(x) = cos(x)-3 sin(x)

x1=1

f (1)
x2=x1-
f ´ (1)
−sin ( 1 ) +3 cos ⁡( 1)
x2=1
cos ( 1 ) −3 sin ⁡( 1)
x2= 2.24104868

−sin ( 2.24 )+ 3 cos ⁡( 2.24)

x3=2.24
cos ( 2.24 )−3 sin ⁡( 2.24)
x3=1.877855073

n x f(x) fp(x)
1 2.2410486805 -1.0798863700 -2.9721785660
2 1.8777170886 0.0468942030 -3.1619299381
3 1.8925479684 -0.0000034381 -3.1622776602
4 1.8925468812 0.0000000000 -3.1622776602

1.8925468812 es una aproximacion a la raiz con un error maximo de 0.0001000000

3. f(x)=x2-5x+4 Donde: x0 =6

f´(x) = 2x-5

x1=6
 f (6)
x2=x1-
f ´ (6)

−(6)2−5 ( 6 ) +4
x2=2
2(6)−5
x2= 1.42

(1.42)2−5 ( 1.42 )+ 4
x3= 1.42−
2(1.42)−5
x3=0.918

n x F(x) Fr(x)
1 4.5714285714 2.0408163265 4.1428571429
2 4.0788177340 0.2426654372 3.1576354680
3 4.0019673693 0.0059059785 3.0039347387
4 4.0000012885 0.0000038655 3.0000025770
5 4.0000000000 0.0000000000 3.0000000000

4.0000000000 es una aproximacion a la raiz con un error maximo de 0.0010000000

METODO DE LA BISECCION
 CODIGO

4. f(x)=3x+sen(x)- e x Intervalo:(0,1)
 (1−0)
x1=1+
2
x1=1.5

f(x1) =3(1.5)+sen(1.5)- e 1.5

f(x1) =1.015

(1.015−1)
x2=1+
1.015
x2=1.014

f(x2) =3(1.014)+sen(1.014)- e 1.014

f(x2) =1.134

k a b c F(c) error
1.0000 0 1.0000 0.5000 0.3307 0.5000
2.0000 0 0.5000 0.2500 -0.2866 0.2500
3.0000 0.2500 0.5000 0.3750 0.0363 0.1250
4.0000 0.2500 0.3750 0.3125 -0.1219 0.0625
5.0000 0.3125 0.3750 0.3438 -0.0420 0.0313
6.0000 0.3438 0.3750 0.3594 -0.0026 0.0156
7.0000 0.3594 0.3750 0.3672 0.0169 0.0078
8.0000 0.3594 0.3672 0.3633 0.0071 0.0039
9.0000 0.3594 0.3633 0.3613 0.0023 0.0020
10.0000 0.3594 0.3613 0.3604 -0.0002 0.0010

f(c)=-0.00018
error= 0.00098

5. f(x)=x3-1 Intervalo:(-1,2)

(2+1)
x1=1+
2
x1=2.5

f(x1) =(2.5)3-1

f(x1) =14.625

(14.625−1)
x2=1+
14.625
x2=2.068

f(x2) = (2.069)3-1

f(x2) =7.8568
 k a b c F(c) error
1.0000 -1.0000 2.0000 0.5000 -0.8750 1.5000
2.0000 0.5000 2.0000 1.2500 0.9531 0.7500
3.0000 0.5000 1.2500 0.8750 -0.3301 0.3750
4.0000 0.8750 1.2500 1.0625 0.1995 0.1875
5.0000 0.8750 1.0625 0.9688 -0.0909 0.0938
6.0000 0.9688 1.0625 1.0156 0.0476 0.0469
7.0000 0.9688 1.0156 0.9922 -0.0233 0.0234
8.0000 0.9922 1.0156 1.0039 0.0118 0.0117
9.0000 0.9922 1.0039 0.9980 -0.0058 0.0059
10.0000 0.9980 1.0039 1.0010 0.0029 0.0029
11.0000 0.9980 1.0010 0.9995 -0.0015 0.0015
2.0000 0.9995 1.0010 1.0002 0.0007 0.0007
13.0000 0.9995 1.0002 0.9999 -0.0004 0.0004
14.0000 0.9999 1.0002 1.0001 0.0002 0.0002
15.0000 0.9999 1.0001 1.0000 -0.0001 0.0001

f(c)=-0.00009
error= 0.00009

6. f(x)=4x-7 Intervalo:(1,2)

(2−1)
x1=1+
2
x1=1.5

f(x1) =41.5-7

f(x1) =1

(1−1)
x2=1+
1
x2=0

f(x2) =40-7

f(x2) =1

k a b c F(c) error
1.0000
2.0000
3.0000
4.0000
5.0000
6.0000
7.0000
8.0000
9.0000
10.0000
7. f(x)=4x-13.9 Intervalo:(1,2)

(2−1)
x1=1+
2
x1=1.5

f(x1) =41.5-13.9

f(x1) =-5.9

(1.5−1)
x2=1+
1.5
x2=1.33

f(x2) =41.33-13.9

f(x2) =-2.5863

FALSA POSICION
 CODIGO
8. f(x)= 2x3+x-1 Intervalo (0,1)
f(a)= 2(0)3+(0)-1=0
f(b)= 2(1)3+(1)-1=2
0 ( 2 )−1(0)
x1=
2−(0)
x1=0
f(x1) =2(0)3+(0)-1
f(x1) =-1
0 ( 3 )−1(−1)
x2=
2−(−1)
x2=0.3333
f(x2) = 2(0.33)3+(0.33)-1
f(x2) =-0.592612
9. f(x)= sin(x)+3x-5 Intervalo (2,4)
f(a)= sin (2) +3(2)-5=1.909297
f(b)= sin (4) +3(4)-5=6.243197
2 ( 6.24 )−4 (−1.90)
x1=
6.24−(−1.90)
x1=2
f(x1) = sin (2) +3(2)-5
f(x1) =1.9193
2 ( 6.24 )−4 (−1.91)
x2=
6.24−(−1.91)
x2=2.000
f(x2) = sin (2) +3(2)-5
f(x2) =1.9193

METODO DE LA SECANTE

 CODIGO
10. f(x)=xlog(x)-10

xi=8

xi+1=9

f(x1)= (8)log(8)-10=-2.7752801

f(x2)= (9)log(9)-10=-1.4118174

−1.41 ( 9−8 )
x2=9- =10.03676471
−1.41−(−2.77)
f(x1)= (8)log(8)-10=-2.7752801

f(x2)= ( 10.03676471 )log

10.03676471 )-10=0.05276073
¿
x2= 10.03676471 -

(10.03676471−9) ( 0.05276073 )
=9.999436
0.05276073−(−1.4118174 )
11. f(x)= x10-1

xi=0

xi+1=1.3

f(x1)= (0)10-1=-1

f(x2)= (1.3)10-1=12.785849

x2=1.3-
12.785849 ( 1.3 )
=0.0943
12.785849−(−1)
f(x1)= (0)10-1=-1

f(x2)= ( 0.0943 )10-1=-0.99999

(0.09431−1.3) (−0.99999 )
x2= 0.09431 - =7988.768701
−0.99999−(−12.785849)
 METODO DEL PUNTO FIJO

 CODIGO
 12. f ( x )=¿ cos(x)
cos ( x )=x
valor inicial=0
 13. f ( x )=x 2−2 x−3
3
g ( x) =
x−2
valor inicial=4
EJERCICIOS CHAPRA

5.1. Determine las raíces reales de f ( x )=−0.5 x 2+2.5 x +4.5 :

a) Gráficamente
b) Empleando la fórmula cuadrática
c) Usando el método de bisección con tres iteraciones para determinar la raíz más
grande. Emplee como valores iniciales x l=5 y x u=10. Calcule el error
estimado Ea y el error verdadero Et para cada iteración.

PRIMERA INTERACION

5+10
x r=
2
x r=7.5
2
f ( x 1 )=−0.5 ( 7.5 ) +2.5 ( 7.5 ) +4.5
f ( x 1 )=−4.875
 SEGUNDA INTERACION

5+7.5
x r=
2
x r=6.25
2
f ( x 1 )=−0.5 ( 6.25 ) +2.5(6.25)+ 4.5
f ( x 1 )=0.59375

ERROR

|
Ea =
0.59375−(−4.875)
0.59375 | 100%= 7.21%
 5.2.Determine las raíces reales de f ( x )=5 x 3−5 x 2+ 6 x −2 :
a) Gráficamente
b) Utilizando el método de bisección para localizar la raíz más pequeña. Use los valores
iniciales x l=0 y x u=1 iterando hasta que el error estimado Ea se
encuentre debajo de ES =10 .
5.3. Determine las raíces reales de f ( x )=−25182 x −90 x2 + 44 x 3−8 x 4 +0.7 x 5
a) Gráficamente
b) Usando el método de bisección para localizar la raíz más grande con es
ES =10 . Utilice como valores iniciales x l=0.5 y x u=1.0
c) Realice el mismo cálculo que en b), pero con el método de la falsa posición y
ES =0.2 .
 5.4. Calcule las raíces reales de f ( x )=−12−21 x+ 18 x 2−2.75 x 3
a) Gráficamente
b) Empleando el método de la falsa posición con un valor es correspondiente a tres cifras
significativas para determinar la raíz más pequeña.
 LA RAIZ MAS PEQUEÑA ES: 2.2198

5.5. Localice la primera raíz no trivial de sen x =x2 , donde x está en radianes. Use una
técnica gráfica y bisección con un intervalo inicial de 0.5 a 1. Haga el cálculo hasta que
Ea sea menor que ES =2 . Realice también una prueba de error sustituyendo la
respuesta final en la ecuación original.
5.6. Determine la raíz real de ln x 2=0.7
a) Gráficamente
b) Empleando tres iteraciones en el método de bisección con los valores iniciales
x l=0.5 y x u=2 .
 c) Usando tres iteraciones del método de la falsa posición, con los mismos valores iniciales
de b).

*METODO DE LA BISECCION
*METODO DE LA PALSA POSICION
ALUMNO: Luiggi Gonsalo Quispe Rodriguez

1. METODO DE LA BESECCION : PROGRAMA

clear all % se borra texto del comand pront

format short;
a=input('Introduzca el valor de a: '); % se asigna el intervalo inferior
b=input('Introduzca el valor de b: '); % se asigna el intervalo superior
cont=input('Introduzca el número de iteraciones cont: '); % se asigna el numero de iteraciones
fun=input('Introduzcal a funcion f(x)=','s'); % se asigna la funcion
f=inline(fun);
for k=1:cont % aqui es donde se realiza la iteracion
c=(a+b)/2; % aqui se calcula el valor del punto medio
e=abs((b-a)/2); % aqui se calcula el error absoluto estimado
 A(k,:)=[k a b c f(c) e]; % aqui se imprime los valores de la iteracion
if f(a)*f(c)<0
b=c; % si es que es menor a 0 se encuentra en el intervalo
else
a=c; % si es es mayor que 0 se encontro la raiz
end
end
%las siguientes lineas son para imprimir la tabla de valores
fprint('\n \tk \ta \tb \tc \t(c) \terror \n')
disp(A)
fprint('Solución:\n c=%8.5f\n',c)
fprint('f(c)=%8.5f\n',f(c))
fprint('error=%8.5f\n',e)

EJEMPLO 1:

0.5*exp(x/3)-sin(x) (1;2) 0.001

-) Primera iteración

1+ 2
xi= =1.5
2

-) Segunda iteración

2+ 1.5
xi= =1.75
2

Introduzca el valor de a:
1
Introduzca el valor de b:
2
Introduzca el número de iteraciones cont:
7
Introduzcal a funcion f(x)=
0.5*exp(x/3)-sin(x)

Solución:
c= 1.91406
f(c)= 0.00471
error= 0.00781

EJEMPLO 2:

x^4+x^2+1 (1;2) 0.001

-)Primera iteración

1+ 2
xi= =1.5
2

-)Segunda iteración

1.5+2
xi= =1.75
2

Introduzca el valor de a:
1
Introduzca el valor de b:
2
Introduzca el número de iteraciones cont:
8
Introduzcal a funcion f(x)=
x^4+x^2+1

Solución:
c= 1.99609
f(c)=20.85976
error= 0.00391

2. FALSA POSICION:PROGRAMA

Fx=input('Ingrese la funcion: ','s'); % se ingresa la funcion

a=input('Ingrese a : '); %se ingresa el intervalo inferior
c=input('Ingrese c : '); %se ingresa el intervalo superior
e=input('Ingrese el error : '); % se ingresa la tolerancia

x=a; % se asigna la variable a -> x

Fa=eval(Fx); %se evalua la funcion para ver si tiene algun error con el intervalo inferior
x=c; % se asigna la variable a -> c
Fc=eval(Fx);%se evalua la funcion para ver si tiene algun error con el intervalo superior
fprint('\n %6s %7s %8s %10s %8s %8s %8s \n ','A','B','C','F(a)','F(b)','F(c)','|c-a|'); %se genera
las cabeceras de la tabla
while abs(c-a)>e % se genera la iteracion hasta el valor sea mayor a la tolerancia
b=(c*Fa-a*Fc)/(Fa-Fc); %Una la funcion evaulada en el intervalo inferior con el superior
x=b;
Fb=eval(Fx);%Evalua la funcion generada en b
fprint('\n %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f \n',a,b,c,Fa,Fb,Fc,abs(c-a)); %se
generan las iteraciones generadas
if abs(Fc)<e % si el valor absoluto de la funcion c menor que la tolerancia
break;
else
if Fa*Fb<=0 %si la funcion de a por b es menor e igual a 0 se vuelve a evaluar
c=b;
Fc=Fb;
else %si no se genera la nueva iteracion
a=b;
Fa=Fb;
end
end
end
fprint('\nEl resultado sera %.4f',b);
ezplot(Fx);%graficamos la funcion
grid on;

EJEMPLO 1:

X^2-2 (0;2) 0.001

Ingrese la funcion:
x^2-2
Ingrese a :
0
Ingrese c :
2
Ingrese el error :
0.001

-) Primera iteración
−1−8−2
0 (−1 )−1(−2)
x 1= ¿=1
¿
-) Primera iteración
1 (−0.22 )−1.33(−1)
x 2= =1.33
−0.22−(−1)

El resultado sera 1.4142

EJEMPLO 2:

exp(-1)-log x (1;2) 0.001

Ingrese la funcion:
exp(-1)-log(x)
Ingrese a :
1
Ingrese c :
2
Ingrese el error :
0.001

-) Priemra iteración
1 (−0.05 )−1.53 ( 0.36 )
x 2= =1.00
−0.05−(0.36)

-) Segunda iteración
 1 ( 1.45 )−1.45(−1)
x 2= =1.00
−0.09−( 0.36)

El resultado sera 1.4447

3. METODO DE LA SECANTE – PROGRAMACION

function xs = secante(fun,a,b,tol)
fprint('Método de la secante\n\n');
i = 1; %contador
fa = feval(fun, a); %se evalua la funcion de intervalo inferior
fb = feval(fun, b);%se evalua la funcion de intervalo superior
xs = b - ((b - a) / (fb - fa))*fb; % se calcula la derivada de la funcion
error = abs(b - a); % se genera la aproximacion de error absoluto

%se genera la cabecera de la tabla de valores

fprint('Iter. \t a \t \t b \t \t Xs \n');
fprint('%2i \t %f \t %f \t %f \n', i, a, b, xs);

while error >= tol % se genera la iteracion mientras que el error sea mayor o igual a la
tolerancia
b = a;
a = xs;
fb = feval(fun,b);%se evalua la funcion de intervalo inferior
fa = feval(fun,a);%se evalua la funcion de intervalo superior
xs = b - ((b - a)/(fb - fa))*fb;% se calcula la derivada de la funcion
error = abs(b - a); % se genera la aproximacion de error absoluto
i = i + 1;% se genera la siguiente iteracion
fprint('%2i \t %f \t %f \t %f \n', i, a, b, xs); %se generan los datos de la tabla de valores
end
w = feval(fun,xs);
fprint('\n La mejor aproximación a la raiz tomando una tolerancia de %f es \n x = %f con \n
f(x)= %f\n y se realizaron %i iteraciones\n',tol, xs, w, i);
end

EJEMPLO 1:

f(x)=xlog(x)-10

xi=8
xi+1=9
f(x1)= (8)log(8)-10=-2.7752801
f(x2)= (9)log(9)-10=-1.4118174
 −1.41 ( 9−8 )
x2=9- =10.03676471
−1.41−(−2.77)
Primera iteracion
f(x1)= (8)log(8)-10=-2.7752801
f(x2)= ( 10.03676471 )log 10.03676471 )-
¿
10=0.05276073
Segunda iteracion
x2= 10.03676471 -
(10.03676471−9) ( 0.05276073 )
=9.999436
0.05276073−(−1.4118174 )

Iter a b Xa
1 5.000000 9.000000 5.886429
2 5.886429 8.000000 5.738306
3 5.738306 5.886429 5.728972
4 5.728972 5.738306 5.728926
5 5.728926 5.728972 5.728926

Ejemplo 2

f(x)= x10-1

xi=0

xi+1=1.3
f(x1)= (0)10-1=-1

f(x2)= (1.3)10-1=12.785849

Primera iteracion:

12.785849 ( 1.3 )
x2=1.3- =0.0943
12.785849−(−1)

f(x1)= (0)10-1=-1

f(x2)= ( 0.0943 )10-1=-0.99999

Segunda iteracion:

(0.09431−1.3) (−0.99999 )
x2= 0.09431 - =7988.768701
−0.99999−(−12.785849)

Iter a b Xa
1 0.000000 1.300000 0.094320
2 0.094300 0.000000 169525.6
3 0.094358 0.094089 0.094356
4 0.093000 0.169595 0.094300
5 0.092125 0.094300 182696.265

4. METODO DE NEWTON - PROGRAMACION

x0=input('Ingrese el valor inicial: ');
tol=input('Ingrese el porcentaje de error: ');
f=input('Ingrese la función: ');
i=1;
fx(i)=x0;

syms x;
f1=subs(f,x,fx(i));
z=diff(f);
d=subs(z,x,fx(i));
ea(1)=100;

while abs(ea(i))>=tol;
fx(i+1)=fx(i)-f1/d; f1=subs(f,x,fx(i+1)); d=subs(z,x,fx(i+1));
ea(i+1)=abs((fx(i+1)-fx(i))/fx(i+1)*100);
i=i+1;
end
fprint('i fx(i) Error aprox (i) \n');
for j=1:i;
fprint('%2d \t %11.7f \t %7.3f \n',j-1,fx(j),ea(j));
end

EJEMPLO 1:
x^3-3*x+4

f´(x) = 3x^2-3
x1=0
f (0)
x2=x1-
f ´ (0)
−(2)2−2(2)−3
x2=2
2(2)−2
x2=1.33
(1.33)2−2(3.5)−3
x3= 1.33−
2(1.33)−2
x3=0.317

Ingrese el valor inicial:

0
Ingrese el porcentaje de error:
0.01
Ingrese la función:
x^3-3*x+4
EJEMPLO 2:
x^2-16*x-16

f´(x) = 2x-16
x1=-7
f (−7)
x2=x1-
f ´(−7)
2
−(−7) −−7(2)−3
x2=2
2(−7)−2
x2= -2.16
(−2.16)2−2(3.5)−3
x3= 3.5−
2(−2.16)−2
x3=-1.017
Ingrese el valor inicial:
-7
Ingrese el porcentaje de error:
0.001
Ingrese la función:
x^2-16*
x-16

5. METODO DEL PUNTO FIJO - PROGRAMACION

%puntofijo('1.14-2*log10*(0.025+3.1167*10^(-4))',1,0.000001)
function puntofijo(g,x0,tol)
g=inline(g);
fprintf('\n it. x g(x) \n')
i=0;
fprintf('%3.0f %10.10f %10.10f \n',i,x0,g(x0))
x1=g(x0);
while abs(x0-x1)>tol
i=i+1;
fprintf('%3.0f %10.10f %10.10f \n',i,x1,g(x1))
x0=x1;
x1=g(x0);
end
fprint('\n La aproximación del punto fijo es %3.10f \n\n',x1)
 14. f ( x )=¿ cos(x)
cos ( x )=x
valor inicial=0

n Columna 1 Columna 2 Columna 3 Columna 4

1 000000000 0000000000 -0.459862 1.000000
2 000000000 0.543020256 0.3179266 0.856662
3 000000000 0.857552366 -0.2032555 0.365259
4 000000000 0.652563265 0.1397156 0.310256
5 000000000 0.792363254 -00625495 0.175426
6 000000000 0.701256632 0.0459222 0.1313132
7 000000000 0.762563645 -0.2892635 0.0821462
8 000000000 0.72655265 0.0192523 0.5642361
9 000000000 0.75124462 -0.0190662 0.0347775
10 000000000 0.739255899 0.0123658 0.2599614
11 000000000 0.745623656 -0.0086595 0.0172655
15.
12 000000000 0.732565562 0.0058962 0.0114269
13 000000000 0.745566263 -0.0039856 0.0075826
14 000000000 0.737532655 0.0025698 0.0053142
15 000000000 0.741325655 -0.0011558 0.0035626

15. f ( x )=x 2−2 x−3

3
g ( x) =
x−2
valor inicial=4
Adicionales
16. f(x)=x2-2x-3 Donde: x0 =2

f´(x) = 2x-2
Primera iteracion:
x1=2
Segunda iteracion:
f (2)
x2=x1-
f ´(2)
−(2)2−2(2)−3
x2=2
2(2)−2
x2= 3.5
Tercera iteracion:
( 3.5)2−2(3.5)−3
x3= 3.5−
2(3.5)−2
x3=3.05

n x F(x) Fp(x)
1 3.50000 2.25000000 5.0000000
2 3.05000 0.20250000 4.1000000
3 3.00060 0.00243956 4.0012195
4 3.00000 0.00000037 4.0000000

17. f(x)=sin(x)+3cos(x) Donde: x0 =1

f´(x) = cos(x)-3 sin(x)

Primera iteracionx1=1
Segunda iteracion
f (1)
x2=x1-
f ´ (1)
−sin ( 1 ) +3 cos ⁡( 1)
x2=1
cos ( 1 ) −3 sin ⁡( 1)
x2= 2.24104868
Tercera iteracion
−sin ( 2.24 )+ 3 cos ⁡( 2.24)
x3=2.24
cos ( 2.24 )−3 sin ⁡( 2.24)
x3=1.877855073

n x F(x) Fp(x)
1 2.2410485 -1.07988562 -2.97217856
2 1.8777175 0.046894585 -3.15192996
3 1.8925547 -0.00003438 -3.16227755
4 1.8954588 0.00000000 -3.16225488

5: f ( x )=−0.5 x 2+2.5 x +4.5 :

PRIMERA INTERACION
5+10
x r=
2
x r=7.5
f ( x 1 )=−0.5 ( 7.5 )2 +2.5 ( 7.5 ) +4.5
f ( x 1 )=−4.875
SEGUNDA INTERACION
5+7.5
x r=
2
x r=6.25
2
f ( x 1 )=−0.5 ( 6.25 ) +2.5(6.25)+ 4.5
f ( x 1 )=0.59375
ERROR

Ea = | 0.59375−(−4.875)
0.59375 | 100%= 7.21%

k a b c F( c) error
1 5.0000 10.0000 7.5000 -4.8750 2.5000
2 5.0000 7.5000 6.2500 0.5938 1.2500
3 6.2500 7.5000 6.8750 -1.9453 0.6250
4 6.2500 6.8750 6.5625 -0.6270 0.3125
5 6.2500 6.5625 6.4063 -0.0044 0.1563
6 6.2500 6.4063 6.3281 0.1474 0.0781
7 6.3281 6.4063 6.3672 0.0717 0.0391
8 6.3672 6.4063 6.3857 0.0337 0.0195
9 6.3867 6.4063 6.3965 0.0335 0.0098
10 6.3956 6.4063 6.4014 0.0147 0.0049

Análisis Numérico Básico de un libro- Luis Rodriguez Ojeda

Pag 32
32: Un círculo de centro O tiene un radio de 50 cm. Dibujar desde el punto P a la
circunferencia exterior de la tangente PA y que une el punto O con el punto de tangencia
A y el punto P; se obtiene un triángulo APO. Sabiendo que el segmento PA es de 120
cm, calcular el área y el perímetro del triángulo. Formule el modelo matemático (una
ecuación no lineal con una incógnita: altura de la perpendicular del centro a cada línea
de corte) y obtenga el ancho de cada uno de los tres cortes.

Datos
Polígono: triángulo AOP
cateto b = PA = 120 cm
cateto c = AO = 50 cm
Solución
Se requiere la hipotenusa de un triángulo rectángulo AOP teniendo:
cateto b = PA = radio AO = 120 cm
cateto c = AO = segmento PA = 50 cm
aplico la fórmula del teorema de Pitágoras:
hipotenusa a = OP = V b² + c² y me pongo:
a = V(120 cm)² + (50 cm)² = 130 cm

Respuesta
La hipotenusa del triángulo rectángulo AOP, que coincide con segmento PO, es 130 cm.
Datos:
Polígono: triángulo POA
PO = 120 cm
OA = 50 cm
AP = 130 cm

Solución
Se requiere el perímetro del triángulo POA cuyos lados son:
PO = 120 cm
OA = 50 cm
AP = 130 cm
aplico la fórmula:
p = PO + OA + AP
y obtengo:
p = 120 cm + 50 cm + 130 cm = 300 cm
Respuesta
El perímetro del triángulo POA es 300 cm

Datos:
Polígono: triángulo rectángulo AOP
base b = AO = 50 cm
altura h = PH = 120 cm

Solución
Calculo el área del triángulo rectángulo AOP que ha:
base b = AO = 50 cm
altura h = PH = 120 cm

aplico la fórmula A = y obtengo:

A= = 3000 cm²

Métodos numéricos – Ricardo Javier Hancco Ancori

Pag 43 – 45
17: Escribir una función CtaDig(num), que reciba un número entero positivo y retorne
cuantos dígitos tiene, así CtaDig(416) debe retornar 3, pues el número 416 tiene 3
dígitos que son 4, 1, 6.
proceso conteoDigitos
definir n,c como entero;
escribir "digite n";
leer n;
c<-0;
mientras n<>0 hacer
n<-trunc(n/10);
c<-c+1;
escribir "despues de la interacion",c,"el numero quedo
en",n;
fin
escribir "el numero tiene",c,"digitos";
fin
ejecución iniciada
digite n
>>416
El numero tiene 3 digitos

22: Escribir una función sumapar(n) que sume todos los números pares desde 2 hasta
n; así sumapar(9), sumará 2+4+6+8 y retornará 20. a) Hágalo usando un for.
b) Luego haga otra versión usando el while.

s=0;
for i=(2+4+6+8)
s=s+i;
end
s
s = 20
function suma=sumapar(n)
suma=0;
i=1;
while i<=n
suma=suma+2*i-1;
i=i+1;
end
>> sumapar(9)
ans = 20
27: Escribir una función que calcule el promedio de todos los números impares que se
encuentren en un vector de números positivos así, si el vector es >>x =[4 5 7 9 3 11 12
15 7 24 1 3 ] ; la función aplicada a este vector deberá devolver 8:75 que es igual a:
(5+7+9+3+11+15+7+13)=8 = 8:75.

s=0;
for i=(5+7+9+3+11+15+7+13)/8
s=s+i;
end
s
>> sumaimpares
s=
8.7500

Pag 19 – 24

2: De un ejemplo de una función de una variable de uso en física, química o ingeniería

(de una aceptable complejidad) y luego escriba las órdenes en Matlab que permitan
trazar su gráfica. Describir la relación de las variables que intervienen e interpretar la
gráfica

√ π 2−9 x 2
f(x) = 2 x
cos ( )
3

clc,clear
x=-1.05:0.01:1.05
y=sqrt(pi.^2-9*x.^2)/cos.^2(x/3)
plot(x,y)
grid on
hold on
17: Escribir los comandos en Matlab que permitan calcular el valor deW para los
siguientes valores de X, Y, Z donde Z es la medida en grados sexagesimales de un
ángulo; luego mostrar los resultados.

f ( c0 ) =1 x 106−1 x 105 (4)−0.25(1 x 106 ) √ 4

f ( c0 ) =100000
c 0 +δ c0 =4 +0.5 ( 4 )=6
f ( c0 + δ c 0 )=1 x 106 −1 x 105 (6)−0.25(1 x 106 ) √ 6
f ( c +δc )=−212372,44
δ c 0 f (c 0)
c 1=c 0−
f ( c 0+ δ c 0 ) −f ( c 0)
( 0.5)( 4)(100000)
c 1=4−
−212372,44−100000
c 1=4+0,64
3
c 1=4,64 g /m

IAN CESPEDES CAMPOS

Problemas creados

METODO DE LA BISECCION
PROBLEMA 1:
biseccion(inline('4*x^2-5*x'),1,1.6,0.01)
k a b c
1 1.0000 1.3000 1.3000
2 1.1500 1.3000 1.1500
3 1.2250 1.3000 1.2250
4 1.2250 1.2625 1.2625
5 1.2438 1.2625 1.2438
6 1.2438 1.2531 1.2531
7 1.2484 1.2531 1.2484

ans =
7.0000 1.2484 -0.0078
PROBLEMA 2:
biseccion(inline('x^2-2*x'),1,2.45,0.0001)
k a b c
1 1.7250 2.4500 1.7250
2 1.7250 2.0875 2.0875
3 1.9063 2.0875 1.9063
4 1.9969 2.0875 1.9969
5 1.9969 2.0422 2.0422
6 1.9969 2.0195 2.0195
7 1.9969 2.0082 2.0082
8 1.9969 2.0025 2.0025
9 1.9997 2.0025 1.9997
10 1.9997 2.0011 2.0011
11 1.9997 2.0004 2.0004
12 1.9997 2.0001 2.0001
13 1.9999 2.0001 1.9999
14 2.0000 2.0001 2.0000

ans =

14.0000 2.0000 -0.0001

METODO DE LA SECANTE

PROBLEMA 1:
ingrese función:'x*log(x-10)'
ingrese primer valor:8
ingrese segundo valor:9
ingrese tolerancia:0.0001
k a b tol
0 8.0000 9.0000
1 8.0000 9.0000 21.7372
2 9.0000 6.8132 62.9432
3 6.8132 5.6274 15.2105
4 5.6274 3.8606 16.7834
5 3.8606 4.8437 14.4546
6 4.8437 4.0514 14.5938
7 4.0514 4.0949 13.9895
8 4.0949 3.8568 13.8475
9 3.8568 3.8478 13.6477
10 3.8478 3.7719 2.7735
11 3.7719 -0.3606 0.3621
12 -0.3606 -0.0372 0.6602
13 -0.0372 0.0476 0.0016
14 0.0476 0.0003 0.0026
ANS:
15 0.0003 0.0007 0.0000

PROBLEMA 2:
ingrese funcion:'x^3+2*x^2+10*x-20'
ingrese primer valor:0
ingrese segundo valor:1
ingrese tolerancia:0.001
k a b tol
1.0000 3.7597
2 1.0000 1.5385 0.3881
3 1.5385 1.3503 0.0188
4 1.3503 1.3679 0.0001
ANS:
5 1.3679 1.3688 0.0000

METODO DE LA FALSA POSICION

PROBLEMA 1:

f_falsapos('x.^3+2*x.^2+10*x-20',1,2,0.0001)
N a c b error

0 1.0000 1.3043 2.0000 100.0000

1.0000 1.3043 1.3579 2.0000 1.3348

2.0000 1.3579 1.3670 2.0000 0.2291

3.0000 1.3670 1.3685 2.0000 0.0386

4.0000 1.3685 1.3688 2.0000 0.0065

5.0000 1.3688 1.3688 2.0000 0.0011

 6.0000 1.3688 1.3688 2.0000 0.0002
Ans
1.3688

PROBLEMA 2:
f_falsapos(‘exp^(-x)-x’,0,1,0.0001)
k a b
0 0 1
1 0 0.61270
2 0 0.57218
3 0 0.56770
4 0 0.56721
Ans
0.56715

NEWTON RAPHRON
PROBLEMA 1:
>> newton_r
Ingrese el valor inicial: 1
Ingrese el porcentaje de error: 0.0001
Ingrese la función: 'exp^x+x^3+2*x^2+10*X-20'
I error
1 0.217145
1.203263 0.013882
1.203200 0.000063
1.203200 10^-9

PROBLEMA 2:
>> newton_r
Ingrese el valor inicial: 0.6
Ingrese el porcentaje de error: 0.0001
Ingrese la función: ‘ln(x^2+1)-exp^(x/2)*cos(pi*x)’
I error
0.6 0.14141
0.45256 0.00603
0.45253 2.65*10^-5
0.45253 0

METODO PUNTO FIJO

PROBLEMA 1:

ingrese funcion f:'cos(x)-3*x'

ingrese funcion g:'cos(x)/3'
ingrese primer valor:pi/8
ingrese tolerancia:0.00001
n xD error
0.0000 0.3927 ------
1.0000 0.3080 0.0291

2.0000 0.3177 0.0030

3.0000 0.3167 0.0003

4.0000 0.3168 0.0000

5.0000 0.3167 0.0000

PROBLEMA 2:
ingrese funcion f:'x^3-x-1'
ingrese funcion g:'x^2+2'
ingrese primer valor:0
ingrese tolerancia:0.00001
n xD error
0.0000 0.0000 ------
1 2 5

2 6 209

3 38 54833

1.0e+09 *

0.0000 0.0000 3.0235

1.0e+18 *

0.0000 0.0000 9.1414

1.0e+37 *

0.0000 0.0000 8.3565

1.0e+75 *

0.0000 0.0000 6.9830

1.0e+151 *

0.0000 0.0000 4.8763

1.0e+303 *

0.0000 0.0000 2.3778

1.0e+202 *

0.0000 1.7815 Inf

11 Inf NaN
Análisis Numérico Básico de un libro- Luis Rodriguez Ojeda

5.7 Determine la raíz real de f(x) = (0.8 – 0.3x)/x:

a) Analíticamente
b) Gráficamente
c) Empleando tres iteraciones en el método de la falsa posición, con valores
iniciales de 1 a 3. Calcule el error aproximado ea y el error verdadero et en
cada iteración.
Clearall
formatshort;
fun=input('Introduzcal a funcion f(x)=','s');
a=input('Introduzca el intervalo menor: ');
b=input('Introduzca el intervalo maor: ');
cont=input('Introduzca el n!mero de iteraciones : ');
f=inline(fun);
for "=1:cont
c=(a#b)$%;
e=abs((b&a)$%);
(",:)=" a b c f(c) e;
if f(a)*f(c)+
b=c;
else
a=c;
end
end
fprintf('-n -t" -ta -tb -tc -tf(c) -terror -n')
disp()
fprintf('oluci/n:-n c=0.2f-n',c)
fprintf('f(c)=0.2f-n',f(c))
fprintf('error=0.2f-n',e)
ezplot(f);
gridon;a)
5.8 Calcule la raíz cuadrada positiva de 18 usando el método de la falsa
posición con es = 0.5%. Emplee como valores iniciales xl = 4 y xu = 5.
>> B = sqrt(18)

B=

4.2426

5.9 Encuentre la raíz positiva más pequeña de la función (x está en radianes)

x2 | cos x | = 5 usando el método de la falsa posición. Para localizar el intervalo
en donde se encuentra la raíz, grafique primero esta función para valores de x
entre 0 y 5. Realice el cálculo hasta que ea sea menor que es = 1%.
Compruebe su respuesta final sustituyéndola en la función original.
COMANDO
%funcion sobre metodo de la falsa posicion.
function y=f_falsapos(h,a,b,tol)
%h es la funcion como cadena.
%[a,b]son los limites del intervalo.
%tol es la tolerancia deseada.
f=inline(h);
c==(a*f(b)-b*f(a))/(f(b)-f(a));n=0;error=100;
fprintf('\t n\t\ta \t\t c \t\t b \t\terror\n')
%adicional, para graficar el intervalo:
x=a-5:0.02:b+5;
y=f(x);
z=zeros(size(x));
plot(x,y)
grid on,hold on
plot(x,z,'r')
while(error>tol)
c=(a*f(b)-b*f(a))/(f(b)-f(a));
disp([n,a,c,b,error])
if(f(a)*f(c)<0)
b=c;
else a=c;
end
error=abs(f(c));
n=n+1;

f_falsapos(‘x^2|cos x|= 5’,0,1,0.001)

Ingrese la funcion:
x^2|cos x|= 5
Ingrese a :
0
Ingrese c :
1
Ingrese el error :
0.001

-) Primera iteración

−1−8−2
0 −1 )−1(−2)
(
x 1= ¿=1
¿
-) Primera iteración

1 (−0.22 )−1.33(−1)
x 2= =1.33
−0.22−(−1)

A B C F(a) F(b) F(c) |c-a|

0.0000 1.0000 2.0000 -2.0000 -1.0000 2.0000 2.0000

1.0000 1.3333 2.0000 -1.0000 -0.2222 2.0000 1.0000

1.3333 1.4000 2.0000 -0.2222 -0.0400 2.0000 0.6667

1.4000 1.4118 2.0000 -0.0400 -0.0069 2.0000 0.6000

1.4118 1.4138 2.0000 -0.0069 -0.0012 2.0000 0.5882

1.4138 1.4141 2.0000 -0.0012 -0.0002 2.0000 0.5862

1.4141 1.4142 2.0000 -0.0002 -0.0000 2.0000 0.5859

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 -0.0000 2.0000 0.5858

1.4142 1.4142 2.0000 -0.0000 0.0000 2.0000 0.5858

El resultado sera 1.4142

5.10 Encuentre la raíz positiva de f(x) = x4 – 8x3 – 35x2 + 450x –1001,

utilizando el método de la falsa posición. Tome como valores iniciales a xl = 4.5
y xu = 6, y ejecute cinco iteraciones. Calcule los errores tanto aproximado
como verdadero, con base en el hecho de que la raíz es 5.60979. Emplee una
gráfica para explicar sus resultados y hacer el cálculo dentro de un es = 1.0%.
%funcion sobre metodo de la falsa posicion.
function y=f_falsapos(h,a,b,tol)
%h es la funcion como cadena.
%[a,b]son los limites del intervalo.
%tol es la tolerancia deseada.
f=inline(h);
c==(a*f(b)-b*f(a))/(f(b)-f(a));n=0;error=100;
fprintf('\t n\t\ta \t\t c \t\t b \t\terror\n')
%adicional, para graficar el intervalo:
x=a-5:0.02:b+5;
y=f(x);
z=zeros(size(x));
plot(x,y)
grid on,hold on
plot(x,z,'r')
while(error>tol)
c=(a*f(b)-b*f(a))/(f(b)-f(a));
disp([n,a,c,b,error])
if(f(a)*f(c)<0)
b=c;
else a=c;
end
 error=abs(f(c));
n=n+1;
end

Ingrese la funcion:
X^4 – 8*x^3 – 35*x^2 + 450*x –1001
Ingrese a :
4.5
Ingrese c :
6
Ingrese el error :
0.001

5.11 Determine la raíz real de x3.5 = 80:

a) En forma analítica.
b) Con el método de la falsa posición dentro de es = 2.5%. Haga elecciones
iniciales de 2.0 a 5.0.
COMANDO FALSA POSICION
%funcion sobre metodo de la falsa posicion.
function y=f_falsapos(h,a,b,tol)
%h es la funcion como cadena.
%[a,b]son los limites del intervalo.
%tol es la tolerancia deseada.
f=inline(h);
c==(a*f(b)-b*f(a))/(f(b)-f(a));n=0;error=100;
fprintf('\t n\t\ta \t\t c \t\t b \t\terror\n')
%adicional, para graficar el intervalo:
x=a-5:0.02:b+5;
y=f(x);
z=zeros(size(x));
plot(x,y)
grid on,hold on
plot(x,z,'r')
while(error>tol)
c=(a*f(b)-b*f(a))/(f(b)-f(a));
disp([n,a,c,b,error])
if(f(a)*f(c)<0)
b=c;
else a=c;
end
error=abs(f(c));
n=n+1;
end

f(x)= 2x3+x-1 Intervalo (0,1)

f(a)= 2(0)3+(0)-1=0
f(b)= 2(1)3+(1)-1=2
0 ( 2 )−1(0)
x1=
2−(0)
x1=0
f(x1) =2(0)3+(0)-1
f(x1) =-1
0 ( 3 )−1(−1)
x2=
2−(−1)
x2=0.3333
f(x2) = 2(0.33)3+(0.33)-1
f(x2) =-0.592612

5.12 Dada f(x) = –2x6 – 1.5x4 + 10x + 2 Use el método de la bisección para
determinar el máximo de esta función. Haga elecciones iniciales de xl = 0 y xu
= 1, y realice iteraciones hasta que el error relativo aproximado sea menor que
5%.
COMANDO BISECCION
xai=input('Ingrese el intervalo inferior: ');
xbi=input('Ingrese el intervalo superior: ');
tol=input('Ingrese el porcentaje de error: ');
syms x;
f=input('Ingrese la funciòn: ');
i=1;

f1=subs(f,x,xai);
f2=subs(f,x,xbi);

ea(i)=100;
if f1*f2 < 0
xa(i)=xai; f1=subs(f,x,xa(i));
xb(i)=xbi; f2=subs(f,x,xb(i));
xr(i)=(xa(i)+xb(i))/2; f3=subs(f,x,xr(i));

fprintf('It. Xa Xr Xb Error aprox \n');

fprintf('%2d \t %11.7f \t %11.7f \t %11.7f \n',i,xa(i),xr(i),xb(i));
while abs(ea(i)) >= tol,
if f1*f3<0
xa(i+1)=xa(i);f1=subs(f,x,xa(i+1));
xb(i+1)=xr(i);f2=subs(f,x,xb(i+1));
end
if f1*f3> 0
xa(i+1)=xr(i);f1=subs(f,x,xa(i+1));
 xb(i+1)=xb(i);f2=subs(f,x,xb(i+1));
end
xr(i+1)=(xa(i+1)+xb(i+1))/2; f3=subs(f,x,xr(i+1));
ea(i+1)=abs((xr(i+1)-xr(i))/(xr(i+1))*100);
fprintf('%2d \t %11.7f \t %11.7f \t %11.7f \t %7.3f \n',...
i+1,xa(i+1),xr(i+1),xb(i+1),ea(i+1));
i=i+1;
end
else
fprintf('No existe una raíz en ese intervalo');
end

1 Metodos Numéricos RICARDO JAVIER HANCCO ANCORI-3

problemas
17.Escribir una función CtaDig(num), que reciba un número entero positivo y retorne cuantos
dígitos tiene, así CtaDig(416) debe retornar 3, pues el número 416 tiene 3 dígitos que son 4, 1,
6.

proceso conteoDigitos
definir n,c como entero;
escribir "digite n";
leer n;
c<-0;
mientras n<>0 hacer
n<-trunc(n/10);
c<-c+1;
escribir "despues de la interacion",c,"el numero quedo
en",n;
fin
escribir "el numero tiene",c,"digitos";
fin

ejecución iniciada
digite n
>>416
El numero tiene 3 digitos

22. Escribir una función sumapar(n) que sume todos los números pares desde 2 hasta n, así
sumapar(9), sumará 2+4+6+8 y retornará 20. a) Hágalo usando un for. b) Luego haga otra
versión usando el while
sumapares
s=0;
for i=(2+4+6+8)
s=s+i;
end
s

s = 20

function suma=sumapar(n)
suma=0;
i=1;
while i<=n
suma=suma+2*i-1;
i=i+1;
end
>> sumapar(9)

ans = 20

27. Escribir una función que calcule el promedio de todos los números impares que se
encuentren en un vector de números positivos así, si el vector es >>x = [ 4 5 7 9 3 11 12 15 7 24
1 3 ] ; la función aplicada a este vector deberá devolver 8.75 que es igual a:
(5+7+9+3+11+15+7+13)/8 = 8.75.
s=0;
for i=(5+7+9+3+11+15+7+13)/8
s=s+i;
end
s

>> sumaimpares

s=

8.7500
Alumno: Leon Garay Renato Fabrizio

NEWTON RAPHSON:

EJERCICIO 1:

f(x)=x3-x-1 Donde: x0 =1

f´(x) = 3x2-1

x1=1

f (1)
x2=x1-
f ´ (1)

−(1)3−(1)−1
x2=1
3(1)2−1
x2= 1.5

(1.5)3−(1.5)−1
x3= 1.5−
3 (1.5)2 −1
x3=1.34783

i Xi f(xi) f´(xi) |xi-xi-1|

0 1 -1 2 -
1 1.5 0.875000 5.75 0.5
2 1.34783 0.100682 4.449937127 0.15217
3 1.32520 0.002058 4.26846512 0.02263
4 1.32472 0.000001 4.264649235 0.00048
EJERCICIO 2:

f(x)=cos(x)-x Donde : x0=1

f´(x)=-1-sen(x)

x1=1

f (1)
x2=x1-
f ´ (1)
−cos ( 1 )−1
x2=1
−1−sen(1)
x2=0.7503638679

−cos ( 0.7503638679 )−1

x3=1
−1−sen( 0.7503638679)
x3=0.73911568745

i Xi f(xi) f´(xi) |xi-xi-1|

0 1 -0.45969 -1.841470985 -
1 0.75036 0.018923 -0.3180978762 0.24964
2 0.73911 0.000046 -0.3263695924 0.01125
3 0.73909 0.000000 -0.3263843739 0.00002
MÉTODO DE LA BISECCIÓN:

EJERCICIO 3:

f(x)=x3+4x2-10 Intervalo:(1,2)

(2−1)
x1=1+
2
x1=1.5

f(x1)=1.53+4(1.5)2-10

f(x1)=2.375

(1.5−1)
x2=1+
1.5
x2=1.25

f(x2)=1.253+4(1.25)2-10

f(x2)=-1.796875

I A C B Error
0 1 1.5 2 0.5
1 1 1.25 1.5 0.25
2 1.25 1.375 1.5 0.1250
3 1.25 1.3125 1.3750 0.0625
4 1.3125 1.3438 1.3750 0.0313
5 1.3438 1.3594 1.3750 0.0156
6 1.3594 1.3672 1.3750 0.0078

EJERCICIO 4:

f(x)=4x2 -5x Intervalo:(1,1.6)

(1.6−1)
x1=1+
1.6
x1=1.375

f(x1)=4(1.375)2+5(1.375)

f(x1)=2.375

(1.375−1)
x2=1+
1.375
x2=1.2727

f(x2)= 4(1.2727)2+5(1.2727)

f(x2)=-1.796875}

I A C B Error
0 1 1.3 1.6 0.3
1 1 1.15 1.3 0.15
2 1.15 1.225 1.3 0.0750
3 1.2250 1.2625 1.3 0.0375
4 1.2250 1.2438 1.2625 0.0188
5 1.2438 1.2531 1.2625 0.0094
6 1.2438 1.2484 1.2531 0.0047

MÉTODO DE LA FALSA POSICIÓN:

EJERCICIO 5:

f(x)=x2-1 Intervalo (0,2)

f(a)=02-1=1
f(b)=22-1=3
0 ( 3 )−2(−1)
x1=
3−(−1)
x1=0.5
f(x1)=0.52-1
f(x1)=-0.75
0.5 ( 3 )−2(−0.75)
x2=
3−(−0.75)
x2=0.8
f(x2)=0.82-1
f(x2)=-0.36

I xa vc vb Error
0 0 0.5 2 100
1 0.5 0.8 2 0.75
2 0.8 0.9286 2 0.36
3 0.9286 0.9756 2 0.1375
4 0.9756 0.9918 2 0.0482
5 0.9918 0.9973 2 0.0163
6 0.9973 0.9991 2 0.0055

EJERCICIO 6:
cos(x) intervalo (1,3)

f(a)=cos(1)=0.54

f(b)=cos(3)=-0.989

1 (−0.989 )−3 (0.54)

x1=
−0.989−0.54
x1=1.7061

f(x1)=cos(1.7061)

f(x1)=-0.1348

1 (−0.1348 )−1.706 (0.54)

x2=
−0.1348−0.54
x2=1.5650
 f(x2)=cos(1.5650)

f(x2)=-0.0057962

I xa vc vb Error
0 1 1.7061 3 100
1 1 1.5650 1.7061 0.1349
2 1.5650 1.5708 1.7061 0.0058
3 1.5650 1.5708 1.5708 0.0000

METODO DE LA SECANTE

EJERCICIO 7:

f(x)=x2-3x-4
xi=7
xi+1=5
f(x1)=72-3(7)-4=24
f(x2)=52-3(5)-4=6
24 ( 5−7 )
x2=7- =4.33
6−24
f(x1)=72-3(7)-4=24
f(x2)=4.332-3(4.33)-4=1.76
1.76 ( 7−4.33 )
x2=7- =4.12
24−1.76
I Ri Ri +1 Ri+2
0 1 1.5 2
1 1 1.25 1.5
2 1.25 1.375 1.5
3 1.25 1.3125 1.3750
4 1.3125 1.3438 1.3750
5 1.3438 1.3594 1.3750

EJERCICIO 8:

f(x)=x3+2x2+10x-20

xi=0

xi+1=1

f(x1)= 03+2(0)2+10(0)-20=-20

f(x2)= 13+2(1)2+10(1)-20=-7

1−0 (−7 )
x2=1- =1.53846153846
−7−(−20)
I Ri Ri +1 Ri+2
0 0 1 -
1 0 1 1.538461
2 1 1.53846 1.39737
3 1.53846 1.39737 1.40632
4 1.39737 1.40632 1.406515
5 1.40632 1.406515 1.406514

MÉTODO DEL PUNTO FIJO:

EJERCICIO 9:

f(x)= cos(x)-3x

I X0 Error
0 0.3080 -
1 0.3080 0.0847
2 0.3177 0.0097
3 0.3167 0.0010
4 0.3168 0.0001
5 0.3167 0.0000
EJERCICIO 10:

F(X) =cos(x) -22x

I X0 Error
0 0.3927 -
1 0.0420 0.3507
2 0.0454 0.0034
3 0.0454 0.000
4 0.0454 0.000
EJERCICIO 11:
Problema 6.15
El balance de masa de un contaminante en un lago bien mezclado se expresa así:
dc
V =W −Qc−KV √ c
dt
Dados los valores de parámetros V =1 x 106 m3 , Q=1 x 105 m 3 /año y
6 0.5
W =1 x 10 g /año , y K=0.25m /año , use el método de la secante modificado
para resolver para la concentración de estado estable. Emplee un valor inicial
c=4 g/m3 y δ=0.5 . Realice tres iteraciones y determine el error relativo
porcentual después de la tercera iteración.
 dc
=0=f (c 0 )
dt
W −Q c 0−KV √ c0
f ( c0 ) = =0
V
f ( c0 ) =W −Qc−KV √ c0

PRIMERA ITERACION

f ( c0 ) =1 x 106−1 x 105 (4)−0.25(1 x 106 ) √ 4

f ( c0 ) =100000

c 0 +δ c0 =4 +0.5 ( 4 )=6

f ( c0 + δ c 0 )=1 x 106 −1 x 105 (6)−0.25(1 x 106 ) √ 6

f ( c +δc )=−212372,44
δ c 0 f (c 0)
c 1=c 0−
f ( c 0+ δ c 0 ) −f ( c 0)

( 0.5)( 4)(100000)
c 1=4−
−212372,44−100000
c 1=4+0,64
3
c 1=4,64 g /m

SEGUNDA ITERACION
6 5 6
f ( c1 ) =1 x 10 −1 x 10 (4,64)−0.25 (1 x 10 ) √ 4,64

f ( c1 ) =−2516,48

c 1+ δ c 1=4+0.5 ( 4,64 ) =6,96

f ( c1 +δ c1 ) =1 x 106−1 x 105 (6,96)−0.25(1 x 106 ) √6,96

f ( c +δc )=−355545,98
δ c 1 f (c 1 )
c 2=c 1−
f ( c1 +δ c 1) −f (c 1)

(0.5)(4,64)(−2516,48)
c 2=4,64−
−355545,98−(−2516,48)
c 2=4,64−0,02
3
c 2=4,62 g/m
TERCERA ITERACION
6 5 6
f ( c2 ) =1 x 10 −1 x 10 (4,62)−0.25(1 x 10 ) √ 4,62

f ( c2 ) =645,37

c 2+ δ c 2=4,62+0.5 ( 4,62 )=6,93

f ( c2 +δ c 2) =1 x 106−1 x 10 5 (6,93)−0.25(1 x 106 ) √6,93

f ( c2 +δ c 2) =−351122,33

δ c 2 f (c2 )
c 3=c 2−
f ( c 2+ δ c 2 )−f (c 2 )

( 0.5)(4,62)(645,37)
c 32=4,64−
−351122,33−(645,37)
−3
c 3=4,62−4,23 x 10
3
c 3=4,62 g/m

Error de aproximación = | A . Actual− A . Antertior

A . Actual |x 100
Error de aproximación = 0 Porque se trabajó con 2 decimales.
Los resultados obtenidos se mostraran en la siguiente tabla.
i xi xi +δxi f(xi) f(xi +δxi) Xi+1
1 4 6 100000 -212372,44 4,64
2 4,64 6,96 -2516,48 -355545,98 4,62
3 4,62 6,93 645,37 -351122,33 4,62

La solución para resolver la concentración de estado estable mediante el uso del método
de la secante modificado es c 3=4,62 g/m3 con un error porcentual del 0%.

EJERCICIO 12
Problema 6.16
Para el problema 6.15, la raíz puede localizarse con iteración de punto fijo como
WKV √ c 0
C=
Q
ó
 W −Q c0
C=( )=0
kV

De las que solo una convergerá para valores iniciales de 2 < c < 6. Seleccione la que sea
correcta y demuestre por qué siempre lo será.
Sabiendo que:
dc
=0=f (c 0 )
dt
W −Q c 0−KV √ c0
f ( c0 ) = =0
V
f ( c0 ) =W −Qc−KV √ c0

Entonces:
W −Q c 0−KV √ c0
f ( c0 ) = =0
V
Representa a:
W −Q c0
C=( )=0
kV
para valores iniciales de 2 < c < 6 ya que convergerá.

EJERCICIO 13
Problema 6.17
Desarrolle un programa amigable para el usuario para el método de Newton-Raphson,
con base en la figura 6.4 y la sección 6.2.3. Pruébelo por medio de repetir el cálculo del
ejemplo 6.3.
Programación:
clear ,clc
cf=input('ingrese funcion a evaluar: ');
symsx
f=inline(cf);
derivada=diff(cf,x);
df=inline(derivada);
tol = input('ingrese tolerancia: ');
error = 50;
x=input('ingrese un valor inicial: ');
n =0;
disp(' n xi error')
while (error>tol)
fprintf('\t%i\t%3.5f\t%f\n', n, x, error);
n=n+1;
x=x-f(x)/df(x);
error= abs(f(x));
end
Calculos del problema 6.3

EJERCICIO 14
Problema 6.18
Desarrolle un programa amigable para el usuario para el método de la secante, con base
en la figura 6.4 y la sección 6.3.2. Pruébelo con la repetición de los cálculos del ejemplo
6.6.
Programacion:
fprintf('METODO DE LA SECANTE\n'); %titulo
symsf(x)%variable simbolica
f(x)=input('Ingrese la funcion: ');
r0=input('Ingrese el primer aproximado: ');
r1=input('Ingrese el segundo aproximado: ');
tol=input('Ingrese la tolerancia maxima de error: ');
n=input('Ingrese el numeromaximo de iteaciones: ');
i=0; %conteo inicial de iteraciones
error=100; %para iniciar el bucle
ri=r0; x %un cambio de variable
para darle forma
rim1=r1; %ri+1=r1
fprintf('iteraciones\t\tri\t\t\t\t\t\t\tri+1\t\t\t\t\t\tri+2\t\t\t\t\t
\terror\n'); %titulo o cabecera de lo que se imprimira
fprintf('%i\t\t\t\t%4.11f\t\t\t\t
%4.11f\t\t\t\t---------\t\t\t\t\t---------\n',i,ri,rim1); %valores
iniciales
while(error>=tol && i<=n)
i=i+1;
rim2=rim1-(f(rim1)*((rim1-ri))/(f(rim1)-f(ri))); %aplicando la
formula de la secante para la nueva raiz
error=abs(((rim1-rim2)/rim2)*100); %formula para
hallar el error
fprintf('%i\t\t\t\t%4.11f\t\t\t\t%4.11f\t\t\t\t%4.11f\t\t\t\t
%4.11f\n',i,ri,rim1,rim2,error);
ri=rim1;
rim1=rim2;
end
fprintf('La raiz es: %4.11f con un error de %4.11f\n',rim2,error);
Cálculos del ejemplo 6.6

EJERCICIO 15
Problema 6.19
Haga un programa amigable para el usuario para el método de la secante modificado,
con base en la figura 6.4 y la sección 6.3.2. Pruébelo con la repetición del cálculo del
ejemplo 6.8.
Programacion:
fprintf('METODO DE LA SECANTE\n'); %titulo
symsf(x)%variable simbolica
f(x)=input('Ingrese la funcion: ');
r0=input('Ingrese el primer aproximado: ');
r1=input('Ingrese el segundo aproximado: ');
tol=input('Ingrese la tolerancia maxima de error: ');
n=input('Ingrese el numeromaximo de iteaciones: ');
i=0; %conteo inicial de iteraciones
error=100; %para iniciar el bucle
ri=r0; x %un cambio de variable
para darle forma
rim1=r1; %ri+1=r1
fprintf('iteraciones\t\tri\t\t\t\t\t\t\tri+1\t\t\t\t\t\tri+2\t\t\t\t\t
\terror\n'); %titulo o cabecera de lo que se imprimira
fprintf('%i\t\t\t\t%4.11f\t\t\t\t
%4.11f\t\t\t\t---------\t\t\t\t\t---------\n',i,ri,rim1); %valores
iniciales
while(error>=tol && i<=n)
i=i+1;
rim2=rim1-(f(rim1)*((rim1-ri))/(f(rim1)-f(ri))); %aplicando la
formula de la secante para la nueva raiz
error=abs(((rim1-rim2)/rim2)*100); %formula para
hallar el error
fprintf('%i\t\t\t\t%4.11f\t\t\t\t%4.11f\t\t\t\t%4.11f\t\t\t\t
%4.11f\n',i,ri,rim1,rim2,error);
ri=rim1;
rim1=rim2;
end
fprintf('La raiz es: %4.11f con un error de %4.11f\n',rim2,error);

Calculo con el ejercicio 6.8

EJERCICIO 16
Problema 6.20
Desarrolle un programa amigable para el usuario para el método de Newton-Raphson
para dos ecuaciones, con base en la sección 6.5. Pruébelo con la solución del ejemplo
6.10.
Programacion
clear ,clc
cf=input('ingrese funcion a evaluar: ');
symsx
f=inline(cf);
derivada=diff(cf,x);
df=inline(derivada);
tol = input('ingrese tolerancia: ');
error = 50;
x=input('ingrese un valor inicial: ');
n =0;
disp(' n xi error')
while (error>tol)
fprintf('\t%i\t%3.5f\t%f\n', n, x, error);
n=n+1;
x=x-f(x)/df(x);
error= abs(f(x));
end
Calculo en 6.10

EJERCICIO 17

√
2
9−x
27) Dada la función f(x) = 2
x +1
a) Determinar su dominio

√f ( x) f(x) ≥ 0
 √ 9−x 2
2
x +1
=
( x+ 3)(x−3)
x 2 +1

=-3≤x≤3

Notación del intervalo: [-3;3]

b) Escribir los comandos de matlab que permiten graficarlo

clc,clear
x=-3:0.01:3
y=sqrt(9-x.^2/x.^2+1)
plot(x,y)
gridon
holdon

EJERCICIO 18

√ π 2−9 x 2
32) Dada la función f(x) = 2 x
cos ( )
3
a) Escribir los comandos de matlab que permiten graficarlo
clc,clear
x=-1.05:0.01:1.05
y=sqrt(pi.^2-9*x.^2)/cos.^2(x/3)
plot(x,y)
gridon
holdon

b) Hacer un bosquejo de la grafica

EJERCICIO 19

37) Graficar Simultáneamente las siguientes funciones:

 2
x
f(x)= y g(x) = x( x 2+1 ¿
√1−x 2

METODOS NUMERICOS
METODO DE LA BISECCIÓN
EJERCICIO N°1
Aplique Bisección para encontrar soluciones, con un error de 0.01, siendo f(x)=
x3-7x2+14x-6, en a = (0,1); b= (1,3.2); c = (3.2, 4); d= (5,10).
El número de iteraciones se calcula con:

Para esto I =7.78, que debe redondearse a 8 iteraciones:

 Ai = 1 y Bi = 3.2
 Sea i=1, f(A1)=? y f(B1)=?.
 Aquí: f(1)=2 y f(3.2)=-0.112
 Si f(A1)* f(B1) 0. Aquí: f(1)* f(3.2)= 2*(-0.112) = - 0.224
 Si i<o= I, calcule pi usando:
  Determine f(Pi) y Si I=1, calcule f(Ai) y vaya a lo siguiente:
I Ai Bi Pi F(Pi) F(Ai) F(Pi)*F(Ai Error
)
1 1 3.2 2.1 1.791 2 3.582
2 2.1 3.2 2.65 0.55212 1.791 0.988856 0.55
5
3 2.65 3.2 2.925 0.08582 0.5521 0.047387 0.275
8 2
4 2.925 3.2 3.0625 -0.05444 0.0858 -0.00467 0.1375
2
5 2.925 3.0625 2.9937 0.00632 0.0858 0.000543 0.0687
5 7 2 5
6 2.9937 3.0625 3.0281 -0.02652 0.0063 -0.00016 0.0343
5 3 2 7
7 2.9937 3.0281 3.0109 -0.01069 0.0063 -0.00006 0.0171
5 3 4 2 8
8 2.9937 3.0109 3.0023 -0.00233 0.0063 -0.00005 0.0085
5 4 4 2 9

 Raíz = 3.0023

EJERCICIO N°2:
2
Hallar las raíces de la ecuación f(x)= 4 x −5 x en [1; 1,6]
1er Paso:
2do Paso: Aplicando el teorema de bolzano
3er Paso:
Primera interacción:
a=1 fa= f(1)= -1
b=1.6 fb=f(1,6)=2.24
a+b 1+1.6
c1= c1= = 1.3
2 2
Segunda interacción:
a=1 f(1)= -1
c1=1.3 f(1.3)= 0.26
1+1.3
c2= =1.15
2

Tercera interacción:
c2= 1.15 f(1.15)= -0.46
c1= 1.3 f(1.3)= 0.26
1.15+1.3
c3= =1.225
2
Cuarta interacción:
c3= 1.225 f(1.225)= -1.225
c1=1.3 f(1.3)= 0.26
1.225+1.3
c4= = 1.2625
2
Quinta interacción:
c3=1.225 f(1.225)= -1.225
c4=1.2625 f(1.2625)= 0.063125
1.225+1.2625
c5= =1.24375
2
f(1.24375)=-0.03109375
i A c b f(a) f(c) f(b) Error
1 1 1.3 1.6 -1 0.26 2.24 0.26

2 1 1.15 1.3 -1 -0.46 0.26 0.46

3 1.15 1.2225 1.3 -0.46 -1.225 0.26 1.225

4 1.225 1.2625 1.3 -1.225 0.063125 0.26 0.0631225

5 1.225 1.24375 1.2625 -1.225 - 0.06312 0.03109375

0.03109375 5

METODO DE PUNTO FIJO

EJERCICIO N°3:
Hallar el punto fijo de: e x-1 -1.5x = 0
1. F(x) = 0
ex-1 -1.5x = 0

Grafica
3.5
3
2.5
2
1.5
1
0.5
0
-3 -2 -1 -0.5 0 1 2 3 4
-1
2. Tabular y graficar F(x) = 0
x ex-1 -1.5x
-2 3.05
-1 1.64
0 0.37
1 -0.5
2 -0.28
3 2.89

3. Sumar en x en F (x):
F(x) + x =x
G(x) = e x-1 -1.5x + x
G(x) = e x-1 + 0.5x
4. Xo = (Xn + Xn+1)/2
Xo = (0+1)/2 = 0.5
 5. G(x)´ < 1
G(x) ´ = e x-1 + 0.5
G(x) ´ = 0.107 <1
6. Iterar:
X1= g(x) = 0.357
X2 = g (0.357) = 0.347
X3 = g (0.347) = 0.347
EJERCICIO N°4:
Para la ecuación f(x)= x 3 + x =6 considerando [1; 2] punto medio 1.5
Despejando x:
1) f(x)= x 3 + x =6
g(x)= (6-x)1/3
2) f(x)= x 3 + x =6
g(x)= (6- x 3 )
Derivamos:
1) g(x)= (6-x)1/3
1
g’(x)= -
3 (6−x)2/ 3
2) g(x)= (6- x 3 )
g’(x) = -3 x 2

Reemplazamos:
1
g’(x)= -
3 (6−1.5)2 /3

g’(x)= - 0.12293601

Si j=0 X1=g(Xo)= g(1.5)= (6-1.5)1/3 = 1.650969

Si j=1 X2=g(X1)= g(1.560969) =(6-1.650969)1/3= 1.632291
Si j=2 X3=g(X2)= g(1.632291)= (6-1.632291)1/3= 1.634624
Si j=3 X4=g(X3)=g(1.634624)= (6-1.634624)1/3 =1.634330
Si j=4 X5=g(X4)=g(1.634330)= (6-1.634330)1/3 =1.636932
Si j=5 X6=g(X5)=g(1.636932)= (6-1.636932)1/3=1.634364
 1.650969−1.5
E1= x100% =7.77%
1.650969

1.632291−1.650969
E2= x100% =1.14%
1.632291

1.634624−1.632291
E3= x100% =0.08%
1.634624

1.634330−1.634624
E4= x100%=0.017%
1.634330

1.636932−1.634330
E5= x100%= 1.79%
1.636932

METODO DE LA SECANTE

EJERCICIO N° 5

Hallar una buena aproximación a la función f(x) = x 3+2x2+10x-20 con xo= 0 y x1

=1

(Xi-1 – X i-2) Para i =2,3,4…

Xi = X i-1 F(Xi-1)
F(Xi-1)–F( X i-2)

 Para i = 2

(X1– X 0)
X2 = X 1 F(X1)
F(X1)–F( X 0)

Antes:

F(X 0) = F(0) = (0)3+2(0)2+10(0)-20 = -20

F(X 1) = F(1) = (1)3+2(1)2+10(1)-20 = -7

 Entonces reemplazamos:

X2 = 1 (1– 0) (-7) = 1.5385

 (-20)– (-
7)

 Evaluando:

F(X i) = F(x2) = F(1.5385) = (1.5385)3+2(1.5385)2+10(1.5385)-20 = 3.749

 Factor:

(1.5385 -1) = 0.538

 Tabla:

i X0 X1 F(X0) F(X1) Xi F(Xi) (Xi - Xi-1)

2 0 1 -20 -7 1.538 3.748 0.536

3 1 1.538 -7 3.758 1.350 -0.384 4.142

4 1.538 1.350 3.756 -0.384 1.367 -0.016 0.403

5 1.350 1.367 -0.384 -0.019 1.357 -0.228 0.248

6 1.367 1.357 -0.019 -0.229 1.372 -0.073 0.155

 I = i +1

i=2+1=3

 Una buena aproximación seria de x =1.3

EJERCICIO N°6

F(x)= x 2 - e−x considerando Xo=-1, X1=2

Primera interacción:
f(Xo)= 1- e−1
f(Xo)=-1.718281
f(X1)= 4- e−2
f(X1)=3.864664
(3.864664 )(2−(−1 ))
X2= 2-
(3.864664 )−(−1.718281)
X2= -0.076680

Segunda interacción:
f(X1)=3.864664
2
f(X2) = −0.076680 −e
−0.076680

f(X2)= -1.073816
(−1.073816)(−0.076680−2)
X3= -0.076680 -
−1.073816−3.864664
X3=0.374872

Tercera interacción:
f(X2)= -1.073816
f(X3)= 0.3748722−e−0.374872
f(X3)= -0.546848
(−0.546848)(0.374872−0.076680)
X4= 0.374872−
−0.546848−(−1.073816)
X4=0.843458
Cuarta interacción:
f(X3)= -0.546848
f(X4)= 0.8434582−e−0.843458
f(X4)=0.281201
(0.281201)(0.843458−0.374872)
X5= 0.843458 -
0.281201−(−0.546848)
X5=0.684328

Quinta interacción:
f(X4)= 0.281201
f(X5)= 0.6843282−e0.684328
f(X5)= -1.514134
(−1.514134)(0.684328−0.374872)
X6= 0.684328−
−1.514134−( 0.843458)
X6= 0.485583

METODO DE FALSA POSICIÓN

EJERCICIO N° 7

Usar el método de la regla falsa para aproximar la raíz de ,

comenzando en el intervalo y hasta que .

 Sabemos que f(x) es continua en el intervalo dado y que toma signos
opuestos en los extremos de dicho intervalo. Por lo tanto podemos
aplicar el método de la regla falsa.
 Calculamos la primera aproximación:

 Puesto que solamente tenemos una aproximación, debemos seguir con

el proceso.
 Así pues, evaluamos

 Y hacemos
nuestra tabla de signos:

 De donde vemos que la raíz se encuentra en el intervalo .

Calculamos la nueva aproximación:

 En este momento, podemos calcular el primer error aproximado:

 Puesto que no se cumple el objetivo seguimos con el proceso.

Evaluamos , y hacemos la
 tabla de signos:

De donde vemos que la raíz se encuentra en el intervalo

, con el cual, podemos calcular la nueva aproximación:

EJERCICIO N° 8

Usar el método de la falsa posición para f(x)= 2xcosx-3 comenzando en el

intervalo [1;2]
Xa=1 f(1)= 2(1)cos(1)-3= -1.00030
Xb=2 f(2)= 2(2)cos(2)-3= 0.99756
Primera interacción:
f ( xa)(xb−xa)
Xr= Xa -
f ( xb )−f (xa)
(−1.00030)(2−1)
X1= 1-
0.99756−(−1.00030)
X1=1-(-0.50068)
X1=1.50068
f(1.50068)=0.00033
1.50068−1
E% = ¿ | x100%
1.50068
E%=3.3%

Segunda interacción:
Xa= 1.50068
Xb= 2
(0.00033)(2−1.50068)
X2=1.50068-
0.99756−(0.00033)
X2=1.50068-0.00016
X2=1.50052
f(1.50052)=0.00001
1.50052−1.50068
E%= ¿ |x100%
1.50052
E%=0.03%
Tercera interacción:
Xa=1.50052
Xb=1.50068
(0.00001)(1.50068−1.50052)
X3=1.50052 -
0.00033−(0.00001)
X3=1.50052 – 0.000005
X3=1.500515
f(1.500515)=0.0000009
1.50015−1.50052
E%= ¿ ∨¿ x 100%
1.50015
E%=0.026%

Cuarta interacción:
Xa=1.500515
Xb=1.50052
(0.0000009)(1.50052−1.500515)
X4=1.500515 -
0.00001−(0.0000009)
X4= 1.500515 – 0.0000004
X4=1.500514
f(1.500514)= 0.000001
1.500514−1.50015
E%= ¿ ∨¿ x100%
1.500514
E%=0.024%

METODO DE NEWTON RAPSON

EJERCICIO N° 9
Encontrar una buena proximacion a una raíz de las siguientes funciones,
usando el método de Newton Raphson. Tomar como punto de partida a x = 1
F(x) = x3-x-1 = 0

 Sabemos:

F(Xn)
Xn+1 = Xn Si F´(Xn) no es igual a 0
F´(Xn)

 Entonces:
F´ (Xn) = 3x2-1
Se sabe que x1 = 1
 Por ello:

F(X1)
X2 = X 1 -
F´(X1)

(-1)
X2 = 1- = 1.5
(2)

F(X2)
X3= X 2 -
F´(X2)

F (1.5)
X3= 1.5- = 1.34782
F
´(1.5)
 F(X3)
X4= X 3 -
F´(X3)

F (1.34782)
X4= 1.34782 - = 1.32520
F´(1. 34782)

F(X4)
X5= X 4 -
F´(X4)

F (1.32520)
X5= 1.32520- = 1.3247
F´(1.32520)

 Raíz mas aproximada es : 1.32

EJERCICIO N° 10

Resolver considerando como punto inicial X0 =0

Donde f(x)=2x -cosx – 3
Xn=0
f(X)= 2(0) -cos(0)-3
f(X)= -4
f’(X)= 2+ sen(0)
f’(X)= 2

Primer interacción:
f (x 0)
X1=X0- ( )
f ´ ( x 0)
−4
x1= 0- ( )
2
X1=2
f(x1)= 2(2) -cos(2)-3
f(x1)=0.000609
f’(x1)= 2+sen(2)
f’(x1)= 2.034899
2−0
E%= x100
2
E%=1%

Segunda interacción:
0.000609
X2=2- ( )
2.034899
X2=1.999700
f(X2)=2(1.999700)-cos(1.999700)-3
f(X2)=0.00000899
f’(X2)= 1.999700+sen(1.999700)
f´(X2)=2.034594
1.999700−2
E%= X100
1.999700
E%=0.0150%

Tercera interacción:
0.00000899
X3= 1.999700 −( )
2.034594
X3=1.999695
f(X3)=2(1.999695)-cos(1.999695)-3
f(X3)= -0.00000101
f’(X3)= 1.999695+sen(1.999695)
f’(X3)=1.964800
1.999695−1.999700
E%= X100
1.999695
E%=1.00001%
Cuarta interacción:
0.00000101
X4=1.9996956- −( )
1.964800
X4= 1.999695514
f(X4)= 2(1.999695514)-cos(1.999695514)-3
f(X4)= 0.00000005623
f’(X4)= 1.999695514+ sen(1.999695514)
f’(X4)= 2.0345897
1.999695514−1.999695
E%= X 100
1.999695514
E%=2.570%
5.13 La velocidad v de un paracaidista que cae está dada por v=gmc(1-e-
(c/m)t)
donde g=9.8 m/s2. Para un paracaidista con coeficiente de arrastre de c=15
kg/s, calcule la masa m de modo que la velocidad sea v=35 m/s en t=9s. Utilice
el método de la falsa posición para determinar m a un nivel de Es=0.1%.

Raíz aproximada =
59

m f(m)

54 -2.61

56 -1.69

58 -0.80

60 0.068

Raíz aproximada = 60
xl=58 xu=60
xr1=60-(0.068)(58-60)-0.8+0.068=59.81420765
xr2=59.8142-(-0.0115)(58-59.81)-0.8+0.0115=59.84079423
Ea=59.84079423-59.8142076559.84079423(100%)=0.044%

Iteración xl xu xr Ea

1 58 60 59.81420765 -

2 58 59.81420765 59.84104239 0.044%

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.

i xl xu xr a
1 2.5 3 2.75 -----
2 2.5 2.75 2.625 4.761
3 2.5 2.625 2.5625 2.439
4 2.5625 2.625 2.59375 1.204
5 2.59375 2.625 2.609375 0.598
6 2.59375 2.609375 2.6015625 0.300
7 2.3015625 2.609375 2.60546875 0.149
8 2.60546875 2.609375 2.607421875 0.074
9 2.607421875 2.609375 2.608398438 0.0374
1
2.608398438 2.609375 2.608886719 0.018
0
11 2.608398438 2.608886719 2.608642579 0.009

xr=2.5+32=2.75 a=2.625-2.752.625×100= 4.761

xr=2.5+2.752=2.625 a=2.5625-
2.6252.5625×100=2.439
xr=2.5+2.6252=2.5625 a=2.59375-2.56252.59375×100=
1.204

xr=2.5625+2.6252=2.59375 a=2.6093-
2.593752.6093×100=0.598

xr=2.59375+2.6252=2.609375 a=2.6015625-
2.593752.6015625×100= 0.300

xr=2.59375+2.6093752=2.6015625 a=2.60546875-
2.60156252.60546875×100=0.149

xr=2.6015625+2.6093752=2.60546875 a=2.607421875-
2.605468752.607421875×100=0.074

xr=2.60546875+2.6093752=2.607421875 a=2.608398438-
2.6074218752.608398438×100=0.0374

xr=2.607421875+2.6093752=2.608398438 a=2.608886719-
2.6083984382.608886719×100=0.018

xr=2.608398438+2.6093752=2.608886719 a=2.608642579-
2.6088867192.608642579×100=0.009

xr=2.608398438+2.6088867192=2.608642579

5.15 Por un canal trapezoidal fluye agua a una tasa de Q=20m3s2. La

profundidad crítica y para dicho canal satisface la ecuación
0=1-Q2gA3B
Donde g=9.81ms2, A=área de la sección transversalm2, y B=ancho del canal
en la superficiem. Para este caso, el ancho y el área de la sección transversal
se relacionan con la profundidad y por medio de
B=3+y y A=3y+y22
Resuelva para la profundidad crítica con el uso de los métodos a) gráfico, b)
bisección y, c) falsa posición. En los incisos b) y c), haga elecciones iniciales de
Xl=0.5 y Xu=2.5, y ejecute iteraciones hasta que el Ea caiga por debajo de 1%
o el número de iteraciones supere a 10. Analice sus resultados.
a) Gráfico
 fy=40.774-y2+6y38y+24

b) Bisección

N Xl Xu Xr fxf(Xr) Ea %
1 0. 2.5 1.5 + -
5
2 1. 2.5 2 - 33.3%
5
3 1. 2 1.75 - 14.3%
5
4 1. 1.75 1.625 - 7.7%
5
5 1. 1.625 1.5625 - 4%
5
6 1. 1.5625 1.53125 - 2.04%
5
7 1. 1.53125 1.516075 - 1%
5
8 1. 1.516075 1.5080375 + 0.53%
5
5.16

VENTANA DE COMANDOS
>> f=inline('((9*pi*h.^2-pi*h.^3)/3)-30') f = Inlinefunction: f(h) = ((9*pi*h.^2-
pi*h.^3)/3)-30
>> h=0:0.1:15;
>>plot(h,f(h)),grid

Encontramos dos raices una entre 2.02 y 2.04

Y la otra entre 8.6 y 8.65

Para la raíz 1 tenemos a=2.02, b=2.04 y n=8. n a r b

a = 8.61250000000000
>> r=(a+b)/2 r = 8.61875000000000
>> f(a) ans = 0.09946626904214
>> f(b) ans = -0.78686987550585
>> f(r) ans = -0.34301227886772
>> e=(b-a)/2 e = 0.00625000000000
>> b=r b = 8.61875000000000
>> r=(a+b)/2 r = 8.61562500000000
>> f(a) ans =

0.09946626904214
>> f(r) ans = -0.12160071969568
>> f(b) ans = -0.34301227886772
>> e=(b-a)/2 e = 0.00312500000000
>> b=r b = 8.61562500000000
>> r=(a+b)/2 r = 8.61406250000000
>> f(a) ans = 0.09946626904214
>> f(r) ans = -0.01102416600664
>> f(b) ans = -0.12160071969568
>> e=(b-a)/2 e = 0.00156250000000
>> b=r b = 8.61406250000000
>> r=(a+b)/2 r = 8.61328125000000
>> f(a) ans = 0.09946626904214
>> f(r) ans = 0.04423181484981
>> f(b) ans = -0.01102416600664
>> e=(b-a)/2 e = 7.812500000001776e-004
>> a=r a = 8.61328125000000
>> r=(a+b)/2 r = 8.61367187500000
>> f(a) ans = 0.04423181484981
>> f(r) ans = 0.01660651544186
>> f(b) ans = -0.01102416600664
>> e=(b-a)/2 e = 3.906250000005329e-004
>> a=r a = 8.61367187500000
>> r=(a+b)/2 r = 8.61386718750000
>> f(a) ans = 0.01660651544186
>> f(r) ans = 0.00279184749608
>> f(b) ans = -0.01102416600664
>> e=(b-a)/2 e = 1.953124999998224e-004
>> r=(a+b)/2 r = 8.61386718750000
Alumna: Benavente Gomez Melany Helein

METODO DE LA BISECCION
 CODIGO

f(x)= -0.4x2+2.2x+4.7 Intervalo:(5,10)

error: 0.1

(10+5)
x1=
2
x1=7.5

f(x1) = −0.4(7.5)2 +2.2 ( 7.5 )+ 4.7

f(x1) =1.3

(7.5+5)
x2=
2
x2=6.25
 2
f(x2) = −0.4(6.25) +2.2 ( 6.25 )+4.7
f(x1) =2.8

k a b c F(c) Error
1.0000 5.0000 10.0000 7.5000 -1.3000 2.5000
2.0000 5.0000 7.5000 6.2500 2.8250 1.2500
3.0000 6.2500 7.5000 6.8750 0.9188 0.6250
4.0000 6.8750 7.5000 7.1875 -0.1516 0.3125
5.0000 6.8750 7.1875 7.0313 0.3934 0.1563
6.0000 7.0313 7.1875 7.1094 0.1233 0.0781
7.0000 7.1094 7.1875 7.1484 -0.0135 0.0391
8.0000 7.1094 7.1484 7.1289 0.0551 0.0195
9.0000 7.1289 7.1484 7.1387 0.0208 0.0098
10.0000 7.1387 7.1484 7.1436 0.0037 0.0049
11.0000 7.1436 7.1484 7.1460 -0.0049 0.0024
12.0000 7.1436 7.1460 7.1448 -0.0006 0.0012
13.0000 7.1436 7.1448 7.1442 0.0015 0.0006
14.0000 7.1442 7.1448 7.1445 0.0005 0.0003
15.0000 7.1445 7.1448 7.1446 -0.0001 0.0002
16.0000 7.1445 7.1446 7.1445 0.0002 0.0001
17.0000 7.1445 7.1446 7.1446 0.0001 0.0000
18.0000 7.1446 7.1446 7.1446 -0.0000 0.0000
19.0000 7.1446 7.1446 7.1446 0.0000 0.0000
20.0000 7.1446 7.1446 7.1446 -0.0000 0.0000

Grafica:

f(x)= 4x3-1.2x2+7x+3.3 Intervalo:(4,8)

error: 0.2

(8+ 4)
x1=
2
x1=6

f(x1) =4(6)3-1.2(6)2+7(6)+3.3

f(x1) =866.1

(4 +6)
x2=
2
x2=5

f(x2) = 4(5)3-1.2(5)2+7(5)+3.3

f(x2) =508.3

k a b c F(c) Error
0.0010 0.0040 0.0080 0.0060 0.8661 0.0020
0.0020 0.0060 0.0080 0.0070 1.3655 0.0010
0.0030 0.0070 0.0080 0.0075 1.6758 0.0005
0.0040 0.0075 0.0080 0.0077 1.8474 0.0003
0.0050 0.0077 0.0080 0.0079 1.9375 0.0001
0.0060 0.0079 0.0080 0.0079 1.9836 0.0001
0.0070 0.0079 0.0080 0.0080 2.0070 0.0000
0.0080 0.0080 0.0080 0.0080 2.0187 0.0000
0.0090 0.0080 0.0080 0.0080 2.0246 0.0000
0.0100 0.0080 0.0080 0.0080 2.0275 0.0000
0.0110 0.0080 0.0080 0.0080 2.0290 0.0000
0.0120 0.0080 0.0080 0.0080 2.0298 0.0000
0.0130 0.0080 0.0080 0.0080 2.0301 0.0000
0.0140 0.0080 0.0080 0.0080 2.0303 0.0000
0.0150 0.0080 0.0080 0.0080 2.0304 0.0000
0.0160 0.0080 0.0080 0.0080 2.0305 0.0000
0.0170 0.0080 0.0080 0.0080 2.0305 0.0000
0.0180 0.0080 0.0080 0.0080 2.0305 0.0000
0.0190 0.0080 0.0080 0.0080 2.0305 0.0000
0.0200 0.0080 0.0080 0.0080 2.0305 0.0000

Grafica:
 FALSA POSICION
 CODIGO

f(x)= x3-2x-2 Intervalo (1,2)

tolerancia: 0.001

a=1 f(a)= (1)3-2(1)-2=-3

b=2 f(b)= (2)3-2(2)-2=2
 1 ( 2 )−2(−3)
x1=
2−(−3)
x1=1.6
f(x1) =(1.6)3-2(1.6)-2
f(x1) =-1.104
1.6 ( 2 ) −(2)(−1.104)
x2=
2−(−1.104)
x2=1.7423
f(x2) =(1.7423)3-2(1.7423)-2
f(x2) =-0.1959

A B C F(a) F(b) F(c) |c−a|

1.0000 4.0000 2.0000 -3.0000 6.0000 -2.0000 1.0000
1.0000 2.0000 4.0000 -3.0000 -2.0000 6.0000 3.0000
2.0000 2.5000 4.0000 -2.0000 -0.7500 6.0000 2.0000
2.5000 2.6667 4.0000 -0.7500 -0.2222 6.0000 1.5000
2.6667 2.7143 4.0000 -0.2222 -0.0612 6.0000 1.3333
2.7143 2.7273 4.0000 -0.0612 -0.0165 6.0000 1.2857
2.7273 2.7308 4.0000 -0.0165 -0.0044 6.0000 1.2727
2.7308 2.7317 4.0000 -0.0044 -0.0012 6.0000 1.2692
2.7317 2.7320 4.0000 -0.0012 -0.0003 6.0000 1.2683
2.7320 2.7320 4.0000 -0.0003 -0.0001 6.0000 1.2680
2.7320 2.7320 4.0000 -0.0001 -0.0000 6.0000 1.2680
2.7320 2.7320 4.0000 -0.0000 -0.0000 6.0000 1.2680
2.7320 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2680
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679
2.7321 2.7321 4.0000 -0.0000 -0.0000 6.0000 1.2679

Grafica:
 f(x)= x3+2x2+10x-20 Intervalo (1,2)

tolerancia: 0.0001

a=1 f(a)= (1)3+2(1)2+10(1)-20 =-7

b=2 f(b)= (2)3+2(2)2+10(2)-20 =16
1 ( 16 ) −2(−7)
x1=
−7−16
x1=1.3043
f(x1) = (1.3043)3+2(1.3043)2+10(1.3043)-20
f(x1) =1.3348
1.3043 ( 16 ) −2(1.3348)
x2=
16−1.3348
x2=1.3579
f(x2) = (1.3579)3+2(1.3579)2+10(1.3579)-20
f(x2) =0.2291

A B C F(a) F(b) F(c) |c−a|

1.0000 1.3043 2.0000 -7.000 -1.3348 16.0000 1.0000
1.3043 1.3579 2.0000 -1.3348 -0.2291 16.0000 0.6957
1.3579 1.3670 2.0000 -0.2291 -0.0386 16.0000 0.6421
1.3670 1.3685 2.0000 -0.0386 -0.0065 16.0000 0.6330
1.3685 1.3688 2.0000 -0.0065 -0.0011 16.0000 0.6315
1.3688 1.3688 2.0000 -0.0011 -0.0002 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0002 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312
1.3688 1.3688 2.0000 -0.0000 -0.0000 16.0000 0.6312

Grafica:

NEWTON RAPHSON

 CODIGO
f(x)= -0.9x2+1.7x+2.5 Donde: x0 =5

error limite: 0.01

f´(x) = -1.8x+1.7

x1=5

f (2)
x2=x1-
f ´(2)

−−0.9(5)2 +1.7 ( 5 ) +2.5

x2=5
−1.8 ( 5 )+1.7
x2= 3.4

−0.9(3.4 )2+1.7 ( 3.4 ) +2.5

x3= 3.4−
−1.8 (3.4 ) +1.7
x3=2.9

n x F(x) Fp(x)
1 3.4246575342 -2.2335334960 -4.4643835616
2 2.9243569966 -0.2252705651 -3.5638425940
3 2.8611469757 -0.0035959561 -3.4500645562
4 2.8601046891 -0.0000009777 -3.4481884403

Grafica:
f(x)= x2-2x+1 Donde: x0 =0.5

error límite:0.01

f´(x) = 2x-2

x1=0.5

f (1)
x2=x1-
f ´ (1)

−(0.5)2−2 ( 0.5 )+ 1
x2=0.5
2 ( 0.5 )−2
x2= 0.75

−(0.75)2−2 ( 0.75 )+ 1
x3=0.75
2 ( 0.75 )−2
x3=0.875

n x F(x) Fp(x)
1 0.7500000000 0.0625000000 -0.5000000000
2 0.8750000000 0.0156250000 -0.2500000000
3 0.9375000000 0.0039062500 -0.1250000000
4 0.9687500000 0.0009765625 -0.0625000000
5 0.9843750000 0.0002441406 -0.3125000000
6 0.9921875000 0.0000610352 -0.0156250000

Grafica:
 METODO DE LA SECANTE

 CODIGO

f(x)=exp(-x)-x e=0.05

xi=0

xi+1=1

f(x1)= exp(-0)-0=1

f(x2)= exp(-1)-1=-0.6321
 1 (1−0 )
x2=1- =0.613
1−(−0.6321)
f(x2)= exp(-0.613)-0.613= -0.0712

(−0.6321) ( 0.613−1 )
x3= 0.613 - =¿ 0.564
−0.6321−(1)

f(x)= x3+2 x2+10x-20 e=0.000001

xi=0

xi+1=1

f(x0)= (0)3+2 (0)2+10(0)-20= -20

f(x1)= (1)3+2 (1)2+10(1)-20=-7

−7 ( 1−0 )
x2=1- =1.5385
−7−(−20)
f(x2)= (1.5385)3+2 (1.5385)2+10(1.5385)-20= 3.7597

(3.7597) (1.5385−1 )
x3= 1.5385 - =0.3881
3.7597−(−7)
METODO DE LA SECANTE

f(x)= cos(x)-3*x

g(x)=cos(x)/3

primer valor: pi/8

n X0 error
0.0000 0.3927 --------
1.0000 0.3080 0.0847
2.0000 0.3177 0.0097
3.0000 0.3167 0.0010
4.0000 0.3168 0.0001
5.0000 0.3167 0.0000
6.0000 0.3168 0.0000

LIBRO RICARDO HANCO ANCORI (PAG.19-24)

12. Dada la expresión:

a) Escribir todos los comandos Matlab necesarios para calcular W si x = 30◦ , y = e, z = 4/7
(x está en grados sexagesimales, e es el número de Euler o base del logaritmo natural).

b) Mostrar el resultado.

Solución:

 x= [ 30 ]
 y= [ exp ]
 z= [ 4 /7 ]
( x ) + √ ( 1+ z )
 w=
(( ( y .3 ) /( y + z .2)−log ⁡( y ).2). /( x .∗sin ¿ ./(z + y ))) ¿
W=

17. Escribir los comandos en Matlab que permitan calcular el valor de W para los siguientes
valores de X, Y, Z donde Z es la medida en grados sexagesimales de un ángulo; luego
mostrar los resultados.

Solución:

 x= [ √(3)ex p(2) pi ]
 y= [ 2.24 3.56−4.67 ]
 z= [ 30 45 60 ]
( z)
 w= log ⁡( x .2+ y ).2
z .∗cos ¿ ./¿
¿
W=

Escribir los comandos necesarios para graficar la función y=√ 4−| x|

Solución:

En primer lugar, resolveremos la inecuación 4− |x| ≥ 0 para determinar el dominio de la

función pedida. Luego, los puntos de referencia son: x = − 4 ∨ x =4 , por tanto x ∈ − 4,4.

 x= -4:0.01:4
 y=sqrt(4-abs(x)

 plot(x,y),grid

LIBRO CHAPRA METODOS NUMERICOS

5.19. Desarrolle un subprograma para el método de bisección que minimice las

evaluaciones de la función, con base en el seudocódigo que se presenta en la figura 5.11.
Determine el número de evaluaciones de la función (n) para el total de iteraciones. Pruebe
el programa con la repetición del ejemplo 5.6.

5.20. Desarrolle un programa amigable para el usuario para el método de la falsa posición.
La estructura del programa debe ser similar al algoritmo de la bisección que se bosquejó en
la figura 5.10. Pruebe el programa con la repetición del ejemplo 5.5.
5.21. Desarrolle un subprograma para el método de la falsa posición que minimice las
evaluaciones de la función en forma similar a la figura 5.11. Determine el número de
evaluaciones de la función (n) para el total de iteraciones. Pruebe el programa por medio
de la duplicación del ejemplo 5.6.

5.22. Desarrolle un subprograma amigable para el usuario para el método de la falsa

posición modificado, con base en la figura 5.15. Pruebe el programa con la determinación
de la raíz de la función del ejemplo 5.6. Ejecute corridas hasta que el error relativo
porcentual verdadero esté por debajo de 0.01%. Elabore una gráfica en papel
semilogarítmico de los errores relativo, porcentual, aproximado y verdadero, versus el
número de iteraciones. Interprete los resultados.
 6.1. Utilice la iteración simple de punto fijo para localizar la raíz de

x
√¿
¿
f ( x)=2 sen ¿
Haga una elección inicial de x0 = 0.5 e itere hasta que ea ≤ 0.001%.

Compruebe que el proceso converge en forma lineal según se describió en el recuadro 6.1.

Raíz Error aproximado

0.7993 37.4432
0.7699 5.1788
0.7710 1.4336
0.7680 0.3846
0.7655 0.1042
0.7686 0.0282
0.7607 0.0076
0.7686 0.0022
0.7656 0.0006

Ejercicio 19

6.2. Determine la raíz real más grande de

f(x) = 2x3 – 11.7x2 + 17.7x – 5

a) En forma gráfica.
b) Con el método de Newton-Raphson (tres iteraciones, x0 = 3, d = 0.001).

f(x)=2 x3-11.7 x2+17.7x-5

f´(x0)= 6x2-23.4x+17.7

x0=3

f (x 0)
x1=x0-
f ´ ( x 0)
2( 3)3−11.7(3) 2+17.7(3)−5
x1=3-
6( 3)2−23.4 (3)+17.7
x1=5.13

f (x 1)
x2=x1-
f ´ ( x 1)
2( 5.13) 3−11.7 (5.13)2+ 17.7(5.13)−5
x2=5.13-
6( 5.1)2−23.4 (5.13)+17.7
x2=4.2678

f (x 2)
x3=x2-
f ´ ( x 2)
2( 4.2)3−11.7 (4.2)2+17.7 (4.2)−5
x3=4.2-
6(4.2)2−23.4( 4.2)+17.7
x3=11.91

c) Con el método de la secante (tres iteraciones x–1 = 3, x0 = 4).

f(x)=2 x3-11.7 x2+17.7x-5

xi=3

xi+1=4

f(x0)= 2 (3)3-11.7 (3)2+17.7(3)-5= -3.2

f(x1)= 2 (4)3-11.7 (4)2+17.7(4)-5= 6.6

 6.6 ( 3−4 )
x2=4-
−3.2−6.6
x2=3.3266

f(x2)= 2 (3.3)3-11.7 (3.3)2+17.7(3.3)-5= -2.19

−2.19 ( 4−3.3 )
x3=3.3-
6.6−(−2.19)
x3=3.5

f(x3)= 2 (3.5)3-11.7 (3.5)2+17.7(3.5)-5= -0.6

−(−0.6) ( 3.3−3.5 )
x3=3.3-
−2.19−(−0.6)
x3=3.5
 Programación metod numéricos …ejercicios
creados
Alumno : pedro bedoya cutipa
Bisección

clear all % se borra texto del comando pront

format short ;
a=input('introdusca el valor de a: ') ; % se asigna el intervlo
inferior
b=input('introdusca el valor de b: ') ; % se asigna el intervlo
superior
cont=input('introdusca el numero de iteraciones cont:'); % se asigna
el numero de iteraciones
fun=input('introdusca la funcion f(x)=','s'); % se asigna la funcion
f=inline(fun);
for k=1: cont% aqui es donde se realiza la iteracion
c=(a+b)/2% aqui es donde se calcula el valor del punto medio
e=abs((b-a)/2);% aqui se calculael error absoluto estimado
A(k,:)=[k a b c f(c) e];% aqui se imprime los valores de iteracion
if f(a)*f(b)<0
b=c;% si es menor a 0 se encuentra en el intervalos;
else
a=c;% si es mayor que 0 se encontro la raiz;
end
end
fprintf('\n\tk\ta\tb\tc\tf(c)\terror\n')
disp(A)
fprintf('solucion:\n c=%8.5f\n',c);
fprintf('f(c)=%8.5f\n',f(c))
fprintf('error=%8.5\n',e)
 Método de la bisección

f(x)= sin(2*x^2-6) =0 ; en (a1,b1)=(-2,2)

f(a)*f(b)>0

a+b
c= f(a)=2; f(b)=2
2
c −2+2 ; f(a)*f(b)>0 ,c1=a2
1= =0
2

c 0+2 ;f(a)*f(b)>0 , c2=a3

2= =1
2

c 1+2 ; f(a)*f(b)>0 .c3=a4

3= =1.5
2

……..

>> Bisección en matlab

Introduzca el valor de a: -2

Introduzca el valor de b: 2

Introduzca el número de iteraciones cont.: 6

Introduzca la función f(x)=sin(2*x^2-6)

k a b c f(c) error
1.0000 -2.0000 2.0000 0 0.2794 2.0000
2.0000 0 2.0000 1.0000 0.7568 1.0000
3.0000 1.0000 2.0000 1.5000 -0.9975 0.5000
4.0000 1.5000 2.0000 1.7500 0.1247 0.2500
5.0000 1.5000 1.7500 1.6250 -0.6584 0.1250
6.0000 1.5000 1.6250 1.5625 -0.8989 0.0625
Solucion: f(c)=-0.89887

c = 1.56250 Error=0.1

Método de bisección
x /2
f ( x )=cos ¿
)-1
a+b
c=
2
c −2+2 ; f(a)*f(b<0 fj(a)= -0.41615, f(b)=1 ,c1=b2
1= =0
2

c −2+0 ;f(a)*f(b)>0 , c2=b3

2= =−1
2

c 1+2 ; f(a)*f(b)>0
3= =1.5
2

>> Bisección en matlab

Introduzca el valor de a: -2

Introduzca el valor de b: 2

Introduzca el número de iteraciones cont:6

Introduzca la función f(x)=cos(x/2-1)

K a b c error
f(c)
1.000 - 2.000 0 0.2794 2.000
0 2.0000 0 0
2.000 0 2.000 1.000 0.7568 1.000
0 0 0 0
3.000 1.0000 2.000 1.500 -0.9975 0.500
0 0 0 0
4.000 1.5000 2.000 1.750 0.1247 0.250
0 0 0 0
5.000 1.5000 1.750 1.625 -0.6584 0.125
 0 0 0 0
6.000 1.5000 1.625 1.562 -0.8989 0.062
0 0 5 5

solución: f(c)=-0.20890

c = -1.56250 ; Error=>> 0.1

Newton Rapshon

function newton
f=input('ingrese la funcion:','s')%ingrese la funcion
f=inline(f);
fp=input('ingrese la derivada dela funcion:','s') % se ingresa ala
deribada dela funcion
fp=inline(fp);
xi=input('escriba el punto de inicio: ');%se escribe el punto de
inicio dela funcion
tol=input('ingrese el error maximo admisible: ')%se escribe la
tolerancia admitida por el programa
n=input('ingrese el numero maximo de iteraciones permitidas: ')% es el
numero de iteraciones permitidas
%se crean dos nuevas variables para guardar la funcion y la derivada
fx=f(xi);
fpx=fp(xi);
i=1;%sera el contador de iteraciones realizadas
error=tol+1;%sera el margen de error que utilizaremos
fprintf('\n n x f(x) fp(x) \n')%ess la
cabesera de la tabla
while i<=n && fx ~=0 && error>tol && fpx~=0 % esta sera la forma de
iterar los valore en la funcion
x=xi-(fx/fpx); %se calcula la susecion de aproximaiones
fx=f(x);% se asigna el valor dela suesion a fx
fpx=fp(x);%se le asiganel valor dela sucesionde x a fpx
fprintf('%1.0f %10.10f %10.10f %10.10f\n',i,x,f(x),fp(x))%se
imprime todas las iteraciones que se calcula en x
error=abs(x-xi);% determinamos el error absoluto estimado
xi=x;
 i=i+1; % segenera la proxima iteracion
end

if fx==0
fprint('\n La raizes: %1.10f\n\n',xi)
else if error<tol
fprintf('\n %1.10f es una aproximacion ala raiz con un error
maximo de %1.10f\n',xi,tol)
else if fpx==0
fprintf(' %1.10 es una posible raix mutiple.',xi)
else
fprintf('\n el mtodo fallo en %0.10 iteraciones\n\n',n)
end
end
end

Newton Rapshon

>> newton

f (0)
x2=x1-
f ´ (0)

ingrese la funcion:x.^3*sin(x)

f =x.^3*sin(x)

ingrese la derivada dela funcion:3*x.^2*sin(x)+x.^3*cos(x)

fp =3*x.^2*sin(x)+x.^3*cos(x)

escriba el punto de inicio: 2.5

ingrese el error maximo admisible: 0.001

tol =1.0000e-03

ingrese el numero maximo de iteraciones permitidas: 15

n = 15
9.4247779631 es una aproximacion ala raiz con un error maximo de 0.0010000000

Newton rapshon

>> newton

f (0)
x2=x1-
f ´ (0)

ingrese la funcion: sin(2*x.^2-6)

f =sin(2*x.^2- n x F(x) Fp(x) 6)

ingrese la derivada dela

1 9.7125027018 -259.9931024891 -958.8500005080
funcion:4*x

fp =4*x 2 9.4413517373 -13.9477491750 -845.9101369873

escriba el 3 9.4248632847 -0.0714325577 -837.2149420818 punto de

inicio: 2
4 9.4247779631 -0.0000019398 -837.1694716030
ingrese el error maximo
admisible: 0.001

tol = 1.0000e-03

ingrese el numero maximo de iteraciones permitidas: 15

n =15
 n x F(x) fp(x)
1 1.886337821 0.898587927 7.545351286
6 3 6
2 1.767246225 0.243835185 7.068984902
6 3 5
3 1.732752562 0.004862865 6.931010249
3 1 2
4 1.732050952 0.000001003 6.928203809
4 7

1.7320509524 es una aproximacion ala raiz con un error maximo de 0.0010000000

Falsa posición

Fx=input('ingrese la función :','s');%se ingresa la funcion

a=input('ingrese a : ');%se ingresa el intervalo inferior
c=input('ingrese c : ');% se ingresa el intervalo superior
e=input('ingrese el valor del error : '); % se ingresa la tolerancia

x=a;%se asigna la variable a->x

fa=eval(Fx);% se evalúa la función para ver si tiene algún error con
el intervalo inferior
x=c;% sea asigna la variable a->c
Fc=eval(Fx);
fprintf('\n%6s %7s %8s %10s %8s %8s %8s
\n','A','B','C','F(a)','F(b)','F(c)','|c-a|');%cabeseras de tabla
while abs(c-a)>e% se genera la iteración hasta el valor sea mayor a la
tolerancia
b=(c*Fa-a*Fc)/(Fa-Fc);%una función evaluada en el intervalo
superior con el inferior
x=b;
Fb=eval(Fx);%evalúa la función generada en b
fprintf('\n %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f
\n',a,b,c,Fa,Fb,Fc,abs(c-a));% se genreran la iteraciones
 if abs(Fc)<e
break;
if Fa*Fb<=0 %si la función de a por b es menor igual a 0 se
vuelve a evaluar
c=b;
Fc=Fb;
else %sino se genera la nueva iteración
a=b;
Fa=Fb;
end
end
end
fprintf('n El resultado será %.4f',b);
ezplot(Fx);%graficamos la función
grid on;

Falsa posición
Ingrese la funcion: sin(x.^2-4)

Ingrese a : -1

Ingrese c : 2

Ingrese el error : 0.001

A B C F(a) F(b) F(c) |c-a|

-1.0000 2.0000 2.0000 -0.1411 0.0000 0.0000 3.0000

El resultado sera 2.0000>>

>>falsaprop

Ingrese la función: sin(x)-csc(x+1)

Ingrese a: 0.5

Ingrese c: 0.7

Ingrese el error :0.00001

i A B C Fa Fb Fc (c-a)
1 0.5 0,700000 0,673667 -0,606404 0,091947 0,020923 0,020923
2 0.5 0,673667 0,667875 -0,606404 0,020923 0,004642 0,004642
3 0.5 0,667875 0,666600 -0,606404 0,004642 0,001024 0,001024
4 0.5 0,666600 0,666319 -0,606404 0,001024 0,000226 0,00022

, método dela secante codigo

clear
clc
f=input('Ingrese la funcion :','s');
x1=input('Ingrese el punto xi-1;');
x2=input('Ingrese el punto xi;');
err=input('Porcentaje de error;');
ezplot(f)
grid on
f=inline(f)
ea=100;
i=0;
fprintf('Iteracion Raiz\n')
while ea>err
xi=x2-(f(x2)*(x1-x2)/(f(x1)-f(x2)));
ea=abs(((x1-x2)*x1)/100);
fprintf('%f %8.5f\n',i,xi)
x1=x2;
x2=xi;
 i=i+1;
end
fprintf('\nRaiz de la funcion:%8.3f\nCalculada en
%4fIteraciones\n',xi,i)

Método dela secante

Ingrese la funcion :sin(x)-csc(x)+1

Ingrese el punto xi-1;0.5

Ingrese el punto xi;0.7

Porcentaje de error;0.00001

f = Inline function:

f(x) = sin(x)-csc(x)+1
 Iteración raíz
0.000000 0.67367
1.000000 0.66591
2.000000 0.66624
3.000000 0.66624

Raiz de la funcion: 0.666

Calculada en 4.000000Iteraciones

Método dela secante

Ingrese la funcion :ln(x.^2-1)

Ingrese el punto xi-1;-2

Ingrese el punto xi;2

Porcentaje de error;0.1

Error using inlineeval (line 14)

Error in inline expression ==> ln(x.^2-1)

Undefined function 'ln' for input arguments of type 'double'.

Error in inline/feval (line 33)

INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);

Error in ezploteval (line 51)

 z = feval(f,x(1));

Error in ezplot1ezplot1 (line 472)

[y, f, loopflag] = ezploteval(f, x);

Error in ezplot (line 144)

[hp, cax] = ezplot1(cax, f{1}, vars, labels, args{:});

Error in secante (line 7)

ezplot(f)

Método del a secante

Ingrese la funcion :x.^3-2*x.^2

Ingrese el punto xi-1;2

Ingrese el punto xi;4

Porcentaje de error;0.0001

f =Inline function:

f(x) = x.^3-2*x.^2

iteración raíz
0.00000
0 2.00000
1.00000
0 2.00000
2.00000 NaN
0

Raiz de la funcion: NaN

Calculada en 3.000000Iteraciones

Problemas libro de chapra métodos numericos

6.21 Use el programa que desarrolló en el problema 6.20 para resolver los problemas 6.12 ,
con una tolerancia de es = 0.001%.

Programacion;
clear , clc
cf=input('ingrese funcion a evaluar: ');
syms x
f=inline(cf);
derivada=diff(cf,x);
df=inline(derivada);
tol = input('ingrese tolerancia: ');
error = 50;
x=input('ingrese un valor inicial: ');
n =0;
disp(' n xi error')
while (error>tol)
fprintf('\t%i\t%3.5f\t%f\n', n, x, error);
n=n+1;
x=x-f(x)/df(x);
error= abs(f(x));
end

>> newton

ingrese la funcion:-x.^2+x+0.75

f =-x.^2+x+0.75

ingrese la derivada dela funcion:-2*x+1

fp =-2*x+1

escriba el punto de inicio: 1.2

ingrese el error maximo admisible: 0.001

tol = 1.0000e-03

ingrese el numero maximo de iteraciones permitidas: 15

n x f(x) fp(x)
1 1.5642857143 -0.1327040816 -2.1285714286
2 1.5019415149 -0.0038867992 -2.0038830297
3 1.5000018811 -0.0000037622 -2.0000037622
4 1.5000000000 -0.0000000000 -2.0000000000

1.5000000000 es una aproximacion ala raiz con un error maximo de 0.0010000000

6.22 El antiguo método de dividir y promediar, para obtener una apoximación de la raíz
cuadrada de cualquier número positivo, a, se formula del modo siguiente;

x=(( x−a)/ x )/2

Demuestre que éste es equivalente al algoritmo de Newton-Raphson.

x=√ a=; x=a2

f ( x )=x 2−a=0
f ' ( x )=2 x
Sustituyendo en l formula de newton
 2
x −a
)/2x
( x )=x−¿
Resolviendo termino tenemos:

x 2−a
x∗x−¿ ))/2
( x )=¿
x=(( x−a)/ x )/2
a
x 2−
x
x=
2

6.23 a) Aplique el método de Newton-Raphson a la función f(x) = tanh (x2 – 9) para evaluar su
raíz real conocida en x = 3. Use un valor inicial de x0 = 3.2 y haga un mínimo de cuatro
iteraciones. b) ¿Converge el método a su raíz real? Bosqueja la gráfica con los resultados para
cada iteración que obtenga

F ( X ) =tanh ( x 2−9 ) ; f ' ( x )=2 xsenh ( x 2−9 ) ;

x 0=3.2

F ( X ) =tanh ( 3.22−9 )=0.8454

 F ( X ) =2 xsenh ( x2 −9 ) =192.8278

0.8454
1=¿ 3.2− ( 192.8278 )=3.1956
x¿
f x 1=0.966

f ' x =−6.8808
1

0.0966
2=¿ 3.1956− ( −6.8808 )=3.1956
x¿

f x 2=0.9974

f ' x 2=39.5281

0.9974
3=¿3.3404− ( 39.5281 )=3.3151
x¿
Si converdge ala raíz real

6.24 El polinomio f(x) = 0.0074x4 – 0.284x3 + 3.355x2 – 12.183x + 5 tiene una raíz real entre 15
y 20. Aplique el método de Newton-Raphson a dicha función con valor inicial x0 = 16.15.
Explique sus resultados
 f ( x)=0.0074 x 4 – 0.284 x 3+ 3.355 x 2 – 12.183 x+5

f ' ( x )=0.0296 x 3−0.852 x2 +6.71 x−12.153

x 0=16.15

f ( x )=−9.5745

f ' ( x )=0.6463

x 2=16.15− ( −9.5745
0.6463 )
=30.9643

fx 2=1215.6767

f ' x 2=257.4707

x 3=30.9643− ( 1215.6767
257.4707 )
=2602426

fx 3=¿ 372.7825

fx 3=112.1034

x 4=26.2426− ( 372.7825
257.4707 )
=22.917

fx 4=110.7548

f ' x 4 =50.3903

110.7548
x 5=22.917−( ) =20.7119
50.3903

Va convergiendo hacia su raíz real

 X +1
6.25 Emplee el método de la secante con la función del círculo ¿ =16 , a fin de encontrar
¿
¿
una raíz real positiva. Haga que el valor inicial sea xi = 3 y xi–1 = 0.5. Aproxímese a la solución
del primer y cuarto cuadrantes. Cuando resuelva para f(x) en el cuarto cuadrante, asegúrese de
tomar el valor negativo de la raíz cuadrada. ¿Por qué diverge la solución?

X +1
¿
¿
¿

X +1
¿
¿
¿2
16−¿
F( X )=√ ¿
1RA ITERACION

I −1=¿ 0.5
X¿
I =¿
3
X¿

X I+1=1.6516

2DA ITERACION

I =¿ 1.6514
X¿
I −1=¿ 0.5
X¿
X I +1=4.142

LA SOLUCION DIVERGE
6.26 suponga el lector qe esta diseñando untanque
esférico fig.6.26 de almacenamiento de agus para un
poblado pequeño para unpais en desarroloo

3 R−h
V =π h2 ( )
3
Donde v= volumen h=profundidad del tanque en
el agua
R=radio del tanque
Si r=3m; ¿a que profundidaddebe llenare
eltanque de modo que contenga 30m3?
Aplicando newton rapshon

f (h)=(π /3) h3−V

X i +1=
πr x i2 − ( π3 ) x −V
i
3

2 πr x i−π x i2

X i +1=
π 10 xi2− ( π3 ) x −1000
i
3

2 πr 10i −π x i2

Ite xi x f (x) ¿ e ∨¿
r ¿
i
f¿
¿
0 10 1094.395 314.1593
1 6.516432 44.26917 276.0353 53.458%
2 6.356057 0.2858 272.4442 2.523%
3 6.355008 1.26E-05 272.4202 0.017%

Después de tres iteraciones, se determina que la raíz

es 6.355008 con un error relativo aproximado de
0.017%.