Add the output of delta_T + change the output of the volume error + chose delta_T=delta_x we can halve delta t untill delta t < dx^2