Testing properties of integrand

Might give insight into whats happening here. Seems like its mainly 1.0 - intgrand of int(c(r))

src/thermodynamics/thermo.rs

@@ -187,7 +187,7 @@ impl<'a> TDDriver<'a> {
                     &self.solutions,
                     &self.wv,
                     self.wu.as_ref().unwrap(),
-                    self.rism_kb_partial_molar_volume_density(),
+                    self.rism_kb_partial_molar_volume_density_2(),
                 );
                 let pressure = self.pressure();
                 let sfe = SFEs::new(&sfed, pressure, rism_kb_pmv, grid.dr);
@@ -373,6 +373,7 @@ impl<'a> TDDriver<'a> {
         let p_sum = self.total_density();
 
         compressibility * (1.0 - p_sum * integrand)
+        //1.0 - integrand
     }
 
     fn total_density(&self) -> f64 {