The AnalyticalApproximations package is powerful Mathematica® toolbox which allows you to find analytical and numerical polynomial approximation for solutions of functions, initial and boundary value problems with symbolic conditions on symbolic approximation interval with singular regular points. The methods applied are numerically - analytical (a-methods by V. K. Dzyadyk). It means that LDE coefficients, boundary or initial conditions and interval of the approximation can be either symbolical or numerical expressions. These methods can be applied to find approximation of some functions with parameters on symbolic approximation interval approximate polynomial for solution of initial value problems (IVP) approximate polynomial for solution of boundary value problems (BVP) approximate polynomial for solution of LDE with regular singular points (RSP) The methods give asymptotically best approximation in Chebyshev metric.