      program plotarch
      implicit doubleprecision (a-h,o-z)
      pi=4*datan(1.0d0)
    
      open(2,file="cc")
      read(2,*)
      read(2,*)
      read(2,*)r
      close(2)
    
      open(2,file="cylinder.pov")
      
      dz=1.0
      nz=11
      dfi=dz/r
      nfi=int(2*pi/dfi)+1
      dfi=2*pi/dble(nfi)
      z1=-5.0
      z2=5.0
     
      do i=1,nfi
       fi=(i-1)*dfi
       x=r*dcos(fi)
       y=r*dsin(fi)
       write(2,'("cylinder{<",f8.4,",",f8.4,",",f8.4,">,<",
     :                             f8.4,",",f8.4,",",f8.4,",> 0.02 ",
     :                             "pigment {Black}}")')
     :    x,y,z1,x,y,z2
      enddo
     
      do i=1,nz 
        z=-5+(i-1)*dz
        nfi=100
        dfi=2*pi/100.0
        do j=1,nfi
          fi1=(j-1)*dfi
          fi2=j*dfi
          x1=r*dcos(fi1)
          y1=r*dsin(fi1)
          x2=r*dcos(fi2)
          y2=r*dsin(fi2)
          write(2,'("cylinder{<",f8.4,",",f8.4,",",f8.4,">,<",
     :                             f8.4,",",f8.4,",",f8.4,",> 0.02 ",
     :                             "pigment {Black}}")')
     :    x1,y1,z,x2,y2,z
        enddo
       enddo

       close(2)
       stop
       end
