Важно, Нужна срочна бло-схема

[Олежка9]

ПОМОГИТЕ написать блок-схему по коду!!!

#include <stdio.h>
#include <math.h>
#include <iostream.h>
#include <conio.h>

float *a, *b, *x, *y, **sums;
//a - неизвестные коэффициенты полинома 
//b - столбец свободных членов (в правой части уравнений) 
//x - координаты, заданные в файле 
//y - координаты, заданные в файле
//sums - суммы степеней x,y при неизвестных коэффициентах полинома
int N, K;
//N - число заданных точек в файле 
//K - степень аппроксимирующего полинома
char filename[100];
FILE* InFile=NULL;

void count_num_lines(){// подсчитывает количество точек, где задана функция.
int sparco=0;        
   do{
       sparco = 0;
       while(fgetc(InFile)!='\n' && !feof(InFile)) sparco=1;
       if(sparco) N++;
   }while(!feof(InFile));
}


void freematrix(){ //освобождает память, которая была выделена ранее.
int i;
   for(i=0; i<K+1; i++){
       delete [] sums[i];
   }
   delete [] a;
   delete [] b;
   delete [] x;
   delete [] y;
   delete [] sums;
}

void allocmatrix(){  //выделяет память для массивов a, b, x, y, sums
int i,j,k;			
   a = new float[K+1];
   b = new float[K+1];
   x = new float[N];
   y = new float[N];
   sums = new float*[K+1];
   if(x==NULL || y==NULL || a==NULL || sums==NULL){
   		cout<<"Ne dostatochno pamyati dlya rasmesheniya: N="<<N<<" , K="<<K<<endl;
       return;
   }
   for(i=0; i<K+1; i++){
       sums[i] = new float[K+1];
       if(sums[i]==NULL){
       cout<<"Ne dostatocno pamyati dlya razmesheniya "<<K+1<<" yravnenii"<<endl;	
       return;
       }
   }
   for(i=0; i<K+1; i++){
       a[i]=0;
       b[i]=0;
       for(j=0; j<K+1; j++){
	   sums[i][j] = 0;
       }
   }
   for(k=0; k<N; k++){
       x[k]=0;
       y[k]=0;
   }
}

void readmatrix(){ //прочитывает из файла координаты точек и вычисляет sumsij
   int i=0,j=0, k=0;
for(k=0; k<N; k++){   
       fscanf(InFile, "%f", &x[k]);
       fscanf(InFile, "%f", &y[k]);
   }
for(i=0; i<K+1; i++){
       for(j=0; j<K+1; j++){
	   sums[i][j] = 0;  //второе обнуление
	   for(k=0; k<N; k++){
	       sums[i][j] += pow(x[k], i+j);
	   }
       }
   }
for(i=0; i<K+1; i++){
       for(k=0; k<N; k++){
	   b[i] += pow(x[k], i) * y[k];
       }
   }
}

void printresult(){//распечатывает получившийся столбец решений
int i=0;
   printf("\n");
   for(i=0; i<K+1; i++){
       printf("a[%d] = %f\n", i, a[i]);
   }
}
void diagonal(){  //делает так, чтобы на главной диагонали не было нулей, чтобы не пришлось
// однажды делить на ноль в процессе приведения системы к треугольному виду.
   int i, j, k;
   float temp=0;
   for(i=0; i<K+1; i++){
       if(sums[i][i]==0){
	   for(j=0; j<K+1; j++){
	       if(j==i) continue;
	       if(sums[j][i] !=0 && sums[i][j]!=0){
		   for(k=0; k<K+1; k++){
		       temp = sums[j][k];
		       sums[j][k] = sums[i][k];
		       sums[i][k] = temp;
		   }
		   temp = b[j];
		   b[j] = b[i];
		   b[i] = temp;
		   break;
	       }
	   }
       }
   }
}

int main(){
   int i=0,j=0, k=0;
   do{
       cout<<"Vvedite pyt k failu s dannimi(Primer C:\\sparco.txt):"<<endl;
       cin>>filename;
       InFile = fopen(filename, "rt");
   }while(InFile==NULL);
   count_num_lines(); 
   cout<<"Kolichestvo tochek(strok) v faile: N="<<N<<endl;
   do{
       cout<<"Vvedite kolichestvo aproksimacii polinoma K<N: ";
       cin>>K;
   }while(K>=N);
   allocmatrix(); 
   rewind(InFile); //Перемещение курсора в начало файла 
readmatrix(); 
diagonal(); 
   fclose(InFile);
for(k=0; k<K+1; k++){
       for(i=k+1; i<K+1; i++){
	   if(sums[k][k]==0){
	   		cout<<"Reshenie ne syshestvyet"<<endl;
	       return 0;
	   }
	   float M = sums[i][k] / sums[k][k];
	   for(j=k; j<K+1; j++){
	       sums[i][j] -= M * sums[k][j];
	   }
	   b[i] -= M*b[k];
       }
   }
for(i=(K+1)-1; i>=0; i--){
       float s = 0;
       for(j = i; j<K+1; j++){
	   s = s + sums[i][j]*a[j];
       }
       a[i] = (b[i] - s) / sums[i][i];
   }
   printresult();
   freematrix();
   getch();// ожидание нажатия клавиши
   return 1;
}

[ACE Valery]

Ну если для вас важно, то могли бы и написать, что должна делать программа.

[Олежка9]

Система для аппроксимации линейной зависимости методом наименьших квадратов