syrte

from numpy import pi, exp
from scipy.special import gamma, gammaincinv, gammainc

# Normalisation constant
def b(n):
    return gammaincinv(2*n, 0.5)

# Total luminosity of a 2D Sérsic profile
def sersic_total_lum(Ie, re, n):
    bn = b(n)

syrte

%%This is often best used in a latex file in two-column mode.  To do that, do:
%%
%%\usepackage{multicol}  % before \begin{document}
%% ... document ...
%%\begin{multicols}{2}
%%\bibliography{filename}{}
%%\bibliographystyle{apj_short_prop}
%%\end{multicols}

syrte

%pylab inline
import numpy as np
import pylab as plt
import GPy
import re

def get_equation(kern):
    match = re.search(r'(math::)(\r\n|\r|\n)*(?P<equation>.*)(\r\n|\r|\n)*', kern.__doc__)
    return '' if match is None else match.group('equation').strip()

syrte