function [dy]=fty(t,y);
dy=2*t/y-y*t;