Нужно программу на Паскале переделать в программу на С++

#include <iostream>
using namespace std;
 
int main()
{
    int i, i1, i2 ,N, z;  // var i, i1, i2 ,N, z: integer; 
    double q,delta,X,Y,a,b,k,Y1,X1; // q,delta,X,Y,a,b,k,Y1,X1: real; 
    double M[2][20], M1[2][20]; // M, M1: array[1..2, 1..20] of real;
    // begin 
    cout << "Input tha count of dots N=";//write('Input tha count of dots N=');
    cin >> N;//readln(N);
    for (i = 0; i <= N; i++) // for i:=1 to N do 
    {
        cout << "Input dot coordinates " << i << " X="; // write('Input dot coordinates ',i,' X=');
        cin >> M[0][i]; // readln(M[1,i]);
        cout << "Input dot coordinates " << i << " Y="; // write('Input dot coordinates ',i,' Y=');
        cin >> M[1][i]; //  readln(M[2,i]);
    }
 
    cout << "Input the aperture value q="; // write('Input the aperture value q='); 
    cin >> q; // readln(q);
  
 
    i1 = 1; i2 = 3;  
    // в массив заносим первую точку
    z=1; M1[0][1]=M[0][1]; M1[1][1]=M[1][1];
    while ((i1<=N-2) && (i2<=N)) // while (i1<=N-2) and (i2<=N) do
    { // begin 
       a = (M[1][i2]-M[1][i1])/(M[0][i2]-M[0][i1]);
       b = M[1][i2]-a*M[0][i2];
       cout << "a=" << a << " ; b=" << b;
       Y1=a*M[0][i2-1]+b;
       delta=abs(M[1][i2-1]-Y1);
       if (delta>q) //then 
       { //begin 
           cout << "Y=" << M[1][i2-1] << "; Y1=" << Y1 << "; delta=" << delta << " ne szimaetsya";
           i1 =i1+1; i2=i1+2;
           z=z+1;
           M1[0][z]=M[0][i2-1] ;
           M1[1][z]=M[1][i2-1];
       } // end
       else
       { //begin
           cout << "Y=" << M[1][i2-1] << "; Y1=" << Y1 << "; delta=" << delta; // { здесь как раз и выводятся не вошедшие в сжатый массив точки }
           i2++; // i2=i2+1;
       } //end;
    }//end;   
  //{записываем в конец массива посл точку }
  z = z+1; M1[0][z]=M[0][N] ; M1[1][z]=M[1][N];
 
  cout << "Array M*:";
  for (i = 0; i <= z; i++)
      cout << " X[" << i << "]=" << M1[0][i] << "; Y[" << i << "]=" << M1[1][i];
 
  cout << "Коэффициент сжатия k = " << N/z;
return 0;
}