[zaliczenie] metoda eliminacji Gaussa w C, PILNE! :|

Hej :-(
Bardzo proszę o pomoc w napisaniu programu w języku C#
Ma być to program rozwiązujący układ 5 równań z 5 niewiadomymi metodą eliminacji Gaussa (Dane wejściowe - liczba równań, współczynniki przy poszczególnych niewiadomych oraz wyjściowe - rozwiązanie układu powinny znaleźć się w pliku tekstowym).
Nie jestem w stanie sobie sama z tym poradzić, a termin goni. :-(
Dziękuję...

Rozwiazanie ukladu rownan liniowych A x = b metoda eliminacji Gaussa

A - kwadratowa macierz wspolczynnikow ukladu N x N
b - wektor prawych stron

for (int i = 1; i <= N-1; i++) {
for (int j = i+1; j <= N; j++) {
for (int k = i+1; k <= N; k++ ) {
A[j-1][k-1] -= A[i-1][k-1]*A[j-1][i-1]/A[i-1][i-1];
}
b[j-1] -= b[i-1]*A[j-1][i-1]/A[i-1][i-1];
}
}

b[N-1] /= A[N-1][N-1];
for (int i = 1; i <= N-1; i++) {
double s = 0;
for (int j = N+1-i; j <= N; j++)
s += A[N-i-1][j-1]*b[j-1];
b[N-i-1] = (b[N-i-1] - s)/A[N-i-1][N-i-1];
}

A jak wyglądałby początek z sekcją VAR?

Masz tu kod w BC++ 3.1, nie wiem jaka jest roznica miedzy c#

#include <iostream.h>
#include <fstream.h>
#include <stdlib.h>

#define N 3

void main()
{
fstream plik;
float A[N][N],b[N];
int i,j;
char str[6];
plik.open("d:\wynik.txt",ios::out);
for(i=0; i<N; i++)
{
for(j=0; j<N; j++)
{
cout<<"A["<<i+1<<","<<j+1<<"] = ";
cin>>A[i][j];
gcvt(A[i][j],6,str);
plik.write(str,6);
}
plik.put('\n');
}
for(i=0; i<N; i++)
{
cout<<"b["<<i+1<<"] = ";;
cin>>b[i];
gcvt(b[i],6,str);
plik.write(str,6);
}
plik.put('\n');

for (i=1; i<=N-1; i++)
for (j=i+1; j<=N; j++)
{
for (int k=i+1; k<=N; k++)
A[j-1][k-1] -= A[i-1][k-1]*A[j-1][i-1]/A[i-1][i-1];
b[j-1] -= b[i-1]*A[j-1][i-1]/A[i-1][i-1];
}

b[N-1] /= A[N-1][N-1];
double s=0;
for (i=1; i<=N-1; i++)
{ s=0;
for (j=N+1-i; j<=N; j++)
s += A[N-i-1][j-1]*b[j-1];
b[N-i-1] = (b[N-i-1] - s)/A[N-i-1][N-i-1];
}

for(i=0; i<N; i++)
{
cout<<b[i]<<endl;
gcvt(b[i],6,str);
plik.write(str,6);
}

plik.close();
}

serdecznie dziękuję :)

Potrafi ktoś przełożyć ten kod z C++ na C#? :|