%NQUAD\_RAM
%
function [n,IR]=nquad_Ram(t,n0,dn,nofdig)   
format long
n=n0-dn; 
err=1;
while err>.5*10^(-nofdig)
  n=n+dn;
  ab=r_jacobi01(n);
  xw=gauss(n,ab); x=xw(:,1); w=xw(:,2);
  s=0;
  for k=1:n
    s=s+w(k)*gam_upinc(x(k),t)/gamma(x(k));
  end
  if n==n0
    IR0=s*exp(t);
  else
    IR=s*exp(t);
    err=abs((IR-IR0)/IR);
    IR0=IR;
  end
end
