- dcostau(j,1,1,i)=fac1*dcosomicron(j,2,1,i-1)+fac3*
- & dcosomicron(j,2,1,i-1)-fac0*(dc_norm(j,i-1)-scalp*
- & (-dc_norm(j,i-2+nres)))/vbld(i-2+nres)
- dtauangle(j,1,1,i)=-1/sing*dcostau(j,1,1,i)
- dcostau(j,1,2,i)=fac1*dcosomicron(j,2,2,i-1)+fac2*
- & dcostheta(j,1,i)+fac3*dcosomicron(j,2,2,i-1)+fac4*
+ dcostau(j,1,1,i)=fac1*dcosomicron(j,2,2,i-1)+fac3*
+ & dcosomicron(j,2,2,i-1)-fac0*(dc_norm(j,i-1)-scalp*
+ & (dc_norm2(j,i-2+nres)))/vbld(i-2+nres)
+ dtauangle(j,1,1,i)=-1.0d0/sing*dcostau(j,1,1,i)
+ dcostau(j,1,2,i)=fac1*dcosomicron(j,2,1,i-1)+fac2*
+ & dcostheta(j,1,i)+fac3*dcosomicron(j,2,1,i-1)+fac4*