For an intense field model, the time-dependent Schrødinger equation with initial and boundary conditions can be discretized into the inhomogeneous linear canonical equation by substituting the symmetric difference quotient for the partial derivative. As the general solution of its homogeneous equation and the particular solution of the inhomogeneous equation can be generalized by the symplectic transformation, it is a reasonable numerical method to use the symplectic scheme. To prove its utility, a simple example is described using the symplectic scheme and RK method, and compared with the exact solution. The results show that the solution using the symplectic scheme can preserve the intrinsic properties of the equations after a long evolution, but RK method cannot.