| 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 |
|
|