GS Interpolaciq Kriva na Bezie

#include <iostream>
#include <graphics.h>
#include <string>
using namespace std;

int main()
{
    int Dx(50),Dy(80),Px(500),Py(400);
    double x[10] = {1,2,3,4,5,6,7,8,9,10},y[10]={1,2,-3,4,6,5,7,13,9,1};
    double xmin = x[0], ymin = y[0],xmax = x[0],ymax = y[0];
    for (int i=1;i<10;i++)
    {
        if (xmin > x[i]) xmin = x[i];
        if (ymin > y[i]) ymin = y[i];
        if (xmax < x[i]) xmax = x[i];
        if (ymax < y[i]) ymax = y[i];
    }
    double Sx = (xmax - xmin )/ Px, Sy = (ymax - ymin) / Py;
    int ki = Px/Dx,G = Py/Dy;

    initwindow(800,600);
    line (200,450,200+Px,450);
    line (200,450,200,450-Py);
    moveto (200,450);
    for (int i=0;i<ki;i++)
    {
        stringstream strs;
        strs << x[i];
        string temp_str = strs.str();
        char* char_type = (char*) temp_str.c_str();
        outtextxy(200+i*Dx,470, char_type);
        line (200+i*Dx,450,200+i*Dx,450+5);
    }
    for (int i=0;i<G;i++)
    {
        stringstream strs;
        strs << y[i];
        string temp_str = strs.str();
        char* char_type = (char*) temp_str.c_str();
        outtextxy(180,450-i*Dy, char_type);
        line (200,450-i*Dy,200-5,450-i*Dy);

    }
    int m=10,n;
    for (int i=0;i<10;i++)
    {
        circle (200+((x[i]-xmin)/Sx),450-((y[i]-ymin)/Sy),2);
        for (float t=0.0f; t<=1.0f;t+=0.007)                //dt = 0.0007
        {


            for (int j=0; j<m;j++)
            {
                                        //Ri = Pi
                n = m;
                while (n > 0)
                {
                    for (int s=0;s<n-1;s++)
                    {

                    }
                    n -= 1;
                    for (int s=0;s<n;s++)
                    {

                    }
                }
            }
        }

       //lineto(200+((x[i]-xmin)/Sx),450-((y[i]-ymin)/Sy));
    }
    getch();
    return 0;
}