function myCalc
clear;
clc;
close all;
opt=odeset('RelTol',1E-5,'AbsTol',[1E-6 1E-6]);
[T Y]=ode45(@getdY,[0 100],[0.02 0.98],opt);
figure;
hold on;
grid on;
box on;
plot(T,dY(:,1),'r')
plot(T,dY(:,2),'b')
return
function dY=getdY(t,Y)
dY=zeros(2,1);
dY(1)=Y(1)*Y(2)-0.3*Y(1);
dY(2)=-Y(1)*Y(2);
return