  We present a novel solution method for It\^o stochastic differential
equations (SDEs). We subdivide the time interval into sub-intervals, then we
use the quadratic polynomials for the approximation between two successive
intervals. The main properties of the stochastic numerical methods, e.g.
convergence, consistency, and stability are analyzed. We test the proposed
method in SDE problem, demonstrating promising results.
