Моделирование полета тела с учетом сопротивления

#include <iostream>
#include <windows.h>
#include <math.h>
#include "resource.h"
using namespace std;
 
int x1, y11, x2, y2; 
int x3[474], y3[474];
int x4[474], y4[474];
int WINAPI DlgProc(HWND hDlg, WORD wMsg, WORD wParam, DWORD)
{
    PAINTSTRUCT ps;
    if (wMsg == WM_CLOSE || wMsg == WM_COMMAND && wParam == IDOK) {
        EndDialog(hDlg, 0);
    }
    else
        if (wMsg == WM_INITDIALOG)
        {
            
            RECT rc;
            GetClientRect(hDlg, &rc);
            int dx = rc.right - rc.left;
            int dy = rc.bottom - rc.top;
            cout << "dx=" << dx << " dy=" << dy << endl;
        }
        else
            if (wMsg == WM_PAINT) {
                BeginPaint(hDlg, &ps);
                HPEN hPen = (HPEN)CreatePen(PS_SOLID, 2, RGB(0, 255, 0));
                HPEN hOldPen = (HPEN)SelectObject(ps.hdc, hPen);
                POINT ptOld;
                hPen = (HPEN)CreatePen(PS_SOLID, 2, RGB(0, 0, 0));
                hOldPen = (HPEN)SelectObject(ps.hdc, hPen);
                MoveToEx(ps.hdc, 1, 296, &ptOld);
                for (int i = 0; i < 474; i++)
                    LineTo(ps.hdc, x4[i], y4[i]);
                SelectObject(ps.hdc, hOldPen);
                DeleteObject(hPen);
                EndPaint(hDlg, &ps);
            }
    return 0;
}
void main()
{
    x1 = 1;  y11 = 1;
    x2 = 474;  y2 = 296;
    cout << "x1 = " << x1 << "\ny11 = " << y11 << "\nx2 = " << x2 <<
        "\ny2 = " << y2 << "\n" << flush;
    double vx, vy;
    double t, ycoord, dt;
    cout << "input velocity: vx,vy:";
    cin >> vx >> vy;
    double k = 0.032, m = 50;
    ycoord = 0;
    t = 0;
    x4[0] = 0;
    y4[0] = 296;
 
    for (int i = 1; i < 474; i++)
    {
        dt = i - (i - 1); //решил брать такой промежуток
        x4[i] = x4[i - 1] + vx * dt;
        ycoord = ycoord + vy * dt;
        vy = vy + 10 * dt + (k / m) * sqrt(vy * vy + vx * vx) * vy * dt;
        vx = vx + (k / m) * vx * sqrt(vy * vy + vx * vx) * dt;
        y4[i] = 296 - ycoord;
    }
 
    DialogBox(NULL, MAKEINTRESOURCE(IDD_DIALOG1), NULL, (DLGPROC)DlgProc);
}