% External Function function y=g(x) y = 3.*exp(x-pi/2.)+sin(x);