function Yout = G_BD2(dt,U,Yn,Yp,Iinj)
 
Yout=-3/2*U+2*Yn-1/2*Yp+dt*F(U,Iinj);

  
