[1] Yabushita K, Yamashita M and Tsuboi K 2007 [12] Bruyn J R and Walsh A M 2004 Can. J. Phys. J Phys A: Math and Theo 40 8403-8416 ISSN 82 439446 1751-8113 [13] Burden R L, Faires J D and Burden A M 2015

[2] Belgacem C H 2014 Eur J Physics 35 055025 Numerical Analysis (Brooks Cole)

[3] Turkyilmazoglu M 2016 Eur J Phys 37 ISSN [14] Hern´andez-Saldan˜a H Energy losses in projec0143-0807 tile motion under constant, linear and quadratic

[4] Morales D A 2016 Acta Mech 227 1593-1607 drag work in progress ISSN 0001-5970 [15] Ferna´ndez-Chapou J L, Salas-Brito A L and

[5] Stewart S M 2012 Eur J Phys 33 149-166 ISSN Vargas C A 2004 Am J Phys 72 1109 ISSN 0002- 0143-0807 9505

[6] Dittrich T, Ha¨nggi P, Ingold G L, Kramer B, [16] Thomas G, Weir M, Hass J and Giordano F 2004 Scho¨n G and Zwerger W 1998 Quantum Trans- Calculus (Addisson-Wesley) port and Dissipation (Wiley-VCH) ISBN 978- [17] Carrillo-Bernal M A, Mancera-Pin˜a P E, 3527292615 Cerecedo-Nun˜ez H H, Padilla-Sosa P, Nun˜ez

[7] Razavy M 2006 Classical and Quantum Dissi- Y´epez H N and Salas-Brito A L 2014 Am. J. pative Systems (Imperial College Press) ISBN Phys. 82 707 ISSN 0002-9505 9781783260393 [18] Stewart S 2006 Int J Math Educ Sci Technol 37

[8] Jones S E, Caipen T L and Butson G J 1991 411-431 International Journal of Applied Engineeering [19] Stewart S M 2011 Eur J Phys 32 L7-L10 ISSN Education 7 321-327 0143-0807

[9] Jones S E 1991 The projectile motion problem [20] Hern´andez-Saldan˜a H 2010 Eur J Physics 31 including the effects of constant thrust and drag 1319-1329 Challenges of a changing world vol 1 and 2 ed [21] In one dimension this term only changes the sign Educ A S E The organization (Univ. of New Or- of the force in order to keep it in opposition to leans, New Orleans, LA: Amer. Soc. Engineering the movement. Education) pp 1839-1841 proceeding paper [22] Since the launching angle remain in the first

[10] MacMillan W 1927 Theoretical Mechanics: quadrant Static and the Dynamics of a Particle (McGraw- [23] Calculations for values of b0 larger than 0.25 reHill) quire of a better selection of the initial condition

[11] Barnes H A 1999 J Non-Newt Fluid Mech 81 as we can be seen for b0 = 0.75 in figure 4(a) 133-178