### Ballistics

December 19, 2016 Digital computers were not the first electronic computers. There was a span of several decades in which analog computers were used to solve come common problems. It's fairly easy to implement the mathematical functions of integration and differentiation using operational amplifiers, the earliest of which were built using vacuum tubes. Notable electrical engineer, Robert Pease (1940-2011), who designed some early operational amplifiers while at George A. Philbrick Researches,[1] has written many articles about their history and design.[2-4]*An operational amplifier integrator.In operation, the capacitor is discharged to make the output voltage zero at time zero. Having the same value resistor in the non-inverting input eliminates drifting caused by input bias currents.(Created using Inkscape.)*

Since I'm a member of the baby boomer generation, my physics education included laboratory exercises in analog computation. The one exercise that I remember most is calculating the trajectory of a projectile, the results of the calculation being displayed on an x-y chart recorder. The interesting part of that exercise was seeing the difference in trajectory between a projectile launched on Jupiter as compared with Earth. The trajectory of a projectile, idealized without atmospheric drag, Coriolis force, etc., has been known since the time of Galileo, and it is derived by a simple addition of the motion arising from the projectile velocity and that of the fall of the projectile arising from gravity,

where (x,y) are the points that define the trajectory, θ is the launch angle from the horizontal, v

_{0}is the initial projectile velocity, v

_{h}is the horizontal component (x component) of the projectile velocity (v

_{h}= v

_{0}cos θ), and g is the gravitational acceleration, 9.80665 m/sec

^{2}. An example trajectory, as generated by a simple C language program (source code here), is shown below,

*An ideal trajectory without air drag and other confounding influences.Maximum range is achieved when the projectile is launched at 45 degrees.(Graphed using Gnumeric.)*

Isaac Newton used ballistics to elucidate the idea of an orbit of one body around another, such as the Moon's orbiting the Earth. In his 1728 book, "A Treatise of the System of the World," Newton imagined the thought experiment, illustrated below, in which you fired cannonballs at increasing velocities, finally reaching a point at which they fly into orbit around the Earth.[6] He then generalized this to the Moon's motion around the Earth.

*Newton's own illustration of his cannonball experiment. (From Isaac Newton, "A Treatise of the System of the World," F. Fayram, 1728.[6]*

The range

**R**of an ideal projectile is given as

**R = (v**, so you can hit a target at any range up to a maximum determined by the projectile's initial velocity

^{2}/g)sin(2θ)**v**by just setting the angle. It can be seen by inspection that the maximum range always occurs at an angle of 45 degrees, since the sine of 90° is one. That angle, however, is not the angle at which the length of the path of projectile travel is maximum. It's easy to calculate the path length of a projectile as a function of angle. You just use the trajectory example presented above and do a piece-wise integration of the path length. My source code for such a calculation can be found here,[6] and the results of this calculation are shown below with the maximum occurring at an angle near 1 radian (57.2958°).

*Trajectory path length as a function of launch angle.(Graphed using Gnumeric.)*

Joshua Cooper and Anton Swifton of the Department of Mathematics of the University of South Carolina have published a paper on arXiv in which they derive an equation to give the angle that produces the maximum projectile path length.[7] This angle, which is independent of the initial speed of the projectile and the value of the gravitational acceleration, is the solution of this interesting equation,

The easiest way to find the solution to an equation such as this is to use one of the many free and open source software computer algebra systems. Cooper and Swifton used SageMath to obtain an angle of about 0.9855 radians (56.47 degrees). I used my favorite method of refined increments (source code here) to obtain the value 0.985514738.

*Intersecting curves, csc(θ) and coth(csc(θ)).(Graphed using Gnumeric.)*

### References:

- GAP/R, George A. Philbrick Researches Archive.
- Application Brief R1, "Practical closed-loop stabilization of solid state operational amplifiers," Philbrick Archive, February 1, 1961.
- Bob Pease, "What's All This Transimpedance Amplifier Stuff, Anyhow? (Part 1)" Electronic Design, January 8, 2001.
- Bob Pease, "What's All This Julie Stuff, Anyhow?" Electronic Design, May 3, 1999.
- Isaac Newton, "A Treatise of the System of the World," F. Fayram, 1728, pp. 5-6.
- Thanks to Anton Swifton, who discovered an error in my original program.
- Joshua Cooper and Anton Swifton, "Throwing a Ball as Far as Possible, Revisited," arXiv, November 8, 2016.

*Permanent Link to this article*

Linked Keywords: Digital; computer; electronic; decade; analog computer; mathematical function; integral; integration; derivative; differentiation; operational amplifier; vacuum tube; electrical engineering; electrical engineer; Robert Pease (1940-2011); George A. Philbrick Researches; operational amplifier integrator; capacitor; discharge; voltage; resistor; non-inverting input; input bias current; Inkscape; baby boomer generation; physics; education; laboratory; analog computer; analog computation; calculation; calculating; trajectory; projectile; x-y chart recorder; Jupiter; Earth; atmospheric drag; Coriolis force; Galileo Galilei; velocity; gravitation; gravity; angle; horizontal plane; gravitational acceleration; C programming language; computer program; source code; trajectory.c; range; Gnumeric; Isaac Newton; ballistics; orbit; orbit of the Moon; Philosophiae Naturalis Principia Mathematica; A Treatise of the System of the World; thought experiment; Newton's cannonball; cannonball experiment; sine; length; arc length; iece-wise integration; radian; Joshua Cooper; Anton Swifton; Department of Mathematics; University of South Carolina; scientific literature; publish; arXiv; equation; independent variable; speed; free and open source software; computer algebra system; SageMath; solver.c; trigonometric function; csc; hyperbolic function; coth.

### Google Search

Free Downloads:

STEM artwork

for your holiday gifts

Latest Books by Dev Gualtieri

Thanks to Cory Doctorow of BoingBoing for his favorable review of Secret Codes!

Other Books

- High Energy Cosmic Rays - November 13, 2017

- Advanced Aluminum Alloys - November 6, 2017

- Joseph Polchinski - October 30, 2017

- Our Magnetic Universe - October 23, 2017

- Cavitation - October 16, 2017

- Pell Numbers - October 9, 2017

- Miniature Antennas - October 2, 2017

- Fizzy Graphene - September 25, 2017

- The First Angiosperm - September 18, 2017

- Noise Thermometry and the Boltzmann Constant - September 11, 2017

- Walking in the Rain - September 4, 2017

- Agitated Atoms - August 28, 2017

- Partial Solar Eclipse at New Jersey - August 24, 2017

- Magnetocapacitive Tunnel Junctions - August 21, 2017

- Tardigrades - August 14, 2017

- Roman Concrete - August 7, 2017

- Solar Spicules - July 31, 2017

- Schroeder Diffuser - July 24, 2017

- Rough Microparticles - July 17, 2017

- Robot Musicians - July 10, 2017

- Walter Noll (1925-2017) - July 6, 2017

- cosmogony - July 3, 2017

- Crystal Prototypes - June 29, 2017

- Voice Synthesis - June 26, 2017

- Refining Germanium - June 22, 2017

- Granular Capillarity - June 19, 2017

- Kirchhoff–Plateau Problem - June 15, 2017

- Self-Assembly - June 12, 2017

- Physics, Math, and Sociology - June 8, 2017

- Graphene from Ethylene - June 5, 2017

- Crystal Alignment Forces - June 1, 2017

- Martian Brickwork - May 29, 2017

- Carbon Nanotube Textile - May 25, 2017

- The Scent of Books - May 22, 2017

- Patterns from Randomness - May 18, 2017

- Terpene - May 15, 2017

- The Physics of Inequality - May 11, 2017

- Asteroid 2015 BZ509 - May 8, 2017

- Fuzzy Fibers - May 4, 2017

- The Sofa Problem - May 1, 2017

### Deep Archive

Deep Archive 2006-2008

**Blog Article Directory on a Single Page**