function dfdx=q3_dfdx(x,p)
mu=p(1);

dfdx=mu-3*x.^2;
end