function [t, y1, y2] = ode45_2(strfunc, tspan, y0)
[t,y] = ode45(@odefunc2, tspan, y0, [], strfunc) ;
y1 = y(:,1);
y2 = y(:,2);
