& dcosomicron(j,1,2,i)-fac0*(dc_norm(j,i-3)-scalp*
& dc_norm(j,i-1+nres))/vbld(i-1+nres)
dtauangle(j,2,3,i)=-1/sing*dcostau(j,2,3,i)
& dcosomicron(j,1,2,i)-fac0*(dc_norm(j,i-3)-scalp*
& dc_norm(j,i-1+nres))/vbld(i-1+nres)
dtauangle(j,2,3,i)=-1/sing*dcostau(j,2,3,i)