function dwdt=vteamRRHS(t,w,p)

Rm=(p.Roff*(w-p.won)+p.Ron*(p.woff-w))/(p.woff-p.won);
vm=p.u(t)*Rm/(Rm+p.R); % memristor voltage
dwdt=0;
if vm>p.voff; dwdt=p.koff*(vm/p.voff-1)^p.aoff; end
if vm<p.von; dwdt=p.kon*(vm/p.von-1)^p.aon; end
if w>=p.woff || w<=p.won; dwdt = 0; end

if p.forceinside==1
  if w>p.woff; dwdt = -1e-6; end
  if w<p.won; dwdt = +1e-6; end
end
