function [t, y] = ode45_1(strfunc, tspan, y0)

[t, y] = ode45(@odefunc, tspan, y0, [], strfunc) ;

