3.233235838169363e-1,6.791528011378658e-1,8.753479805169189e-1,
1.695040627613259e-1,5.420702855281473e-1,7.985688950544638e-1};
double[] xx=buildVector(x),yy=buildVector(y),zz= new double[N],uu= new double[N],vv= new double[N],c=buildVector(N,1.);
for (i=0;i<N;i++) {
zz[i]=gam.gammp(a[i],x[i]); // Test gammp
uu[i]=gam.gammq(a[i],x[i]); // Test gammq
vv[i]=gam.invgammp(zz[i],a[i]); // Test invgammp
}
System.out.printf("Gamma.gammp: Maximum discrepancy = %f\n", maxel(vecsub(zz,yy)));