00001 subroutine xdgamma (x, result) 00002 external dgamma 00003 double precision x, result, dgamma 00004 result = dgamma (x) 00005 return 00006 end