function y=dfunc(x)
y=3*x.^2+8*x;