program iraich5
      dimension f(5),s(5),c(5),g(5)

      pi=4.0*atan(1.0)
      d=pi*4.0/5.0
      do 2 i=1,5
      c(i)=cos(real(i-1)*d)
2     s(i)=sin(real(i-1)*d)
      
88    Write(6,6005)
6005  format(' Input f1 f2 f3 f4 f5 in degs: ',$)
      read(5,*) f(1),f(2),f(3),f(4),f(5)

      do 4 i=1,5
4     f(i)=f(i)*pi/180.0

      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=atan2(s1,-s2)
      if(p.lt.0)p=p+2.0*pi

      s3=0.0
      s4=0.0
      do 3 i=1,5
      s3=s3+f(i)* cos(P+real(i-1)*d)
3     s4=s4+     (cos(p+real(i-1)*d))**2

      fm=s3/s4

      o=0.0
      write(6,*)' old  new'
      do 6 i=1,5
      g(i)=fm*cos(P+real(i-1)*d)
      write(6,6003)f(i)*180.0/pi,g(i)*180.0/pi
6003  format(2f10.2)
6     o=o+(g(i)-f(i))**2
      write(6,*)

      write(6,6001)fm *180.0/pi,p*180.0/pi,o*180.0**2/pi**2
6001  format(' fm = ',f10.3,' P = ',f10.3,' S = ',f15.2)
      write(6,*)
      
      goto 88
      end