# iserstep.m

Lie-group and Runge-Kutta-Munthe-Kaas solvers for the Lie equation. LG methods are implementations of Equation A.1, with minor typos fixed, from (http://dx.doi.org/10.1088/0305-4470/39/19/S07). The key difference from step.m function is that the Liouvillian can depend on the density matrix.

## Syntax

rho_b=iserstep(spin_system,L,rho_a,t,dt,method)

## Arguments

spin_system - Spinach data structure from create.m and basis.m constructors L - a handle to a function L(t,rho) that must take time and state vector, and return the evolution generator (in rad/s) of the Lie equation: d_rho/d_t = -i*L(t,rho)*rho rho_a - state vector at the start of the evolution period t - time at the start of the evolution, seconds dt - evolution time step, seconds method - 'PWCL', 'PWCM', 'RKMK4', or 'LG4', the latter one is recommended

## Outputs

rho_b - state vector at the end of the evolution time step

## See also

Time evolution functions, step.m, isergen.m, evolution.m, krylov.m

*Version 2.8, authors: Ilya Kuprov, Anupama Acharya, Andrew Graham*