Разложение фунции, Разложение функции в ряд Гауссианов Опции

[MaxPIC]

Как можно выполнить разложение функции на Гауссы в заданном диапазоне [a,b].
Функция Гаусса Gauss(x,m,s)=1/(sqrt(2*Pi)*s)*exp(-(x-m)^2/(2*s^2)).
Ряд должен иметь следующую форму:
F(x)=A00+sum(sum(Aij*Gauss(x,m0+dm*i,s0+ds*j),i=0..N),j=0..M),
Причём m0,dm,N,s0,ds,M заданы и конечны. Aij - амплитуды, которые необходимо найти. Как можно выполнить такое разложение?

[gena_dj]

i=0..N),j=0..M)...m0,dm,N,s0,ds,M заданы и конечны

Так и не в ряд, судя по заданию разложить надо, а апроксимировать функцию суммой Гауссоид. Метод наименьших квадратов пробовали?