51 |
integer ifail !=-1 if singular matrix error, =0 if not singular |
integer ifail !=-1 if singular matrix error, =0 if not singular |
52 |
integer jfail !=0 if determinant can be evaluated, =-1 if determinat is probably too small, =+1 if too large |
integer jfail !=0 if determinant can be evaluated, =-1 if determinat is probably too small, =+1 if too large |
53 |
|
|
54 |
parameter (big=1.e4) !just a number greater than C(npoints,3) |
parameter (ibig=1.e4) !just a number greater than C(npoints,3) |
55 |
double precision xxc(big),zzc(big),rrr(big) !centres and radii to be averaged |
double precision xxc(ibig),zzc(ibig),rrr(ibig) !centres and radii to be averaged |
56 |
|
|
57 |
double precision tmp1,tmp2,tmp(npoints) !temp variables |
double precision tmp1,tmp2,tmp(npoints) !temp variables |
58 |
|
|