void zgam(double[],double[],double*,long);