In order to avoid the intensive computation, we would like to push the limits higher to z value more than 0.5 so that we can accommodate the conductive materials with high k values. Implicit and Crank Nicholson methods are used to stabilize that.