Skip to content

Instantly share code, notes, and snippets.

@richpsharp
Created April 12, 2021 20:53
Show Gist options
  • Save richpsharp/fe1bb9030168a925aed2b8b57a4b5975 to your computer and use it in GitHub Desktop.
Save richpsharp/fe1bb9030168a925aed2b8b57a4b5975 to your computer and use it in GitHub Desktop.
is beta/m broken?
import numpy
import matplotlib.pyplot
slope_in_radians_list = numpy.linspace(0.089, 0.785, 100)
beta_list = []
m_list = []
for slope_in_radians in slope_in_radians_list:
beta = ((numpy.sin(slope_in_radians) / 0.0896) /
(3 * numpy.sin(slope_in_radians)**0.8 + 0.56))
m_list.append(beta/(1+beta))
beta_list.append(beta)
matplotlib.pyplot.plot(slope_in_radians_list, beta_list)
matplotlib.pyplot.plot(slope_in_radians_list, m_list)
matplotlib.pyplot.show()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment