00001 subroutine xgamma (x, result) 00002 external gamma 00003 real x, result, gamma 00004 result = gamma (x) 00005 return 00006 end