Курсовая работа: Решение линейных интегральных уравнений

y.setbounds(1, n);

result = true;

for(i = 1; i <= n; i++)

{

k1 = i;

m1 = fabs(smat(i,i));

for(j = i+1; j <= n; j++)

{

if( m1<fabs(smat(j,i)) )

{

m1 = fabs(smat(j,i));

k1 = j;

}

}

if( fabs(m1)>=epsilon )

{

for(j = i; j <= n+1; j++)

{

t = smat(i,j);

smat(i,j) = smat(k1,j);

smat(k1,j) = t;

}

for(k1 = i+1; k1 <= n; k1++)

{

t = smat(k1,i)/smat(i,i);

smat(k1,i) = 0;

for(j = i+1; j <= n+1; j++)

{

smat(k1,j) = smat(k1,j)-t*smat(i,j);