Схема алгоритма программы

Схема алгоритма функции FI (k, x[i]).

Код программы

#include.

#include.

#include.

#include.

#define m 3.

#define n 6.

void Massiv (float *M, char *Name, int x).

{.

int i, j;

printf («Array %s: «, Name);

for (j=0;j.

{.

printf («%5.3f «, M[j]);

}.

printf (««);

}.

void znachenia (float X[n], float Y[n]).

{.

int i;

printf («vvedite: «);

for (i=0;i.

{.

printf («X[%i]=», i+1);

scanf («%f» ,&X[i]);

printf («Y[%i]=», i+1);

scanf («%f» ,&Y[i]);

}.

}.

float F1(int k, float N).

{.

if (k==0) return 1;

if (k==1) return log (N);

else.

return N;

}.

void koefficent (float A[m][m], float B[m], float X[n], float Y[n]).

{.

int i, l, k;

for (i=0;i.

for (l=0;l.

{.

A[i][l] = 0;

for (k=0;k.

A[i][l]=A[i][l] + F1(i, X[k]) * F1(l, X[k]);

}.

for (k=0;k.

{.

B[k] = 0;

for (i=0;i.

B[k]=B[k]+Y[i]*F1(k, X[i]);

}.

}.

int ved (float A[m][m], int i).

{.

int g, h, k;

float MaxA;

h = -1;

MaxA = 0;

for (k=i;k.

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

{.

MaxA = A[k][i];

h = k;

}.

if (h == -1).

{.

printf («Matrix virogdena «);

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.

{.

temp = A[i][j];

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

A[IM][j] = temp;

}.

temp = B[i];

B[i] = B[IM];

B[IM] = temp;

}.

}.

void rkoef (float A[m][m], float B[m], int i, int l).

{.

int j; float Q;

Q=A[l][i] / A[i][i];

A[l][i] = 0;

for (j=i+1;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, IM, l, j;float Q;

for (i=0;i.

{.

IM = ved (A, i);

perest (A, B, i, IM);

for (l=i+1;l.

rkoef (A, B, i, l);

}.

for (i=0;i.

{.

for (j=0;j.

printf («%f «, A[i][j]);

printf (««);

}.

}.

void obrgauss (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.

Sum = Sum — A[k][j]*C[j];

C[k] = Sum / A[k][k];

}.

}.

void apr (float C[m], float X[n], float Y[n], float Y1[n], float D[n], float *Kr).

{.

int i;

*Kr=0;

for (i=0;i.

{.

Y1[i]=C[0]*F1(0,X[i])+C[1]*F1(1,X[i])+C[2]*F1(2,X[i]);

D[i]=fabs (Y[i]-Y1[i]);

*Kr=*Kr+D[i]*D[i];

}.

}.

void krappr (float D[n], float *Dmax, int *IM).

{.

int i;

• *Dmax=D[0];
• *IM=0;

for (i=0;i.

if (fabs (D[i])>fabs (*Dmax)).

{.

• *Dmax=D[i];
• *IM=i;

}.

}.

void vivodvsego (float C[m], float X[n], float Y1[n], float D[n], float Kr, float Dmax, int IM).

{.

Massiv (C, «C», m);

Massiv (Y1, «Y1», n);

Massiv (D, «D», n);

printf («Dmax = %5.3f pri X[%i] = %5.3f «, Dmax, IM+1,X[IM]);

printf («Kr = %5.4f «, Kr);

}.

int main ().

{.

float X[n], Y[n], A[m][m], B[m], C[m], Y1[n], D[n], Dmax, Q, Kr, Sum;

int i, j, IM;

znachenia (X, Y);

koefficent (A, B, X, Y);

printf («Massiv A «);

for (i=0; i.

{.

for (j=0; j.

printf («%f «, A[i][j]);

printf (««);

}.

Massiv (B, «B», m);

prgauss (A, B);

obrgauss (A, B, C);

apr (C, X, Y, Y1, D, &Kr);

krappr (D, &Dmax, &IM);

vivodvsego (C, X, Y1, D, Kr, Dmax, IM);

getch ();

return 0;

}.