A Newton–Cotes method for quantum stochastic differential equations

We present an implicit Simpson’s rule for the numerical computation of quantum stochastic differential equations. We derive the method, present convergence results and a numerical example. Our method avoid the high cost of high-order methods as well the cost of obtaining starting values for multi-step methods.