ばねの運動が上手くいかないです。
rx,ry,vx,vy,ax,ayを位置、速度、加速度ベクトル、
cx,cyをばねが自然長になった時の位置としました。

case WM_TIMER:
rx=(rx+vx)%width;ry=(ry+vy)%height;
vx+=ax;vy+=ay;
ax= (cx-rx)/50;ay= (cy-ry)/50;
InvalidateRect(hw,NULL,TRUE);
break;

こうすると、単振動の振幅が増大していってしまうんです。
どうやったら解決できますか？