Получить график спирали для физической задачи: удержание позитрона в магнитном поле

#include "stdafx.h"
#include <iostream>
#include <cmath> 
#include <iomanip>
#include <fstream>
using namespace std;
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale(LC_ALL, "Russian");
    FILE *fp;
    if ((fp=fopen("cppstudio.txt", "w")) == NULL)
    {
        return 1;
    }
        const double m=5, q=0.5, B=5, E=10, dt=0.005, r=0.05;
    double v, F, Fx, Fy, x, y, vx, vy, ax, ay, t;
    x=0; y=120; vx=30; vy=0; t=0;
    printf("%9s   %9s", "x", "y\n");
    fprintf(fp,"%9s  %9s", "x", "y\n");
    do
    {
        F=B*v*q; v=sqrt(vx*vx+vy*vy);
        Fx=F*vy/v; Fy=-F*vx/v;
        ax=(Fx-r*vx)/m; ay=(Fy-r*vy)/m;
        vx=vx+ax*dt; vy=vy+ay*dt;
        x=x+vx*dt; y=y+vy*dt;
        
        printf("%9.3f %9.3f\n",x,y);
        fprintf(fp,"%9.3f %9.3f\n",x,y);
        t=t+dt;
    }
    while (t<50)
    system("pause");
    return 0;
}

do
    {
        v=sqrt(vx*vx+vy*vy);
        F=B*v*q;