应该没有吧。我曾经根据一本参考书的算法，编了一个求INCOMPLETE GAMMA 的小程序。因为和网络计算器给出的结果不一样（只是部分相同），也就没敢放到blog里去。
京剧

研究了一下发现：虽然SAS中没有直接的不完全Gamma函数，
但可以变通地得到不完全gamma函数．介绍二种方法．
一是借助cdf函数，再有就是用积分函数．

proc iml;
  t=t({1 5 10 20});
  a=3;  lambda=4;
  yt=cdf("gamma", t, a, lambda);
  Igamma1=lambda**a * gamma(a) * yt;
  
start fun(x) global(a, lambda);
  v=x**(a-1) * exp(-x/lambda);
  return(v);
finish fun;
  call quad(Igamma2, "fun", {0}//t);
  Igamma2=cusum(Igamma2);
  print Igamma1 Igamma2;
  quit;