texto: metodos numericos para ecuaciones diferenciales ordinarias

xlabel(’tiempo (s)’);
ylabel(’velocidades(m/s)’);
grid;
for i=1:1:length(Y)
fprintf(’%i ’,i-1);fprintf(’%4.1f ’,t(i));
fprintf(’%8.4i ’,Y(1,i));fprintf(’%8.4i ’,Y(3,i));
fprintf(’%8.4i\n’,Y(5,i));
end
for i=1:1:length(Y)
fprintf(’%i ’,i-1);fprintf(’%4.1f ’,t(i));
fprintf(’%8.4i ’,Y(2,i));fprintf(’%8.4i ’,Y(4,i));
fprintf(’%8.4i\n’,Y(6,i));
end
Los desplazamientos de cada masa es
n ti x1 x2 x3
0 0.0 0000 0000 0000
1 0.1 4.0492e-003 1.8553e-004 1.5409e-007
2 0.2 1.2129e-002 1.7013e-003 9.8981e-005
3 0.3 1.9227e-002 5.4582e-003 7.9289e-004
4 0.4 2.3799e-002 1.0704e-002 2.8566e-003
5 0.5 2.7024e-002 1.5859e-002 6.7719e-003
6 0.6 3.0385e-002 2.0116e-002 1.2256e-002
7 0.7 3.4343e-002 2.3911e-002 1.8408e-002
8 0.8 3.8664e-002 2.8111e-002 2.4317e-002
9 0.9 4.3220e-002 3.3079e-002 2.9572e-002
10 1.0 4.8187e-002 3.8519e-002 3.4310e-002
11 1.1 5.3675e-002 4.3951e-002 3.8904e-002
12 1.2 5.9448e-002 4.9165e-002 4.3646e-002
13 1.3 6.5079e-002 5.4258e-002 4.8631e-002
14 1.4 7.0314e-002 5.9382e-002 5.3816e-002
15 1.5 7.5231e-002 6.4570e-002 5.9120e-002
16 1.6 8.0093e-002 6.9759e-002 6.4467e-002
17 1.7 8.5104e-002 7.4906e-002 6.9790e-002
18 1.8 9.0293e-002 8.0042e-002 7.5045e-002
40