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

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