anchor.diff:  !                - SparseRow ====> An array formed with values to represent row storage. Includes the index numbers of start of each row (Size of array = MatrixSize + 1)
   
   
  -
  +use diagnostics, only: enter_exit
   use arrays, only:dp
   
   !Paramters:
  
positive.code: integer, intent(out) :: NonZeros
  !Local Variables:
  integer :: i, j, counter
  
  
  allocate(sparsecol(2*nn+3)) !allocation of sparsecol
  NonZeros = 0 !Initialisation for NonZeros
  
negative.code: !           1- As of now it will print out the constants from Fung's paper (Pulmonary microvascular impedance-1972) 
  
  
  use diagnostics, only: enter_exit
  use arrays, only:dp
  
  ! Parameters:
  