PYTHON: We want to find the position, as a function of time, of a damped harmonic oscillator.

The equation of motion is: m d²x/dt² = -kx b dx/dt, x(0) = 0.5, dx/dt (0) = 0

Take m = 0.25 kg, k = 100 N/m, and b = 0.1 N.s/m.

Solve x(t) for t in the interval [0, 10T], where T = 2/, and ² = k/m.

Divide the interval into N = 10⁴ intervals: t_i = i*h, I = 0, 1, 2, , N, where h = 10T/N.