Запись ядра интегрального уравнения в MATLAB.

[tpom6oh]

Интегральное уравнение решается методом конечных сумм,

Используется синтаксис Matlab M

function y = IntEq( l, c, d, N )
h = (d-c)/(N-1);
for i = 1:N
    x(i) = c+(i-1)*h;
    F(i) = l;
end
A(1) = h/2; A(N) = h/2;
for i = 1:N
    A(i) = h;
end
%k(n, s) = 0.5 .* exp(-log(2).*(x/(1.9585)).^2) .* x; тут выдавало ошибку - Undefined function or variable 'r'.

for i = 1:N
    for j = 1:N
        K(i, j) = k(x(i),x(j));
    end
end
for i = 1:N
    for j = 1:N
          H(i,j) = A(j)*K(i,j);
    end
end
y = H^(-1)*F';
plot(x,y);
title('Решение уравнения расхода газа в струе');
grid on;
end

function k = k(n, s)
k(n, s) = 0.5 .* exp(-log(2).*(n/(1.9585)).^2) .* s;
end
% тут ошибка - Attempted to access k(0,0); index must be a positive integer or logical.
 

Ругается на k(n, s). Как ее правильно записать?