/* Incomplete Gamma Function Gam_x(a) (from 0 to x) */
new; cls;
print Gx(2,3);

library pgraph;
graphset;
x=seqa(0.05,0.05,99)';
a=seqa(0.05,0.05,99);
xlabel("x");
ylabel("a");
surface(x,a,Gx(a,x));



proc Gx(a,x);
   retp(cdfgam(a,x).*gamma(a));
endp;