Козенко, Галанина Информатика_ч2

void VVOD (…) // функция ввода исходных данных

…

float FI (…) // функция вычисления значений базовых функций

…

float koeff (…) // функция вычисления коэффициентов системы уравнений

…

int ved ( float A[m][m], int i ) // функция определения ведущего уравнения

{

int j, g,h,k;

float MaxA , temp;

MaxA=0; h = -1;

for (k = i; k < m; k++)

if ( fabs( A[k][i]) >fabs( MaxA ))

{

MaxA = A[k][i]; h = k;

}

if ( h== -1 )

{

printf (" Матрица вырождена \n "); abort() ;

}

return h;

}

void perest ( float A[m][m], float B[m], int i, int IM ) // функцияперестановкиуравнений

// системы

{

float temp; int j; if ( IM != i )

{

for (j = i; j < m; j++)

{

temp = A[i][j] ; A[i][j] = A[IM][j]; A[IM][j] = temp;

}

Рис. 2.5 Фрагмент программной реализации решения задачи (начало)

temp = B[i]; B[i] = B[IM]; B[IM] = temp;

}

}

void rkoef ( float A[m][m], float B[m], int i, int l ) // функцияпересчетакоэффициентов

{

float Q; int j; Q=A[l][i] / A[i][i] ; A[l][i] = 0;

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

A[l][j] = A[l][j] - Q*A[i][j]; B[l] = B[l] - Q*B[i];

}

void prgauss ( float A[m][m], float B[m] ) // функция прямого хода методаГаусса

{

int i, l, IM; float Q ; for (i=0;i<m;i++)

{

IM= ved(A,i); perest(A,B,i,IM); for (l=i+1;l<m;l++)

rkoef(A,B,i,l);

}

}

void obgauss( float A[m][m], float B[m], float C[m] ) // функция обратного хода

// методаГаусса

{

int k, j; float Sum;

C[m-1] = B[m-1] / A[m-1][m-1]; for( k=m - 2; k>=0; k--)

{

Sum = B[k];

for (j=k+1; j<m; j++)

Sum = Sum - A[k][j]*C[j]; C[k] = Sum / A[k][k] ;

}

}

…

Рис. 2.5. Фрагмент программной реализации решения задачи (продолжение)

float approcs (…) // функция вычисления значений аппроксимирующей функции

…

float krappr (…) // функция определения качества аппроксимации

…

void VYVOD (…) // функция вывода результатов

...

Рис. 2.5. Фрагмент программной реализации решения задачи (окончание)

3. ВАРИАНТЫ ЗАДАНИЙ К КУРСОВОЙ РАБОТЕ

Варианты заданий к курсовой работе приведены в табл. 3.1.

Продолжение табл. 3.1