function f=forwardDiff(fcn,x,h)
% compute the forward difference of 
% function fcn at x
f=(feval(fcn,x+h)-feval(fcn,x))/h;
