anchor.diff:          sparserow(i+1) = counter + 1
        end if
    enddo ! do i
  -
  +call enter_exit(sub_name,2)
    
   end subroutine Mat_to_CSR
   
  
positive.code: !           1- As of now it will print out the constants from Fung's paper (Pulmonary microvascular impedance-1972) 
  
  
  use arrays, only:dp
  
  ! Parameters:
  
negative.code: integer :: i
  integer :: factors(8)
  integer :: NonZeros
  character(len=60) :: sub_name
  
      
  sub_name = 'calc_cap_imped'
  call enter_exit(sub_name,1)
  
  
  rep = ha**4 - hv**4
  dx = 1.0/(N_nodes - 1.0)
  h(1) = ha**3
  h(2) = (ha**4 - 2*rep*dx)**(0.75)
  h(3) = (ha**4 - (N_nodes-1)*rep*dx)**(0.75)
  h(4) = (ha**4 - (N_nodes-2)*rep*dx)**(0.75)
  h(5) = (ha**4 - (N_nodes-3)*rep*dx)**(0.75)
  
  n = N_Nodes*2 + 2
  call Matrix(N_nodes, ha, hv, omega, stiff1, stiff2, RHS1, RHS2)
  call Mat_to_CC(stiff1,N_nodes,sparsecol,sparserow,sparseval,NonZeros)
   nrhs = 1
    ldb = n
  