% Written by Min-Hua Huang
% June 9, 2011
% This is a file to calculate the forward difference function.
% input: difference(s,c,order); s is the base parameter; c is the power 
% parameter; order is the order of the finite difference.
function fv=difference(s,c,order)
sum=0;  
    for i=0:1:order
        temp=nchoosek(order,i)*(-1)^(order-i)*h(s+i,c,0);
        sum=sum+temp;
    end
fv=sum;

