function f=dy0(t,x)
f=x^2-t^2;