% Written by Min-Hua Huang
% June 9, 2011
% This is a file to calculate the "h" function with arbitrary precision.
% input: h1(s,c,pre); s:the base parameter; c:the power parameter;
% pre:level of precision (number of singificant digits in matlab's vpa function).
function fv=h1(s,c,pre)
format long;  n=pre;  % n is the number of h's expansion terms. 
start1=vpa(1/(s+1),pre); sum2=vpa(start1,pre);
  for i=1:1:n
      start1=vpa(start1*(-c)/(s+1+i),pre);
      sum2=vpa(sum2+start1,pre);
  end
fv=vpa(sum2,pre); 
