#!/usr/bin/python# -*- кодировка: utf8 -*-импортировать matplotlib.pyplot как pltимпортировать matplotlib как mplимпортировать numpy как npиз импорта математики *code_website = 'http://commons.wikimedia.org/wiki/User:Geek3/mplwp'пытаться: импорт mplwpза исключением ImportError, эр: распечатать 'ImportError:', эр print 'Вам необходимо загрузить mplwp.py с', code_website выход(1)имя = 'mplwp_factorial_gamma_stirling.svg'фиг = mplwp.фиг_стандарт(mpl)xlim = -0,2, 3,2; рис.gca().set_xlim(xlim)илим = 0, 6,5; fig.gca().set_ylim(ylim)mplwp.mark_axeszero(рис.gca())из scipy.special импорт гаммаиз scipy.optimize импорт brentqf1 = лямбда x: гамма(1 + x)x1 = диапазон(4)y1 = [f1(xx) для xx в x1]plt.plot(x1, y1, 'o', label='n!', zorder=3)x2max = brentq(лямбда x: f1(x) - ylim[1], 0, xlim[1])x2 = np.linspace(xlim[0], x2max, 5000)y2 = [f1(xx) для xx в x2]plt.plot(x2, y2, label='$\Gamma(n+1)$', zorder=2, цвет=рис.gca().линии[-1].get_color())f3 = лямбда x: sqrt(2*pi*x)*(x/e)**xx3max = brentq(лямбда x: f3(x) - ylim[1], 0, xlim[1])x3 = np.linspace(0, x3max, 5000)y3 = [f3(xx) для xx в x3]plt.plot(x3, y3, label=ur'$\sqrt{2\pi n}(n/e)^n$', zorder=1)mpl.rc('легенда', borderaxespad=1.3)plt.legend(loc='верхний левый')plt.savefig(имя)mplwp.postprocess(имя)