echo off;
odefun=@frosc;
tspan=0:0.01:10;
x0=[1 1];
options=[];
p=[0 1];
hold on
for k=1:2
    [t x]=ode45(odefun,tspan,x0,options,p(k));
    subplot(1,2,k);
    plot(x(:,1),x(:,2));axis equal;grid on
end
subplot(1,2,1);
legend('R=0');
subplot(1,2,2);
legend('R=1');
hold off