function dy = fun(t,y)
dy = -2*t*y^2;
