Python code

importnumpyasnp
importmatplotlib.pyplotasplt
fromscipyimport stats
frommatplotlib.tickerimport MultipleLocator
plt.rcParams.update({
 "text.usetex": True,
 "font.family": "serif",
})
x = np.linspace(-4, 4, 500)
beta_list = [0.0, 0.25, 0.50, 0.75, 1.00]
y = [stats.levy_stable.cdf(x, 0.5, beta) for beta in beta_list]
fig, ax = plt.subplots()
c = ['black', 'blue', 'lime', 'red', 'magenta']
for i, data in enumerate(y):
 ax.plot(x, data, color=c[i], label=r'$\beta=' + str(beta_list[i]) + r'$')
 
ax.xaxis.set_minor_locator(MultipleLocator(.25))
ax.xaxis.set_major_locator(MultipleLocator(2))
ax.yaxis.set_minor_locator(MultipleLocator(.025))
ax.tick_params(axis='both', which='both', direction='in', top=True, right=True)
ax.set(ylim=(0, 1), xlim=(-4, 4))
plt.plot([], [], ' ', label=r'$\alpha = 0.5$')
plt.plot([], [], ' ', label=r'$c=1$')
plt.plot([], [], ' ', label=r'$\mu=0$')
plt.legend(frameon=False)
fig.savefig('Levyskew_distributionCDF.svg', bbox_inches='tight')