Integration implemented and working.

Results are however wrong, as most of my plots are mostly linear, except
for the pressure. Do not know what is wrong yet.