[Vout]=V_backward(vm,dt,stimVe,StimIin,lamda_dt_matrix,Gamma_dt,Omega_dt)

function [Vout]=V_backward(vm,dt,stimVe,StimIin,lamda_dt_matrix,Gamma_dt,Omega_dt)
%zeng.dummyvm(:,i),dt,stimVe,StimIin,zeng.setup.lamda_dt_matrix,Gamma_dt,Omega_dt

Vout=sparse(dt*(lamda_dt_matrix+diag(Gamma_dt))+diag(ones(size(Omega_dt,1),1)))\(vm+dt*(Omega_dt-lamda_dt_matrix*stimVe)+StimIin);

%matrix_real=sparse(lamda_matrix+ diag(1+Gamma));
%Vout=matrix_real\(vm-lamda_matrix*ve+Omega+Istim);