Fixed wrong assumption with pressure.

I though pressure from rhs was only gas pressure. This is wrong, its the
total pressure. Updated software.

Also updated to the best initial condition in integrateAndPlot.