function z=fnv1p3b(t,y);
z=-(y.^3+y).*((3*(y.^2) +1).*t).^(-1);
