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