Check out the basic functionality with the following commands:
import flm
# Initialize fractional levy motion with trajectories of length 1000, hurst parameter = 0.4, alpha = 1.75
fractional_levy_motion = flm.FLM(1000, 0.4, 1.75)
# generate 10 independent trajectories
fractional_levy_motion.generate_realizations(10)
# plot the marginal distribution of hops
fractional_levy_motion.plot_marginal(show=False, bounds=(-5, 5))
# calculate and plot the autocorrelation function of the increments
fractional_levy_motion.plot_autocorrelation()
# plot the first 2 realizations
fractional_levy_motion.plot_trajectory([0, 1])
# calculate and plot the mean squared displacement of the trajectories
fractional_levy_motion.plot_msd(show=True)
Copyright (c) 2020, Benjamin Coscia
Project based on the Computational Molecular Science Python Cookiecutter version 1.1.