Tarea Ultima Calculo To

UNIVERSIDAD TECNICA DE
MANABI NOMBRE:
TAREA
PREGUNTA 1:
Dada la tabla de altura contra tiempo de un cuerpo
 Calcule el valor aproximado de la velocidad y la aceleración para cada valor de t
 Realice un modelo de regresión para v – t
t (s) y (m) v (m/s) a (m/s^2)

0 500 -4,9 y (m)
600
1 495,1 -14,7 -9,8
2 480,4 -24,5 -9,8 500
f(x) = − 4.9 x² + 3.25188439971586E-14 x + 500
3 455,9 -34,3 -9,8 400 R² = 1
4 421,6 -44,1 -9,8
300
5 377,5 -53,9 -9,8
6 323,6 -63,7 -9,8 200
7 259,9 -73,5 -9,8
100
8 186,4 -83,3 -9,8
9 103,1 -93,1 -9,8 0
0 2 4 6 8 10 12
10 10 1
v (m/s)
20
f(x) = − 21.5016471690217 ln(x) − 15.9330250158396
R² = 0.263804398618296
0
0f(x) = 0.0132817
2 x⁶ − 0.350653
4 x⁵ + 3.474071
6 x⁴ − 15.95671
8 x³ + 33.74278
10x² − 37.4320112x
R² = 0.997289867439918
-20
-40
-60
-80
-100
PREGUNTA 2
1
Resolver la integral ∫ ecos (x) dx empleando los métodos

0
Trapecio con paso h =0.02 2,38171110
x y Area 0,52 2 0,047395
0 2,718281828 0,054360 2,35775213
0,02 2,717738245 0,054338 0,54 8 0,046910
0,04 2,716108362 0,054295 2,33323358
0,06 2,713394787 0,054230 0,56 7 0,046414
0,08 2,709601853 0,054143 2,30818765
0,58 1 0,045908
0,1 2,704735607 0,054035
2,28264672
0,12 2,698803795 0,053906
0,6 7 0,045393
0,14 2,691815832 0,053756
2,25664332
0,16 2,683782777 0,053585 0,62 9 0,044869
0,18 2,674717301 0,053394 2,23021001
0,2 2,664633642 0,053182 0,64 4 0,044336
0,22 2,653547567 0,052950 2,20337930
0,24 2,641476325 0,052699 0,66 9 0,043796
0,26 2,628438591 0,052429 0,68 2,17618364 0,043248
0,28 2,614454413 0,052140 2,14865526
0,3 2,599545154 0,051833 0,7 3 0,042695
0,32 2,583733428 0,051508 2,12082619
0,34 2,567043033 0,051165 0,72 9 0,042136
0,36 2,549498884 0,050806 2,09272816
0,38 2,531126943 0,050431 0,74 9 0,041571
0,4 2,511954145 0,050040 2,06439253
0,76 4 0,041002
0,42 2,492008322 0,049633
2,03585023
0,44 2,471318128 0,049212
0,78 6 0,040430
0,46 2,44991296 0,048777
2,00713174
0,48 2,427822883 0,048329 0,8 AREA TOTAL:3 2,34152669
0,039854
0,5 2,405078545 0,047868 1,97826699
Ir
0,82 7 0,039276
2,34157
1,94928536
0,84 5 0,038695
1,92021559
0,86 4 0,038113
1,89108576
0,88 7 0,037530
1,86192326
0,9 7 0,036947
1,83275473
Error 0,00184943 0,92 9 0,036364
1,80360605
0,94 6 0,035781
Simpson 1/3 y 3/8 dos veces. Obtenga el 1,77450229
error relativo cometido en las tres 0,96 4 0,035200
aproximaciones 1,74546770
x f(x) Simpson 1/3 una vez 0,98 7 0,034620
0 2,71828183 2,3425203 1 1,7165257
VALOR REAL IR: 2,34157
0,5 2,40507854
Error 0,040583
1 1,71652570
VALOR REAL IR: 2,34157
Error 0,002795
x f(x) Simpson 1/3 dos veces

0 2,71828183 2,341635
0,25 2,63507705 Simpson 3/8 una
x0,5 2,40507854
f(x) vez Ir 2,34157
0,75 0 2,07858811
2,71828183 2,341997 Er 0,018228
10,333333 1,71652570
2,57270261
0,666667 2,19435304
1,000000 1,7165257
Ir 2,34157
Er 0,000359
x f(x) Simpson 3/8 dos veces

0,00 2,71828183 2,341578
0,11 2,70157118
0,22 2,65225477
0,33 2,57270289
0,44 2,46662263
0,56 2,33872955
0,67 2,19435395
PREGUNTA 3:
0,78 2,03903197
Sea el problema del valor inicial
{
0,89 1,87812838
dy
1,00 1,71652714 = y (x−5)
dx
y ( 1 )=0.5
y (2 )=?
Aplique los métodos de Euler, Runge Kutta orden 2 Y Runge kutta orden 4. Tome el paso h =
0.01
Compare los resultados con la solución analítica y explique qué método es más efectivo.
Métodos de Euler
>> syms x
>> x=1; y=0.5;
>> for x=0:0.01:1.99
f=y*(x-5)
y=y+0.01*f
disp([x+0.01,y])
end
f=
-4.3547e-03
y=
1.4032e-03
2.0000 0.0002
%El resultado por el método de Euler es y(2)= 0.014032
>> y=dsolve('Dy= y*(x-5)','y(1)=0.5','x')
y=
(exp((x*(x - 10))/2)*exp(9/2))/2
>> y=vpa(subs(y,2),5)
y=
0.015099
%El valor real es y(2)=0.015099
>> Er=abs((0.015099-0.014032)/0.015099)*100
Er =
7.0667
%Error es 7.0667
Runge Kutta orden 2

>> syms x
>> x=1; y=0.5;
>> for x=0:0.01:1.99
k0=y*(x-5); k1=(y+0.01*k0)*((x+0.01)-5);
y=y+0.01/2*(k0+k1);
disp([x+0.01,y])
end
0.0100 0.4756
2.0000 0.0002
%El resultado de RK2 es y(2)= 0.0002
Runge kutta orden 4
>> syms x
>> x=1; y=0.5;
>> for x=0:0.01:1.99
k1=y*(x-5);k2=(y+0.01/2*k1)*((x+0.01/2)-5);k3=(y+0.01/2*k2)*((x+0.01/2)-
5);k4=(y+0.01/2*k3)*((x+0.01)-5);
y=y+0.01/6*(k1+2*k2+2*k3+k4);
disp([x+0.02,y])
end
0.0200 0.4755
2.0000 0.0002
%El resultado de RK4 es y(2)= 0.0002
>> y=vpa(subs(dsolve('Dy=y*(x-5)','y(1)=0.5','x'),2),5)
y=
0.015099
%El valor real es y(2)= 0.015099
Valor Real 0.01509869

Error Euler 6.60632634

Error R-K2 0.0749235

Error R-K4 0.01509869
El método más efectivo es el de Runge Kutta de orden 4, ya que, el resultado

comparado con el de Matlab fue lo mismo.
PREGUNTA 4:
Resuelva la ecuación diferencial dada, empleando los métodos de sustitución y R K 2.
2
d y dy −x dy
2
−4 x − y=5 e con y ( 0 )=1; ( 0 )=−1 Halle y(1) tome h=0.01
dx dx dx
y’’=4x*y’+y+5 e−x
P2’=4t*P2+P1+5exp(-x)
Y’’=P1’’=P2’’
Y(0)=P1(1)
Y’(0)=P2(-1)
Y(1)=?
F de p1=p2
>> syms x
>> p1=1; p2=-1; h=0.01;
>> for x=0:h:1-h
k01=p2; k02=(4*(x)*(p2))+p1+5*exp(-x);
k11=(p2+h*k02); k12=(4*(x+h)*(p2+h*k02))+(p1+h*k01)+5*exp(-x+h);
p1=p1+h/2*(k01+k11);p2=p2+h/2*(k02+k12);
disp([x+h,p1])
end
0.0100 0.9903
1.0000 4.1667

Tarea Ultima Calculo To

Cargado por

Copyright:

Formatos disponibles

Tarea Ultima Calculo To

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tarea Ultima Calculo To

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD TECNICA DE

t (s) y (m) v (m/s) a (m/s^2)

Resolver la integral ∫ ecos (x) dx empleando los métodos

x f(x) Simpson 1/3 dos veces

x f(x) Simpson 3/8 dos veces

>> y=dsolve('Dy= y*(x-5)','y(1)=0.5','x')

Runge Kutta orden 2

Valor Real 0.01509869

Valor Real 0.01509869

Valor Real 0.01509869

El método más efectivo es el de Runge Kutta de orden 4, ya que, el resultado

También podría gustarte