24d23
< #define EPSERR   3.E-3
67,79c66
< 	static double lagrangepoly(const int i, const double z, const double* zpts, const int np)
<     {
< 		double h=1.0;
< 		for(int k=0;k<i; ++k)
< 		{
< 			h = h * (z - zpts[k])/(zpts[i]-zpts[k]);
< 		}
< 		for(int k=i+1;k<np;++k)
< 		{
< 			h = h * (z - zpts[k])/(zpts[i]-zpts[k]);
< 		}
<     	return h;
<     }
---
> 
1374c1361
<         if (fabs(dz) < EPSERR) return lagrangepoly(i, z, zgj, np);
---
>         if (fabs(dz) < EPS) return 1.0;
1450c1437
<         if (fabs(dz) < EPSERR) return lagrangepoly(i, z, zgrj, np);
---
>         if (fabs(dz) < EPS) return 1.0;
1534c1521
<         if (fabs(dz) < EPSERR) return lagrangepoly(i, z, zgrj, np);
---
>         if (fabs(dz) < EPS) return 1.0;
1618c1605
<         if (fabs(dz) < EPSERR) return lagrangepoly(i, z, zglj, np);
---
>         if (fabs(dz) < EPS) return 1.0;