! complx_func.f90  show complex functions

program complx_func
  implicit none
  complex :: z   = (2.0, 3.0);    ! complex constant
  complex :: z2  = cmplx(2.0,3.0) ! construct complex
  real    :: a

  print *, 'complx_func.f90 starting'

  print *, z, ' = z'   ! output 
  print *, z2, ' = z2'
  a = real(z)
  print *, 'real(z)=', a
  a = aimag(z)
  print *, 'aimag(z)=', a
  a = abs(z)
  print *, 'abs(z)=', a
  z2 = cmplx(4.0, 3.0)
  print *, 'cmplx(4.0,3.0)=', z2
  z2 = conjg(z)
  print *, 'conjg(z)=', z2
  z2 = sin(z)
  print *, 'sin(z)=', z2
  z2 = cos(z)
  print *, 'cos(z)=', z2
  z2 = log(z)
  print *, 'log(z)=', z2
  z2 = exp(z)
  print *, 'exp(z)=', z2
  z2 = sqrt(z)
  print *, 'sqrt(z)=', z2

  print *, 'complx_func.f90 finishing'
end program complx_func