- dtauangle(j,3,1,i)=-1/sing*dcostau(j,3,1,i)
- dcostau(j,3,2,i)=fac1*dcosomicron(j,2,2,i-1)+fac2*
- & dcosomicron(j,1,1,i)+fac3*dcosomicron(j,2,2,i-1)+fac4*
+ dtauangle(j,3,1,i)=-1.0d0/sing*dcostau(j,3,1,i)
+ dcostau(j,3,2,i)=fac1*dcosomicron(j,2,1,i-1)+fac2*
+ & dcosomicron(j,1,1,i)+fac3*dcosomicron(j,2,1,i-1)+fac4*