import matplotlib.pyplot as plt
import numpy as np
from leniax.growth_functions import poly_quad4
x = np.linspace(0, 1, 1000)
m = 0.3
s = 0.05
params = np.array([m, s])
y = poly_quad4(params, x)
plt.plot(x, y)
plt.axhline(y=0, color='grey', linestyle='dotted')
plt.title(r'Poly_quad4: m=%.2f, s=%.2f'%(m, s))
plt.show()