program t dimension f(5),s(5),c(5) pi=4.0*atan(1.0) do iu=1,5 if(iu.eq.1)then f(1)=0.0 *pi/180.0 f(2)=25.15 *pi/180.0 f(3)=-40.9 *pi/180.0 f(4)=40.9 *pi/180.0 f(5)=-25.14*pi/180.0 endif if(iu.eq.2)then f(5)=0.0 *pi/180.0 f(1)=25.15 *pi/180.0 f(2)=-40.9 *pi/180.0 f(3)=40.9 *pi/180.0 f(4)=-25.14*pi/180.0 endif if(iu.eq.3)then f(4)=0.0 *pi/180.0 f(5)=25.15 *pi/180.0 f(1)=-40.9 *pi/180.0 f(2)=40.9 *pi/180.0 f(3)=-25.14*pi/180.0 endif if(iu.eq.4)then f(3)=0.0 *pi/180.0 f(4)=25.15 *pi/180.0 f(5)=-40.9 *pi/180.0 f(1)=40.9 *pi/180.0 f(2)=-25.14*pi/180.0 endif if(iu.eq.5)then f(2)=0.0 *pi/180.0 f(3)=25.15 *pi/180.0 f(4)=-40.9 *pi/180.0 f(5)=40.9 *pi/180.0 f(1)=-25.14*pi/180.0 endif d=pi*2.0/5.0 do 2 i=1,5 c(i)=cos(real(i-1)*d) 2 s(i)=sin(real(i-1)*d) s1=0.0 s2=0.0 do 1 i=1,5 do 1 j=1,5 s1=s1+f(i)*(c(i)*s(j)*c(j)-s(i)*c(j)*c(j)) 1 s2=s2+f(i)*(s(i)*s(j)*c(j)-c(i)*s(j)*s(j)) p=atan(-s1/s2) if(p.lt.0)p=p+pi s3=0.0 s4=0.0 do 3 i=1,5 if(iu.eq.1)then s3=s3+f(i)*cos(P+real(i-1)*d) s4=s4+ (cos(p+real(i-1)*d))**2 endif if(iu.eq.2)then s3=s3+f(i)*cos(P+real(i)*d) s4=s4+ (cos(p+real(i)*d))**2 endif if(iu.eq.3)then s3=s3+f(i)*cos(P+real(i+1)*d) s4=s4+ (cos(p+real(i+1)*d))**2 endif if(iu.eq.4)then s3=s3+f(i)*cos(P+real(i+2)*d) s4=s4+ (cos(p+real(i+2)*d))**2 endif if(iu.eq.5)then s3=s3+f(i)*cos(P+real(i+3)*d) s4=s4+ (cos(p+real(i+3)*d))**2 endif 3 continue fm=s3/s4 write(6,*)'P = ',p *180.0/pi write(6,*)'fm = ',fm *180.0/pi write(6,*) enddo end