texto: metodos numericos para ecuaciones diferenciales ordinarias

end
T=[t];
y=[w];
w(j)=c(j)
for i=1:N,
end
for j=1:m,
end
k1(j)=h*fs(t,w(1),w(2),j);
for j=1:m,
end
k2(j)=h*fs(t+h/2,w(1)+(1/2)*k1(1),w(2)+(1/2)*k1(2),j);
for j=1:m,
end
k3(j)=h*fs(t+h/2,w(1)+(1/2)*k2(1),w(2)+(1/2)*k2(2),j);
for j=1:m,
end
k4(j)=h*fs(t+h,w(1)+k3(1),w(2)+k3(2),j);
for j=1:m,
end
w(j)=w(j)+(k1(j)+2*k2(j)+2*k3(j)+k4(j))/6
t=a+i*h;
T=[T;t];
y=[y;w]
input(’pulse’)
fprintf(’%4.4f ’,T(i+1));
fprintf(’%12.6i ’,w(1));fprintf(’%12.6i\n’,w(2));
plot(T,y(:,1))
xlabel(’tiempo (s)’);
ylabel(’posicion (m)’);
grid
Programa fs.m
function y1=fs(t,x1,x2,j)
31