function solve_system

global a b f m
a = .9;
b = .9;
f = .5;
m = .3;
  
% options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-5]);
% [t,x] = ode45(@my_system,[0 100],[0.2 0.3], options);
[t,x] = ode45(@my_system,[0 100],[0.2 0.3]);
plot(t,x);
legend('u','v');
xlabel('time');
ylabel('population density');