function y = f321(x) % % y = f321(x) % % a simple oscillating function a = 3; b = 2; c = 1; y = (x-a).*(cos(b*x) + sin(c*x));