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.
