function Ejercicio3
[t1, yy1] = ode45(@ej3, [0 40], [1 0], [], 5);
plot(t1, yy1(:,1))
hold on
[t2, yy2] = ode45(@ej3, [0 40], [1 0], [], 40);
plot(t2, yy2(:,1))
[t3, yy3] = ode45(@ej3, [0 40], [1 0], [], 200);
plot(t3, yy3(:,1))
legend('c=5','c=40','c=200')
end

function dy=ej3(t,y,c)
dy=[y(2);(-(c*y(2)+20*y(1)))/20];
end

%%Otra forma
% function Ejercicio3
% [t1, yy1] = ode45(@ej3_1, [0 40], [1 0]);
% plot(t1, yy1(:,1))
% hold on
% [t2, yy2] = ode45(@ej3_2, [0 40], [1 0]);
% plot(t2, yy2(:,1))
% [t3, yy3] = ode45(@ej3_3, [0 40], [1 0]);
% plot(t3, yy3(:,1))
% end
%
% function dy=ej3_1(t,y)
% dy=[y(2);(-(5*y(2)+20*y(1)))/20];
% end
%
% function dy=ej3_2(t,y)
% dy=[y(2);(-(40*y(2)+20*y(1)))/20];
% end
%
% function dy=ej3_3(t,y)
% dy=[y(2);(-(200*y(2)+20*y(1)))/20];
% end

Published with MATLAB® R2018b