Crank-Nicolson simulation with starting point (0.5, 0.5, 0.5, 0.5)