#include <iostream>
#include <stdio.h>
#include <math.h>
#include <float.h>
using namespace std;

#define DT 0.001 // 時刻の刻み幅
#define TIMELIMIT  20.0//シミュレーションの打ち切り時刻
#define Rm 0.1 //距離rの最低値
#define N 256 //配置できる電荷数の最大値

struct coordinate
{
    double x;
    double y;
};

struct charge
{
    struct coordinate qxy; //座標の位置
    double q; //電荷
};

inline double getdouble(double x)//実数読み込み
{   
    cin>>x;
    
    return x;
}
int inputq (struct charge qi[])
{int i;
    for(i=0; i<N; i++)
    {   //cout << i<<"番目";
        //cout <<"配置位置x_iを入力してください";
        qi[i].qxy.x=getdouble(0);
        //cout <<"配置位置y_iを入力してください";
        qi[i].qxy.y=getdouble(0);
        //cout <<"電荷の値q_iを入力してください";
        qi[i].q=getdouble(0);
    }
    return i;
}

/*****************************************/
/*        main()関数                     */
/****************************************/

int main()
{
    struct coordinate v,r; //質点の速度と位置
    struct charge qi[N]; //電荷qi
    int nofq; //電荷の個数
    double t=0; //時刻 
    double h=DT; //時刻の刻み幅
    double Rx, Ry, R, Rmin ; //電荷と質点の距離

    cout <<"初期速度v_0xを入力してください\n";
    v.x=getdouble(0);
    cout <<"初期速度v_0yを入力してください\n";
    v.y=getdouble(0);
    cout <<"初期位置x_0を入力してください\n";
    r.x=getdouble(0);
    cout <<"初期位置y_0を入力してください\n";
    r.y=getdouble(0);
    nofq=inputq(qi) ;//電荷の入力

    cout<<"t,x,y,v_x,v_y\n";
    cout<< t <<","<< r.x <<","<<r.y <<","<<v.x <<","<< v.y<< "\n";

    //二次元運動の計算
    while(t<=TIMELIMIT) //打ち切り時間までの計算
    {
        t+=h;
        Rmin=100 ; //距離の最小値を初期化
        for(int i=0; i<nofq; i++)
        {
            Rx=qi[i].qxy.x-r.x; //距離の計算
            Ry=qi[i].qxy.y-r.y;
            R=sqrt(Rx/R/R/R*qi[i].q)*h;
            R=sqrt(Rx*Rx+Ry*Ry);
            if(R<Rmin) Rmin=R; // 距離の最小値を更新
            v.x+=(Rx/R/R/R*qi[i].q)*h ;
            v.y+=(Ry/R/R/R*qi[i].q)*h ;
        }
        r.x+=v.x*h ;
        r.y+=v.y*h ;
        cout<< t <<","<< r.x <<","<<r.y <<","<<v.x <<","<< v.y<< "\n";
        if(Rmin<Rm) break; //電荷に非常に近づいたら終了; 
    }
    return 0;
}