!****************************************************
program Electric_Fields
!****************************************************
 implicit none
 integer,parameter :: P=20     !max number of charges
 real,dimension(P) :: X,Y,Q
 integer           :: N
!-------------  SET CHARGE DISTRIBUTION ----
 N    =  2
 X(1) =  1.0
 Y(1) =  0.0
 Q(1) =  1.0
 X(2) = -1.0
 Y(2) =  0.0
 Q(2) = -1.0
!-------------  DRAWING LINES  -------------
 call eline(0.0, 0.5,X,Y,Q,N)
 call eline(0.0, 1.0,X,Y,Q,N)
 call eline(0.0, 1.5,X,Y,Q,N)
 call eline(0.0, 2.0,X,Y,Q,N)
 call eline(0.0,-0.5,X,Y,Q,N)
 call eline(0.0,-1.0,X,Y,Q,N)
 call eline(0.0,-1.5,X,Y,Q,N)
 call eline(0.0,-2.0,X,Y,Q,N)

end program Electric_Fields
