Решение системы уравнений методом Гаусса. Исправить код

0,3  0,1  2 | 2
1   3   0 |1
0.9  0.3  6 |0

#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#pragma warning(disable: 4996)
using namespace std;
void vvodsohr(); void chtreshsohr(); void nevyazka();
int n, i, j;
 
FILE *ish, *rez;
 
int main()
{
    while (true){
        cout << "MENU:" << endl << "1. Zapis' matricy v file." << endl << "2. Reshenie matricy." << endl << "3. Nevyazka." << endl << "0. Vyhod." << endl;
        cin >> i;
        switch (i){
        case 1: vvodsohr(); break;
        case 2: chtreshsohr(); break;
        case 3: nevyazka(); break;
        case 0: return 0;
        }
    }
}
 
void vvodsohr(){
 
    ish = fopen("ishodnye_dannye.txt", "wb");
    cout << "Vvedite razmer matricy (nxn): ";
    cin >> n;
    fwrite(&n, sizeof(int), 1, ish);
    float **a = new float*[n], *b = new float[n]; //инициализация динамических массивов
    for (i = 0; i < n; i++) a[i] = new float[n];
    for (i = 0; i<n; i++){//заполнение массива и запись в файл
        cout << "Vvodite elementy i svobodny chlen stroki " << i + 1 << ":";
        for (j = 0; j < n; j++){
            cin >> a[i][j];
            fwrite(&a[i][j], sizeof(float), 1, ish);
        }
        cin >> b[i];
        fwrite(&b[i], sizeof(float), 1, ish);
    }
 
    cout << endl << "File sozdan, vvodimye dannye sohraneny." << endl << endl;
    fclose(ish);
    for (i = 0; i < n; i++) delete[]a[i];
    delete[]a;
    delete[]b;
}
 
void chtreshsohr(){
    ish = fopen("ishodnye_dannye.txt", "rb");
    fread(&n, sizeof(float), 1, ish);
    float **a = new float*[n], *b = new float[n];
    for (i = 0; i < n; i++) a[i] = new float[n];
 
    for (i = 0; i < n; i++){//считывание из файла в массив.
        for (j = 0; j < n; j++)
            fread(&a[i][j], sizeof(float), 1, ish);
        fread(&b[i], sizeof(float), 1, ish);
    }
 
 
 
 
    int f, jmax; float k, max, *p = NULL, l; //приводим к стандартному виду
    for (i = 0; i < n; i++){
        max = a[i][i];
        jmax = i;
        for (j = i; j < n; j++)
        if ((a[j][i])>max) { max = a[j][i]; jmax = j; }
 
        p = a[i];
        a[i] = a[jmax];
        a[jmax] = p;
        l = b[i];
        b[i] = b[jmax];
        b[jmax] = l;
 
        for (f = i + 1; f < n; f++){
            k = a[f][i] / a[i][i];
            b[f] = b[f] - k*b[i];
            for (j = i; j < n; j++)
                a[f][j] = a[f][j] - a[i][j] * k;
        }
    }
 
    cout << "Standartnaya:" << endl;
    for (i = 0; i < n; i++){//вывод приведенной матрицы
        for (j = 0; j < n; j++)cout << setw(9) << a[i][j] << " ";
        cout << "|" << b[i] << endl;
    }
 
    for (i = 0; i < n; i++)
    if (a[i][i] == 0){
        cout << "Sistema ne imeet resheniy" << endl;
        cout << "Press any key to exit...";
        getch();
        exit(1);
    }
 
    float s = 0, *x = new float[n];//находим x
    for (i = 0; i<n; i++) x[i] = 0;
    for (i = n - 1; i >= 0; i--){
        s = b[i];
        for (j = n - 1; j>i; j--) s -= x[j] * a[i][j];
        x[i] = s / a[i][i];
    }
 
    for (i = 0; i < n; i++) cout << "x" << i + 1 << "=" << x[i] << " ";
    cout << endl << endl;
    fclose(ish);
    rez = fopen("rezult.txt", "wb");//запись иксов в файл
    for (i = 0; i < n; i++)
        fwrite(&x[i], sizeof(float), 1, rez);
    fclose(rez);
    for (i = 0; i < n; i++) delete[]a[i];
    delete[]a;
    delete[]b;
    delete[]x;
}