function f = ford(t, y)
f = [y(2) y(3) sin(t*y(1))+y(3)*sqrt(1+sqrt(y(2)^2))]';