subroutine prop_fullturn(ie) implicit real*8 (a-h,o-z) include 'tempstor.inc' c tempbeg is the full turn at the start. propogate to ie. tempie= c temp_invXtempbegXtemp dimension S(4,4),temp_inv(4,4) data S/0.0,-1.0,2*0.0,1.0,6*0.0,-1.0,2*0.0,1.0,0.0/ c temp is symplectic so temp_inv=-SXtempXS do i=1,4 do j=1,4 temp_inv(i,j)=0. do k=1,4 do l=1,4 temp_inv(i,j)=temp_inv(i,j)-S(i,k)*temp(l,k)*S(l,j) end do end do end do end do c do i=1,4 do j=1,4 tempie(i,j)=0. do k=1,4 do l=1,4 tempie(i,j)=tempie(i,j)+temp(i,k)*tempbeg(k,l)*temp_inv(l,j) end do end do end do end do return end