With f(v)=bv, we can use an arithmetic mean for bvn a la Crank-Nicolson schemes. un=un−1+Δtvn−12,vn+12=(1+b2mΔt)−1(vn−12+Δtm−1(Fn−12f(vn−12)−s(un))).