Как вычислить двойной бесконечный ряд в MATLAB?

---

Я хочу вычислить двойную бесконечную сумму, определяемую этот Я хочу иметь 100 значений summ, то есть по одному для каждого различного значения 'p', каждый раз используя одно конкретное значение 'deltaHvec'. Я имею в виду для первого первого значения 'summ', p=1 и deltaH=deltaHvec(1). Для второго значения 'summ' 'p=2' и deltaH=deltaHvec(2). Проблема в том, что конечный график "суммирования" должен быть равен или меньше величины "дельтахвека". Но я получаю довольно высокие значения "summ" , что означает, что я добавляю дополнительные значения в "summ" с каждым последующим шагом
.

Что я уже пробовал:

x = 45; % для вычисления Омеги и u
l = 300; % для расчета Омеги
Р = 100;
omega= x/l; deltaHvec = sin (linspace(0.3, 2.6, 100));
deltat = 5; % для расчета u
v = 40; % для вычисления u
u = x/sqrt(v*deltat); suma = 0;
для i = 1: p
deltaH=deltaHvec (i);
для m = 1: i
для n = 1: 1000
e1Arg=(2*(Н-1)*у + у*Омега) / (2*омега*функция sqrt(п-м));
e1=erfc(e1Arg);
e2Arg = (2*n*u-u * omega)/(2*omega*sqrt(p-m));
e2 = erfc(e2Arg);
totalerror = e1+e2;
f= (-1)^(n-1) * deltaH * totalerror;

summ (i) = summ (i) + f;
конец
конец

конец

сумма;
plot (1: 100, summ); hold on; plot (1: 100, deltaHvec,'k');

---