In the same way that resistors are inserted in a simulation scheme by stamps in order to establish MNA matrix, it is necessary to insert memristors and memristive systems in MNA or Tableau matrix.
at an instant, a memristor can be thought as a resistor or a nonlinear resistor. If there is a only a charge dependency, R=M(q). Solve the node voltage matrix by taking the memristor as a resistor, find either charge or flux of the memristor, then using either euler or a runge-kutta method, calculate the charge q(t+dt)=q(t)+i*dt (for a memristive system calculate x(t+dt)=x(t)+f(x,v,t)*dt) then update M(q) as M(q(t+dt) for all memristors and solve the MNA again for the time step. If the memristor model is more complex, a nonlinear solution procedure is needed.