function z=solfun(sol,t)
z=deval(sol,t);
z=z([1 3],:);