Metodos Numericos
Metodos Numericos
p 3
p* 1.732
p*
, donde
P 3 1.732050808
P* 1.732
t
ValorVerdadero ValorAproximado
ValorVerdadero
1.732050808 1.732
100% 0.002933%
1.732050808
p*
2. Determina el mayor intervalo en que debe estar
con un error relativo de
103
P 3 1.732050808
t
ValorVerdadero ValorAproximado
ValorVerdadero
10 3
3 P*
3
p* 3 3 103 1.730318757
p* 3 3 103 1.733782858
Intervalo:
P* 1.730318757,1.733782858
para aproximar
p 3
4e 12
3 2
26.56081949
16
3 2
16
P* 26.56
Et ValorVerdadero ValorAproximado
Et 26.560819 26.56 8.195 104
ValorVerdadero ValorAproximado
ValorVerdadero
26.560819 26.56
100% 0.003085%
26.560819
f x 25 x3 6 x 2 7 x 88
x 1.75
4. Evala el polinomio
en
. Utiliza aritmtica de
cuatro dgitos por redondeo y truncamiento. Evala el error relativo
porcentual verdadero en cada caso. Repite el clculo anterior, pero expresa
el polinomio en su forma anidada.
Valor
x3
25x3
x2
Valor exacto
1.7
5
5.3593
75
133.98437
5
3.062
5
Valor truncamiento
4 cifras
1.7
5
5.359
133.9
3.062
6x 2
18.37
5
-18.37
7x
f(x)
12.2
5
39.85937
5
12.2
5
39.78
1.7
5
5.359
134.0
3.063
-18.38
12.2
5
39.87
39.859375 39.78
100% 0.199137593%
39.859375
t
Truncamiento:
39.859375 39.87
100% 0.026656213%
39.859375
Redondeo:
f(x) ((25 x 6) x 7) x 88
Forma Anidada
39.859375
39.85
39.86
VALOR EXACTO
Truncamiento
Redondeo
Truncamiento:
39.859375 39.85
100% 0.023520188%
39.859375
Redondeo:
39.859375 39.86
100% 0.001568012544%
39.859375
x 0.19
y 0.08
Calculo del error para x:
t
0.1945908156 0.19
100% 2.36%
0.1945908156
0.08305599852 0.08
100% 3.68%
0.08305599852
1
1 x
x
1/(1+x)
es
0.75
0.57142
857
0.00000
05
Iter
0
1
2
3
4
1/(1+x)
1
0.25
0.8125
0.39062
5
0.70703
125
ea
300
69.2307
692
108
44.7513
812
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
0.46972
656
0.64770
508
0.51422
119
0.61433
411
0.53924
942
0.59556
293
0.55332
78
0.58500
415
0.56124
689
0.57906
483
0.56570
137
0.57572
397
0.56820
702
0.57384
473
0.56961
645
0.57278
766
0.57040
925
0.57219
306
0.57085
52
0.57185
86
0.57110
605
0.57167
046
0.57124
715
0.57156
50.5197
505
27.4783
264
25.9584
57
16.2961
675
13.9239
253
9.45551
031
7.63293
225
5.41472
261
4.23294
357
3.07702
125
2.36228
187
1.74086
824
1.32292
395
0.98244
521
0.74230
345
0.55364
529
0.41696
535
0.31174
91
0.23435
979
0.17546
143
0.13176
948
0.09872
954
0.07410
202
0.05554
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
463
0.57132
652
0.57150
511
0.57137
117
0.57147
162
0.57139
628
0.57145
279
0.57141
041
0.57144
219
0.57141
836
0.57143
623
0.57142
282
0.57143
288
0.57142
534
0.57143
1
0.57142
675
0.57142
994
0.57142
755
0.57142
934
0.57142
8
0.57142
9
0.57142
825
0.57142
881
0.57142
839
565
0.04167
66
0.03124
768
0.02344
125
0.01757
785
0.01318
513
0.00988
787
0.00741
645
0.00556
203
0.00417
17
0.00312
867
0.00234
656
0.00175
989
0.00131
993
0.00098
994
0.00074
246
0.00055
684
0.00041
763
0.00031
322
0.00023
492
0.00017
619
0.00013
214
9.9106E
-05
7.433E05
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
0.57142
871
0.57142
847
0.57142
865
0.57142
851
0.57142
861
0.57142
854
0.57142
86
0.57142
855
0.57142
859
0.57142
856
0.57142
858
0.57142
857
0.57142
858
0.57142
857
0.57142
857
0.57142
857
0.57142
857
0.57142
857
5.5747E
-05
4.181E05
3.1358E
-05
2.3518E
-05
1.7639E
-05
1.3229E
-05
9.9218E
-06
7.4414E
-06
5.581E06
4.1858E
-06
3.1393E
-06
2.3545E
-06
1.7659E
-06
1.3244E
-06
9.933E07
7.4498E
-07
5.5873E
-07
4.1905E
-07
f x e x
en
xi 1 1
para
en cada caso.
xi 0.25
. Emplea
xi 1 1
h xi 1 xi 1 0.25 0.75
f ( xi 1 ) f (1) e 1 0.367879441
Valor Verdadero
Para n=0
f ( xi 1 ) e .0.25 0.778800783
0.367879441 0.778800783
100% 111.7%
0.367879441
Para n=1
f ( xi 1 ) f ( xi ) f (') ( xi )h
f ( xi 1 ) e .0.25 (e0.25 )(0.75) 0.1947001958
t
0.367879441 0.1947001958
100% 47.07499955%
0.367879441
Para n=2
f ( xi 1 ) f ( xi ) f (') ( xi ) h
f ('') ( xi )h 2
2!
(e 0.25 )(0.75) 2
0.413737916
2!
0.367879441 0.413737916
100% 12.46562593%
0.367879441
Para n=3
f ( xi 1 ) f ( xi ) f (') ( xi ) h
f ('') ( xi )h 2 f (''') ( xi ) h3
2!
3!
0.3589784856
2!
3!
0.367879441 0.3589784856
100% 2.419530452%
0.367879441
Para n=4
f ( xi 1 ) f ( xi ) f ' ( xi ) h
2!
3!
4!
0.369245879
2!
3!
4!
0.367879441 0.369245879
100% 0.371436359%
0.367879441
Para n=5
f ( xi 1 ) f ( xi ) f ' ( xi ) h
2!
3!
4!
5!
0.36770577
2!
3!
4!
5!
0.367879441 0.36770577
100% 0.047208672%
0.367879441
Para n=6
f ( xi 1 ) f ( xi ) f ' ( xi )h
2!
3!
4!
5!
6!
2!
3!
4!
5!
(e 0.25 )(0.75) 6
0.3678982836
6!
t
0.367879441 0.3678982836
100% 0.005121950%
0.367879441
g x sen x
g x sen x
xi 1
xi
g( xi 1 ) g ( ) sen( ) 0.866025403
3
3
Valor verdadero
h xi 1 xi
3 4 12
Para n=0
g( xi 1 ) g ( xi )
g( xi 1 ) sen( ) 0.707106781
4
0.866025403 0.707106781
100% 18.35034185%
0.866025403
Para n=1
g( xi 1 ) g ( xi ) g (') ( xi )h
xi 1
, en
xi
para
en cada caso.
g( xi 1 ) sen( ) cos( )( ) 0.892226903
4
4 12
g( xi 1 ) g ( ) sen( ) 0.866025403
3
3
Valor verdadero
0.866025403 0.892226903
100% 3.025488618%
0.866025403
Para n=2
g( xi 1 ) g ( xi ) g (1) ( xi )h
g (2) ( xi ) h 2
2!
sen( ) ( ) 2
4
12 0.867994736
g( xi 1 ) sen( ) cos( )
4
4 12
2!
g( xi 1 ) g ( ) sen( ) 0.866025403
3
3
Valor verdadero
0.866025403 0.867994736
100% 0.227398988%
0.866025403
Para n=3
g( xi 1 ) g ( xi ) g (1) ( xi )h
g (2) ( xi ) h 2 g (3) ( xi ) h3
2!
3!
g( xi 1 ) sen( ) cos( )
4
4 12
sen( ) ( ) 2 cos( ) ( )3
4
12
4
12 0.865880079
2!
3!
g( xi 1 ) g ( ) sen( ) 0.866025403
3
3
Valor verdadero
0.866025403 0.865880079
100% 0.016780570%
0.866025403
Para n=4
g( xi 1 ) g ( xi ) g (1) ( xi )h
2!
3!
4!
g( xi 1 ) sen( ) cos( )
4
4 12
g( xi 1 ) g ( ) sen( ) 0.866025403
3
3
Valor verdadero
0.866025403 0.866018482
100% 0.00079916%
0.866025403
Para n=5
g( xi 1 ) g ( xi ) g (1) ( xi )h
2!
3!
4!
5!
sen( ) ( ) 2 cos( ) ( ) 3
4
12
4
12
g( xi 1 ) sen( ) cos( )
4
4 12
2!
3!
sen( ) ( ) 4 cos( ) ( )5
4
12
4
12 0.866025728
4!
5!
g( xi 1 ) g ( ) sen( ) 0.866025403
3
3
Valor verdadero
0.866025403 0.866025728
100% 0.00003752%
0.866025728
Para n=6
g( xi 1 ) g ( xi ) g (1) ( xi ) h
2!
3!
4!
5!
6!
sen( ) ( ) 2 cos( ) ( )3
4
12
4
12
g( xi 1 ) sen( ) cos( )
4
4 12
2!
3!
g( xi 1 ) g ( ) sen( ) 0.866025403
3
3
Valor verdadero
0.866025403 0.866025411
100% 9.237604315 10 7 %
0.866025728
f x ln 1 x
, para valores de
x 2 x3 x 4
...
2
3
4
Solucin:
%........................................................................
%Programa para calcular la expansion por serie de maclaurin para la
funcion
%y=Ln(1+x).
%Elaborado por:Cynthia Susana Fierro Rendon
%MetodosNumericos
%Ingenieria en Electronica y Telecomunicaciones
%Plataforma: MATLAB
%........................................................................
clc
clearall
clear()
disp('
disp('
disp('
f x 1 x
1 x 1 kx
k
k k 1 x 2
2!
k k 1 k 2 x 3
3!
, para valores de
...
Solucin:
%........................................................................
..
%Programa para calcular la expansin por seriedemaclaurin para la funcin
%y=Sqrt(1+x).
%Elaborado por:Cynthia Susana Fierro Rendn
%MetodosNumericos
%Ingenieria en Electronica y Telecomunicaciones
%Plataforma: MATLAB
%........................................................................
..
clc
clearall
clear()
disp('
Universidad Politcnica del Estado de Morelos');
disp('
MetodosNumericos ');
disp('
Calculo de la funcion y=Sqrt(1+x) por la serie de
Maclaurin');
%Se introducen datos por el usuario
x=input('Ingrese el valor de x que sea -1 <x < 1: ');
es=input('Ingrese el valor de estimado (es):
');
%Calcula la serie de Maclaurin
%k=n
Iter=0;
suma0=0;
k=0;
yv=sqrt(1+x); % valor verdadero de la funcion
ea=100; %error aproximado
fprintf('iter
y
ea\n');
%Realiza suma y calculo de terminos de la serie
while (ea> es);
suma=suma0+(((-1)^(k))*((x ^(k))*factorial(2*k)))/((1(2*k))*((factorial(k)^2)*(4^k)));
ea=abs((suma-suma0)/suma)*100;
fprintf('%3d \t %11.7f \t
%11.7f \n', Iter, suma, ea);
if(ea>es);
suma0=suma;
k=k+1;
Iter=Iter+1;
else
fprintf('\n\n El valor aproximado (ea) de la funcion es: %11.7f',suma);
fprintf('\n\n El valor verdadero o real es:
%11.7f',yv);
end;
end;
disp(' ');
disp('