The analytical model [Moreno and Sun, Opt. Express 16, 1808-1819 (2008)] is applied to LED with attached secondary optics. It is shown that a slightly modified model using only three cosine power functions can be used for the cases with symmetric radiation pattern and that good fitting on realistic examples can be achieved with several standard optimization algorithms including local search heuristics and genetic algorithm.