Non capisco dove sta il problema, perché non ti limiti semplicemente a sommare -bv/m all'accelerazione al momento di integrare con Eulero? La v già ce la dovresti avere, se, come dovresti fare, stai riducendo il sistema di equazioni differenziali del second'ordine ad uno vettoriale al prim'ordine come spiegato qui...