PROGRAM dcomplex
c     TEST COMPLEX DOUUBLE PRECISION     
c     PROGRAM POKAZUJE, JAK UZYWAC ARYTMETYKI ZESPOLONEJ
c     ZAPIS LICZB ZESPOLONYCH: (x,y) OZNACZA x+iy
      implicit none
      integer pyt
      double precision ra,rb,rc,rd,ia,ib,ic,id,rx,y
      double complex a,b,c,d,x,z
   10 write(*,*)'PROGRAM DCOMPLEX'
      write(*,*)'PODAJ a,b,c - ZESPOLONE'
      read (*,*) a,b,c
      write(*,*)'complex => real (a=> re(a), im(a))'
      ra=real(a)
      ia=dimag(a)
      write(*,*)'GOTOWE'
      write(*,*)'re(a)=',ra,', im(a)=',ia
      write(*,*)'<ENTER>'
      read (*,*)
      write(*,*)'real=>complex'
      write(*,*)'PODAJ ra, ia -RZECZYWISTE'
      read (*,*)ra,ia
      a=dcmplx(ra,ia)
      write(*,*)'GOTOWE'
      write(*,*)'a=',a
      write(*,*)'<ENTER>'
      read (*,*)
      write(*,*)'DZIALANIA: a+b,a*b,a/b,sqrt(a),|a|'
      write(*,*)'GOTOWE'
      x=b+c
      write(*,*)'b+c=',x
      x=b*c
      write(*,*)'b*c=',x
      x=b/c
      write(*,*)'b/c=',x
      x=sqrt(b)
      write(*,*)'sqrt(b)=',x
      z=x*x
      write(*,*)'sqrt(b)^2=',z
      y=abs(b)
      write(*,*)'|b|=',y
      write(*,*)'<ENTER>'
      write(*,*)' 0 - STOP'
      write(*,*)' 1 - OD POCZATKU'
      read (*,*) pyt
      if(pyt.ne..0)then
      goto 10
      else
      continue
      endif
      stop 100
      end