def Pn(x):
Pn = {}<EOL>Pn['<STR_LIT:0>'] = P0(x)<EOL>Pn['<STR_LIT:1>'] = P1(x)<EOL>Pn['<STR_LIT:2>'] = P2(x)<EOL>Pn['<STR_LIT:3>'] = P3(x)<EOL>Pn['<STR_LIT:4>'] = P4(x)<EOL>Pn['<STR_LIT:5>'] = P5(x)<EOL>Pn['<STR_LIT>'] = P6(x)<EOL>Pn['<STR_LIT>'] = P8(x)<EOL>Pn['<STR_LIT>'] = P10(x)<EOL>Pn['<STR_LIT>'] = P12(x)<EOL>Pn['<STR_LIT>'] = P14(x)<EOL>Pn['<STR_LIT>'] = P16(x)<EOL>Pn['<STR_LIT>'] = P18(x)<EOL>Pn['<STR_LIT>'] = P20(x)<EOL>Pn['<STR_LIT>'] = P22(x)<EOL>Pn['<STR_LIT>'] = P24(x)<EOL>Pn['<STR_LIT>'] = P26(x)<EOL>Pn['<STR_LIT>'] = P28(x)<EOL>return Pn<EOL>
Calculate Legendre polyomials P0 to P28 and returns them
    in a dictionary ``Pn``.

    :param float x:     argument to calculate Legendre polynomials
    :return Pn:         dictionary which contains order of Legendre polynomials
                        (from 0 to 28) as keys and the corresponding evaluation
                        of Legendre polynomials as values.
    :rtype:             dict