из импорта numpy * из импорта matplotlib.pyplot * из импорта mpl_toolkits.axes_grid.axislines SubplotZero pdf = лямбда c , mu , x : sqrt ( c / ( 2. * pi )) * exp ( -0.5 * c / ( x - mu ) ) / ( x - mu ) ** 1.5 fig = figure ( figsize = ( 6 , 5 )) ax = SubplotZero ( fig , 111 ) fig . add_subplot ( ax ) ax . grid ( True ) ax . minorticks_on ( ) x = logspace ( - 1,6 , 1,6 ) ax . loglog ( x , pdf ( .5 , 0. , x ), c = 'r' , lw = 1,5 , метка = r '$c=1/2$' ) #'$c={}^1\!/_2$') ax . loglog ( x , pdf ( 1. , 0. , x ), c = 'k' , lw = 1,5 , метка = r '$c=1$' ) ax . loglog ( x , pdf ( 2. , 0. , x ), c = 'c' , lw = 1,5 , метка = r '$c=2$' ) ax . loglog ( x , pdf ( 4. , 0. , x ), c = 'b' , lw = 1.5 , метка = r '$c=4$' ) ax . loglog ( x , pdf ( 8. , 0. , x ), c = 'g' , lw = 1.5 , метка = r '$c=8$' )ax.set_ylim ( 1e-4 , 1 ) ax.set_xlim ( 0.05 , 15 ) ax.legend ( loc = 4 , frameon = False ) ax.set_xticks ( ( 0.1 , 1.0 , 10. )) ax.set_xticklabels ( ( r ' 0.1 ' , r ' 1.0 ' , r '10 .0 ' )) ax.set_xlabel ( r ' x $-\!\rightarrow $ ' )рис . сохранить рис ( "Levy0_LdistributionPDF.svg" , bbox_inches = "tight" , pad_inches = .15 )