ode solver

% Solve a system of differential equations using ODE45
%minster
% Initial conditions: y1(0) = 2; y2(0) = 3
% Solve on t between 0 and 20.
clear all
close all
[T Y] = ode45(@fty, [0, 30], [4,0,0,0.2]);

plot(T,Y(:,1),’-‘)
hold on
plot(T,Y(:,2),’o’)
xlabel(‘t’); ylabel(‘y)’); title(‘y1 — and y2 ooo’);