Enter below a function using normal fortran notation (e.g., .5*g*x^2, note: "^" = "**" is OK; the mass of 1 amu is 1.66054E-27 kg = 1.66054 × 10-27 kg; trig functions like sin and asin = sin-1 = arcsin are designed for radian measure; sqrt(1/(mu0*epsilon0)))