function X = poliofun(b,tvec)

global initial_conditions
[tvec,Xsoln] = ode45(@polioode,tvec,initial_conditions,[],b);
X = Xsoln(:,3);