made metabolic heat generation input volume specific, and also created a script to compare the odesolver set up with the analytical setup