function dydx=myproblem(y,x)
dydx=x*y.^2+y;