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

int j;

int u;

int k1;

int m;

smat.setbounds(1, n, 1, n+1);

y.setbounds(1, n);

h = (b-a)/(n-1);

i = 1;

do

{

x = a+(i-1)*h;

smat(i,n+1) = f(x);

j = 1;

do

{

smat(i,j) = -h*k(x, a+(j-1)*h);

if( j==1||j==n )

{

smat(i,j) = smat(i,j)/2;

}

if( j==i )

{

smat(i,j) = 1+smat(i,j);

}

j = j+1;

}

while(j<=n);

i = i+1;

}