enddo
be=0.0D0
if (i.gt.2) phi(i+1)=beta(i-2,i-1,i,i+1)
+ if (i.gt.2) tauangle(3,i+1)=beta(i+nres-1,i-1,i,i+nres)
+ if (i.gt.2) tauangle(1,i+1)=beta(i-1+nres,i-1,i,i+1)
+ if (i.gt.2) tauangle(2,i+1)=beta(i-2,i-1,i,i+nres)
omeg(i)=beta(nres+i,i,maxres2,i+1)
theta(i+1)=alpha(i-1,i,i+1)
alph(i)=alpha(nres+i,i,maxres2)