function f=forwardDiffF(fcn,x,h)
% compute the forward difference of 
% function fcn at x
% k. Ming Leung, October 20,2008.
% New version: x can be a matrix, or h a matrix, but not both!

f=(feval(fcn,x+h)-feval(fcn,x))./h;
