diff options
author | Joseph Hunkeler <jhunkeler@gmail.com> | 2021-08-04 15:18:08 -0400 |
---|---|---|
committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2021-08-04 15:18:08 -0400 |
commit | 70fd14e98bc476a84aef352eaf063cfcba44ab83 (patch) | |
tree | 2737f5c1c21c76824f49c6a1b9c7061952b646f5 /Rinteg.f | |
parent | c396a4bc7a3282de9d2191de0e1076fec9df0783 (diff) | |
download | moog-70fd14e98bc476a84aef352eaf063cfcba44ab83.tar.gz |
Use modern equality operators
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)) |