{
smat(k1,j) = smat(k1,j)-t*smat(i,j);
}
}
}
else
{
result = false;
break;
}
}
if( result )
{
i = n;
do
{
y(i) = smat(i,n+1);
j = i+1;
while(j<=n)
{
y(i) = y(i)-smat(i,j)*y(j);
j = j+1;
}
y(i) = y(i)/smat(i,i);
i = i-1;
}
while(i>=1);
}
return result;
}
Данная функция решает интегральное уравнение Фредгольма второго рода, заданное ядром интегрирования K(X,S) и правой частью F(X), на отрезке [A, B] методом итераций.
Результат помещается в массив Y с номерами элементов от 1 до N, где 1 соответствует A, N соответствует B.
Epsilon - малое число, передаваемое для сравнения с нолем в ходе решения получаемой системы уравнений.
Для работы этой функции необходима библиотека ap.h
В заключение данной курсовой хотелось бы отметить, что был составлен алгоритм, и на его основе написана функция для решения линейных интегральных уравнений методом итераций. Эта функция может стать основой для написания целой системы, которая будет решать задачи нахождения решения линейных интегралов.
1. Коробов, Н. М. Теоретико-числовые методы в приближенном анализе/ Н. М. Коробов. –М.: 2003. - 316 с.
2. Коробов, Н. М. О приближенном решении интегральных уравнений/
Н. М. Коробов. –ДАН СССР, 1959.
3. http://alglib.sources.ru