function yprime = firstode(t,y);
yprime = t*y;
