A simple numerical solution to the Townes profile in python using Eulers method.