before modifying space-time integration
error analysis function created
cleaning up
now crank-nicholson time integration