Переделать под с++

#include<stdio.h>
#include<conio.h>
#include<math.h>
#include<stdlib.h>
#include<stdafx.h>
 
void Readstr(char *s, int max) // reads a string
{
    int c;
    while ((c = getchar()) != '\n')
    {
        if (max>1)
            *s++ = c, max--;
        if (max>0)
            *s = 0;
    }
}
 
void ReadMatrix(int flag, double *a, double *c, int sz) //reads the matrix
{
    FILE *f; //the matrix file
    int i, j;
 
    if (flag == 1) //if that file exists
    {
        f = fopen("matrix.txt", "rt"); //opens file on reding and check its existence
        if (f == NULL)                 //if file don't exist
        {
            printf("No file is found. \n");
            return;
        }
 
        int ch;
 
        for (i = 0; i<sz; i++) //reads the matrix A and the free elements matrix from the file
        {
            for (j = 0; j<sz + 1; j++)
            {
                if (((ch = fgetc(f)) != EOF) && j<sz)
                    fscanf(f, "%lf", &a[i*sz + j]);
                else if (j == sz)
                    fscanf(f, "%lf", &c[i]);
            }
        }
 
        fclose(f); //closes the file 
    }
 
    else if (flag == 2) //if you enter a matrix yourself
    {
        for (i = 0; i<sz; i++) //reads the matrix A
        {
            for (j = 0; j<sz; j++)
            {
                printf("Enter an A[%i][%i] element  ", i, j);
                scanf("%lf", &a[i*sz + j]);
            }
        }
 
        for (j = 0; j<sz; j++) //reads free elements
        {
            printf("\nEnter an b[%i] element ", j);
            scanf("%lf", &c[j]);
 
        }
    }
}
 
bool CheckII(double *x, double *prev, int sz, double acc, double &q, double r, double &e, double f)
{                                                       // checking roots by ||x[k] - x[k-1]||<acc
    double n = 0;
    double w;
    int i, j;
 
    for (j = 0; j<sz; j++)
    {
        if (fabs((x[j] - prev[j])) >= acc)
            n++;
    }
 
    if (n == 0)                             // ||x[k] - x[k-1]||
    {
        q = fabs((x[0] - prev[0]));
        for (i = 1; i<sz; i++)
        {
            w = fabs((x[i] - prev[i]));
            if (w >= q)
                q = w;
        }
        e = r - f;                  // f(x[k]) - f(x[k-1])
    }
    if (n)
        return false;
    else
        return true;
}
 
double MPS(double *a, double *c, double *x, int sz, double acc, int &p, double &q, double &e)
{                                                       //minimizes function "f"
    double *prev, g, f, z, r;
    double l = 0.7;     // the first value of step
    int i, j, k, t, y;
 
    prev = new double[sz];
 
    for (i = 0; i<sz; i++)
    {
        for (j = 0; j<sz; j++)
        {
            a[i*sz + j] = a[i*sz + j] / 2;  // (1/2)*A
        }
    }
 
    z = 0;
    y = 0;
 
    for (i = 0; i<sz; i++)              // f = (1/2)*<A*x,x> + <b,x>
    {
        g = 0;
        for (j = 0; j<sz; j++)
        {
            g += a[i*sz + j] * x[j];
        }
        g = g*x[i];
        z += g + c[i] * x[i];
    }
 
    do      //makes and count iterations
    {
    yo:     for (k = 0; k<sz; k++)
    {
        prev[k] = x[k];         // saves previous iteration results
    }
 
            t = 0;
 
            for (k = 0; k<sz; k++)
            {
                x[k] -= l;              // f(x[k-1] - l*e)
                f = 0;
                for (i = 0; i<sz; i++)
                {
                    g = 0;
                    for (j = 0; j<sz; j++)
                    {
                        g += a[i*sz + j] * x[j];
                    }
                    g = g*x[i];
                    f += g + c[i] * x[i];
                }
 
                if (!(f<z))             // f(x[k-1] - l*e) >= f(x[k-1])
                {
                    x[k] += 2 * l;          // f(x[k-1] + l*e)
                    f = 0;
                    for (i = 0; i<sz; i++)
                    {
                        g = 0;
                        for (j = 0; j<sz; j++)
                        {
                            g += a[i*sz + j] * x[j];
                        }
                        g = g*x[i];
                        f += g + c[i] * x[i];
                    }
                    if (!(f<z))         // f(x[k-1] + l*e) >= f(x[k-1])
                    {
                        x[k] = prev[k];
                        t++;
                    }
                    else                // f(x[k-1] + l*e) < f(x[k-1])
                    {
                        r = z; z = f;
                    }       // saves previous value of function "f"
                }
                else            // f(x[k-1] - l*e) < f(x[k-1])
                {
                    r = z; z = f;
                }       // saves previous value of function "f"
            }
            if (t == sz)
            {
                l = l / 2;      // the step subdivision
                goto yo;
            }
            p++;            // a number of steps
    } while (!CheckII(x, prev, sz, acc, q, r, e, f));       // ||x[k] - x[k-1]|| < acc
 
    return f;
    delete prev;
}
void main(void)
{
    double *a, *x, *c, acc; //presents arrays and an accuracy
    double f, q, e;
    int p, sz; //presents the size of your matrix and the number of iterations
    int h, i;
    printf("Enter the size of your matrix: "); // asks for the size of your matrix
    scanf("%i", &sz);
    printf("Enter an accuracy: "); // asks for an accuracy
    scanf("%lf", &acc);
 
    a = new double[sz*sz];
    c = new double[sz];
    x = new double[sz];
 
    do      //asks your wish about the entering your matrix
    {
        printf("Press 1 if you want to read your matrix from the file (matrix.txt)\n");
        printf("Press 2 if you want to enter your matrix yourself\n");
        scanf("%i", &h);
        Readstr(NULL, 0);
    } while ((h != 1) && (h != 2));
    for (i = 0; i<sz; i++) //the first approximation
    {
        printf("Enter the first approximation: x[%i] = ", i);
        scanf("%lf", &x[i]);
    }
    ReadMatrix(h, a, c, sz);  //reads the matrix
    p = 0;
    q = 0;
    e = 0;
    f = MPS(a, c, x, sz, acc, p, q, e);
    for (i = 0; i<sz; i++) //shows need roots by 10-numeral accuracy
    {
        printf("\nx[%i] = %.10lf", i, x[i]);
    }
    printf("\nf = %.10lf", f);          //shows value of function "f"
    printf("\n\n\nThe number of steps: %i\n", p);  //shows the number of steps we need
    printf("\n||x[k] - x[k-1]|| = %.10lf\n", q);    //shows ||x[k] - x[k-1]||
    printf("\nf(x[k]) - f(x[k-1]) = %.10lf\n", e);  //shows f(x[k]) - f(x[k-1])
    delete a;
    delete c;
    delete x;
}