CODE:

typedef struct {
	double x;
	double y;
	double vx;
	double vy;
	double m;
	double r;
} chara_t;

static void dumpcalc(chara_t* ca,chara_t* cb,double e) {
	/*２物体の中心の角度*/
	double angle;
	/*その角度でのsinとcos*/
	double sin_a,cos_a;
	/*分解した速度*/
	/*s：垂直　h：平行*/
	double ca_s_v,ca_h_v,cb_s_v,cb_h_v;
	/*衝突後の分解した平行方向の速度*/
	double ca_h_v_2,cb_h_v_2;
	/*角度を求める*/
	angle=atan2(cb->y-ca->y,ca->x-cb->x);
	/*sinとcosを求めておく*/
	sin_a=sin(angle);
	cos_a=cos(angle);
	/*速度を分解する*/
	/*cos(pi/2-angle)==sin(angle)*/
	/*sin(pi/2-angle)==cos(angle)*/
	ca_s_v=ca->vx*sin_a+ca->vy*cos_a;
	ca_h_v=ca->vx*cos_a-ca->vy*sin_a;
	cb_s_v=cb->vx*sin_a+cb->vy*cos_a;
	cb_h_v=cb->vx*cos_a-cb->vy*sin_a;
	/*近づいていなければキャンセル*/
	if(cb_h_v-ca_h_vm+cb_h_v*cb->m-cb->m*e*(ca_h_v-cb_h_v))/(ca->m+cb->m);
	cb_h_v_2=(ca_h_v*ca->m+cb_h_v*cb->m+ca->m*e*(ca_h_v-cb_h_v))/(ca->m+cb->m);
	/*分解した速度を元に戻す*/
	ca->vx=ca_s_v*sin_a+ca_h_v_2*cos_a;
	ca->vy=ca_s_v*cos_a-ca_h_v_2*sin_a;
	cb->vx=cb_s_v*sin_a+cb_h_v_2*cos_a;
	cb->vy=cb_s_v*cos_a-cb_h_v_2*sin_a;
}