arccos(y(t))'=(a+b* y(t))^2