Školní projekt v C/C++ Maticové operace – C / C++ – Fórum – Programujte.com

Dobrý večer,
mam projekt do školy a už nevím jak dal a tímto bych chtěl vás zde požádat o pomoc sedím nad tím dnem a noci a nemůžu se hnout z místa.

Zadání:
Program přijímá na vstupu matici zadanou po řádcích, a to buď z klávesnice, nebo souboru. Vstup z klávesnice je průběžně validován, a průběžně je plněno dvourozměrné dynamicky alokované pole. Při vstupu ze souboru se kontroly (validace) neprovádí. Nad maticí v paměti lze zavolat Gaussovu eliminační metodu (implementujte ji jako funkci, jejímž konstantním vstupem je ukazatel na vstupní pole, výstupem ukazatel na nově vytvořené pole v trojúhelníkovém tvaru s vypočtenými hodnotami nezávislých proměnných, pokud to hodnost rozšířené matice dovolí). Program umožňuje výstup do HTML souboru, obsahující výslednou matici (použijte HTML tagy, jinak se formátem nijak zvlášť nezabývejte), pod maticí budou hodnoty nezávislých proměnných. Dalším možným výstupem bude výpis matice a hodnoty jejího determinantu, pokud to hodnost matice dovolí.

Můj kod:

// Maticove oprace.cpp : Defines the entry point for the console application.
//

//#include "stdafx.h"
#include <stdio.h>
#include <iostream>
#include<math.h>
#include<iomanip>
using namespace std;

float GEM(float *a[],float *b[], int N, int M)
{
	for (int i = 0; i < N; i++)
	{
		for (int j = 0; j < M; j++)
		{
			b[i][j] = a[i][j];
		}
	}

	int max = 0;
	float pom;
	for (int i = 0; i<N; i++)              //nalezení pivota
	{
		for (int j = i + 1; j<N; j++)
        {
            if (fabs(b[j][i])>fabs(b[i][i]))
                max = j;
        }

		for (int k = i; k<M; k++)
		{
			pom = b[i][k];
			b[i][k] = b[max][k];
			b[max][k] = pom;
		}

		for (int j = i + 1; j<N; j++)
        {
            for (int k = N; k >= i; k--)
            {
                b[j][k] = b[j][k] - b[i][k] * b[j][i] / b[i][i];
            }
        }
	}
/************************************************************************/
    float r[N];
    int res = M-1;
    for(int i=0; i<N; i++)
       r[i]=b[i][res];
    int pr = 0;
    if(pr==0)
    {
     for(int j=N-1; j>=0; j--)
     {
       pom=0;
       for(int k=j+1; k<=N; k++)
           pom=pom+b[j][k]*r[k];
       r[j]=(b[j][N]-pom)/b[j][j];
     }
   }
   for (int j = 0; j < M-1; j++)
		{
			printf("x%d = %5.3f \n",j+1, r[j]);
		}
/************************************************************************/
	return **b;
}
void vypismat(float *a[],int N,int M)
{
    for (int i = 0; i<N; i++)
	{//vypisuje prvky, když dojde na konec řadky, odřádkuje
		for (int j = 0; j<M; j++)
		{
			printf(" %5.3f ", a[i][j]);
			if (j == M - 1)
				printf("\n");

			if (j == N-1)
				printf(" | ");
		}
	}
}
void naplmatici(float *a[],int N,int M)
{
    int res;
    printf("Vyplnte matici:\n");
	for (int x = 0; x < N; x++)
	{
		for (int j = 0; j < M; j++)
		{
			scanf("%f", &a[x][j]);
		}

	}
}
int main(int argc)
{
	int N, M;
	printf("zadejte počet řádku matice:");
	scanf("%d", &N);
	printf("zadejte počet sloupcu matice:");
	scanf("%d", &M);
	/*****************************************************************/
	float**pole;
	pole = new float*[N];
	for (int i = 0; i < N; i++)
	{
		pole[i] = new float[M];
	}
	naplmatici(pole,N,M);
	/************************************************************************/
	float**pole1;
	pole1 = new float*[N];
	for (int i = 0; i < N; i++)
	{
		pole1[i] = new float[M];
	}
	GEM(pole, pole1, N, M);
	/************************************************************/
	vypismat(pole,N,M);
	printf("******************************************************\n");
	vypismat(pole1,N,M);

/*************************************************************/


/*************************************************************/
	for (int i = 0; i<N; i++) {
		delete[] pole[i];
		delete[] pole1[i];
	}
	delete[] pole;
	delete[] pole1;

	return 0;
}

1) Chtěl bych, lepe vyřešit Gaussovu eliminační metodu nebo bych jí spíš chtěl vyřešit podle zadání.

2) Chtěl bych taky načítat s csv souboru což absolutně nevím jak sem začátečník.