function Yout = G_BE(dt,U,Yn,Iinj)
 
Yout=-U+Yn+dt*F(U,Iinj);
  
end

