anchor.diff:    write ( *, '(a)' ) 'D_SAMPLE:'
     write ( *, '(a)' ) '  Normal end of execution.'
     write ( *, '(a)' ) ''
  -!  do i = 1, n
  -!    write ( *, '(2x,g14.6,2x,g14.6)' ) RHS1(i)
  -!  end do
   
   call Mat_to_CC(stiff2,N_nodes,sparsecol,sparserow,sparseval,NonZeros)
    nrhs = 1
  
positive.code:     integer :: MatrixSize,NonZeros,submatrixsize,ngen,i
      real(dp) :: area,Q01_mthrees,sheet_number
      character(len=60) :: sub_name
      real(dp) :: ha,hv,omega
      
      sub_name = 'cap_flow_ladder'
      call enter_exit(sub_name,1)
  
      ha = 4.0
      hv = 1.0
      omega = 10.0
      call calc_cap_imped(ha,hv,omega)
  !     Number of non-zero entries in solution matrix. 
        NonZeros=3
        do i=2,cap_param%num_symm_gen
  
negative.code:   write ( *, '(a)' ) '  Computed solution:'
    write ( *, '(a)' ) ' '
  
  write (*,*) 'C_1 = ',RHS1(2),'+ i(',RHS1(N_Nodes+2),')'
  deriv1(1) = (h(2)*RHS1(3) - h(1)*RHS1(1))/(2*dx)
  !c3 = (h(3) * H1R(3) - h(1) * H1R(1))/(2*dx);
  !c4 = (h(3) * H1I(3) - h(1) * H1I(1))/(2*dx);
  deriv1(2) = (h(2)*RHS1(N_nodes+4) - h(1)*RHS1(N_Nodes+2))/(2*dx)
  write (*,*) 'C_2 = ',deriv1(1),'+ i(',deriv1(2),')'
  
  !
  !  Free memory.
  