перевод кода с C# на C++

/// <summary>
    /// Description of MainForm.
    /// </summary>
    public partial class MainForm : Form
    {
        int N; // кол-во измерений
        double a,b; // диапазон
        double h=1.0; // амплитуда случайной погрешности
        double [] x; // аргументы
        double [] y; // результаты
        Random rand=new Random();
        double k,d; // к-ты линейной функции
        
        // исходная функция
        double F(double x)
        {
            return -5.0*Math.Log(x)+2.0;
        }
        
        // линейная функция - аппроксимация
        double f(double x)
        {
            return k*x+d;
        }
        
        public MainForm()
        {
            //
            // The InitializeComponent() call is required for Windows Forms designer support.
            //
            InitializeComponent();
            
            //
            // TODO: Add constructor code after the InitializeComponent() call.
            //
            Calc();
        }
        
        void Calc()
        {
            N=Convert.ToInt32(edtN.Text);
            a=Convert.ToDouble(edtA.Text);
            b=Convert.ToDouble(edtB.Text);
            h=Convert.ToDouble(edtH.Text);
 
            // Инициализируем массивы X и Y
            x=new Double[N];
            y=new Double[N];
            for(int i=0;i<N;++i) {
                x[i]=a+(b-a)*rand.NextDouble();
                y[i]=F(x[i])+(-1.0+2.0*rand.NextDouble())*h;
            }
            
            // Ищем линейную зависимость f(x)=k*x+d, т.ч. sum (y[i]-(k*x[i]+d))^2 -> min
            // производная по d: sum((y[i]-k*x[i])+d)=0
            // k*sum(x[i])-N*d=sum(y[i])
            // производная по k: sum((d-y[i])*x[i]+k*x[i]^2)=0
            // k*sum(x[i]^2)+d*sum(x[i])=sum(x[i]*y[i])
            // Решаем СЛУ 
            double a11=0,a12=-N,a21=0,a22,b1=0,b2=0;
            for(int i=0;i<N;++i) {
                a11+=x[i];
                b1+=y[i];
                a21+=x[i]*x[i];
                b2+=x[i]*y[i];
            }
            a22=a11;
            double D=a11*a22-a12*a21;
            k=(b1*a22-a12*b2)/D;
            d=(a11*b2-b1*a21)/D;
            
            double e=0;
            for(int i=0;i<N;++i) {
                double ee=k*x[i]+d-y[i];
                e+=ee*ee;
            }
            e/=N;
            
            txtStat.Text="Для функции -5·ln(x)+2 по "+N.ToString()+" случайным точкам в диапазоне [2,20] построенно линейное приближение\r\n"+
                k.ToString()+"·x+"+d.ToString()+"\r\n"+
                "Среднеквадратичная ошибка составляет "+e.ToString();
                
        }
        
        void Panel1Paint(object sender, PaintEventArgs e)
        {
            Graphics g=e.Graphics;
            
            int r=4; // радиус "точек"
            float w=panel1.Width,h=panel1.Height;
            
            double maxy=Math.Abs(y[0]);
            for(int i=1;i<N;++i)
                maxy=Math.Max(maxy,Math.Abs(y[i]));
            double lx=(b-a)*1.2,ly=maxy*1.5;
            float kx=w/(float)lx,ky=0.5f*h/(float)ly;
            float sx=0.1f*w,sy=0.5f*h;
            float dx=(float)Math.Exp(Math.Log(10)*Math.Floor(Math.Log10(0.33*lx)));
            float dy=(float)Math.Exp(Math.Log(10)*Math.Floor(Math.Log10(0.33*ly)));
            
            // рисуем оси координат
            Pen pen0=new Pen(Color.Black,3);
            g.DrawLine(pen0,0,sy,w,sy);
            Brush brush0=new SolidBrush(Color.Black);
            Font font0=new Font("Times New Roman",8);
            // OX
            for(double x1=(a-sx/kx),x2=x1+(double)w/kx,x=Math.Ceiling(x1/dx)*dx;x<=x2;x+=dx) {
                float xx=(float)(x-a)*kx+sx;
                g.DrawLine(pen0,xx,sy-5,xx,sy+5);
                SizeF size=g.MeasureString(x.ToString(),font0);
                g.DrawString(x.ToString(),font0,brush0,xx-size.Width*0.5f,sy+size.Height*0.3f);
            }
            // OY
            for(double y2=(double)h/ky,y1=dy;y1<=y2;y1+=dy) {
                float yy=(float)y1*ky;
                g.DrawLine(pen0,0,yy,5,yy);
                g.DrawLine(pen0,w-5,yy,w,yy);
                SizeF size=g.MeasureString(y1.ToString(),font0);
                g.DrawString(y1.ToString(),font0,brush0,8,-yy+sy-size.Height*0.5f);
                g.DrawString(y1.ToString(),font0,brush0,w-8-size.Width,-yy+sy-size.Height*0.5f);
                g.DrawString((-y1).ToString(),font0,brush0,8,yy+sy-size.Height*0.5f);
                g.DrawString((-y1).ToString(),font0,brush0,w-8-size.Width,yy+sy-size.Height*0.5f);
            }
            
            // сортируем X, это надо для правильной отрисовки линий графика
            for(int i=0;i<N;++i) {
                for(int j=i+1;j<N;++j) {
                    if(x[j]<x[i]) { // обмен i и j
                        double t=x[i];x[i]=x[j];x[j]=t;
                        t=y[i];y[i]=y[j];y[j]=t;
                    }
                }
            }
 
            
            // линии графика
            Pen pen1=new Pen(Color.Green,3),pen2=new Pen(Color.Red,3);
            for(int i=1;i<N;++i) {
                g.DrawLine(pen1,(float)(x[i-1]-a)*kx+sx,(float)-y[i-1]*ky+sy,(float)(x[i]-a)*kx+sx,(float)-y[i]*ky+sy);
                g.DrawLine(pen2,(float)(x[i-1]-a)*kx+sx,(float)-f(x[i-1])*ky+sy,(float)(x[i]-a)*kx+sx,(float)-f(x[i])*ky+sy);
            }
 
            // точки на графике
            Brush brush1=new SolidBrush(Color.Green),brush2=new SolidBrush(Color.Red);
            for(int i=0;i<N;++i) {
                g.FillPie(brush1,(float)(x[i]-a)*kx+sx-r,(float)-y[i]*ky+sy-r,2*r,2*r,0,360);
                g.FillPie(brush2,(float)(x[i]-a)*kx+sx-r,(float)-f(x[i])*ky+sy-r,2*r,2*r,0,360);
            }
        }
        
        void Panel1Resize(object sender, EventArgs e)
        {
            panel1.Refresh();
        }
        
        void Button1Click(object sender, EventArgs e)
        {
            Calc();
            panel1.Refresh();
        }
    }
}