2楼: Originally posted by asaka at 2011-08-30 21:42:16:
每个苯环各取不共线的3个点，然后联立方程组求解。

      Program Dihedral_angel
      implicit double precision(a-h,o-z)
      dimension a(3,3,2),b(3,2)
      read(*,*) a
! calculate normal vectors
      do i=1,2
        do j=1,3
          a(j,2,i)=a(j,2,i)-a(j,1,i)
          a(j,3,i)=a(j,3,i)-a(j,1,i)
        end do
        b(1,i)=a(2,2,i)*a(3,3,i)-a(3,2,i)*a(2,3,i)
        b(2,i)=a(3,2,i)*a(1,3,i)-a(1,2,i)*a(3,3,i)
        b(3,i)=a(1,2,i)*a(2,3,i)-a(2,2,i)*a(1,3,i)
        tmp=sqrt(b(1,i)**2+b(2,i)**2+b(3,i)**2)
        do j=1,3
          b(j,i)=b(j,i)/tmp
        end do
      end do
! calculate dihedral angle
      c=b(1,1)*b(1,2)+b(2,1)*b(2,2)+b(3,1)*b(3,2)
      ang=acos(c)
      write(*,*) "Dihedral angle : ",ang,ang/atan(1.d0)*45.d0
      End

7楼: Originally posted by yqzhang6518 at 2011-09-04 18:47:35:
cartes 1.3.5