function y = fp321(x) % % y = fp321(x) % % The derivative of f123(x) a = 3; b = 2; c = 1; y = (cos(b*x) + sin(c*x)) + (x-a).*(-b*sin(b*x) + c*cos(c*x));