diff options
Diffstat (limited to 'Rinteg.f')
-rwxr-xr-x | Rinteg.f | 4 |
1 files changed, 2 insertions, 2 deletions
@@ -35,7 +35,7 @@ c****************************************************************************** c(n)=0. b(n)=(f(n)-f(n1))/(x(n)-x(n1)) a(n)=f(n)-x(n)*b(n) - if(n.eq.2)return + if(n==2)return do 1 j=2,n1 j1=j-1 d=(f(j)-f(j1))/(x(j)-x(j1)) @@ -49,7 +49,7 @@ c****************************************************************************** c(3)=0. b(3)=(f(4)-f(3))/(x(4)-x(3)) a(3)=f(3)-x(3)*b(3) - if(c(j).eq.0.)go to 2 + if(c(j)==0.)go to 2 j1=j+1 wt=abs(c(j1))/(abs(c(j1))+abs(c(j))) a(j)=a(j1)+wt*(a(j)-a(j1)) |