COMPLEX NUMBERS RCL CIRCUIT

c RLC under steady state  Schaum BEE probelm 3.7
      real L
      data v0,R, f , L, C/100.,100.,79.6,1.,5.E-6   /
      complex Zl , Zc ,ztotal, V , I, Vr, Vl, Vc
      pi=2.*asin(1.)
      radtodeg=180./pi
      omega=2.*pi*f
      v=cmplx(v0,0.)
      ZL=cmplx(0.,omega*L)
      Zc=cmplx(0.,-1./(omega*c))
      Ztotal=R+Zl+Zc
      i=V/ztotal
      vr=I*R
      vL=I*zl
      Vc=I*zc
      print*,'abs(i),angle(i) deg=',abs(i), radtodeg*atan(aimag(i)/
     $real(i))
      print*,'  '
      print*,'abs(vr),angle(vr) deg=',abs(vr), radtodeg*atan(aimag(vr)/
     $real(vr))
      print*,'  '
      print*,'abs(vl),angle(vl) deg=',abs(vL), radtodeg*atan(aimag(vL)/
     $real(vL))
      print*,'  '
      print*,'abs(vc),angle(vc) deg=',abs(vc), radtodeg*atan(aimag(vc)/
     $real(vc))
      print*,'  '
      stop
      end
      
       abs(i),angle(i) deg=  0.706206739 -45.0728798
 abs(vr),angle(vr) deg=  70.6206741 -45.0728798
 abs(vl),angle(vl) deg=  353.203339  44.9271164
 abs(vc),angle(vc) deg=  282.402771  44.9271126