var dist += 1 / curveLength * speed * time_step

spline(ent.x, dist);

spline(nextpoint.x,dist);

while(vec_dist(ent.x,nextpoint.x) > 0)
{
  c_move(...); // speed * time_step
}

dist += 0.1

vec_set(dir,currentNode->next->pos);
vec_sub(dir,currentNode->pos);
vec_normalize(dir,speed*time_step);
while(distToNext>=speed*time_step)
{
	//Interpolate the rotation(Not needed if you don't care about the rotation, cause the c_move is done in absoulute coordination
	var percentDone=1-(distToNext/wholeDist);
	quat_slerp(&quaternInter,&quaternStart,&quaternEnd,percentDone);
	ang_for_quat(my.pan,&quaternInter);
	c_move(me,nullvector,dir,IGNORE_MODELS);
	distToNext=vec_dist(my.x,currentNode->next->pos);
	wait(1);
} 

function quat_slerp(QUATERNION *q, QUATERNION *u, QUATERNION *v, double t) {
	QUATERNION qa, qb;
	quat_set(&qa, u);
	quat_set(&qb, v);
	double cosHalfTheta=qa.w*qb.w + qa.x*qb.x + qa.y*qb.y + qa.z*qb.z;
	if(abs(cosHalfTheta)>=1.0)
		quat_set(&q,&qa);
	else
	{
		double halfTheta=acos(cosHalfTheta);
		double sinHalfTheta=sqrt(1.0 - cosHalfTheta*cosHalfTheta);
		if (abs(sinHalfTheta) < 0.001)
		{ 
			q.w = (qa.w * 0.5 + qb.w * 0.5);
			q.x = (qa.x * 0.5 + qb.x * 0.5);
			q.y = (qa.y * 0.5 + qb.y * 0.5);
			q.z = (qa.z * 0.5 + qb.z * 0.5);
		}
		else
		{
			double ratioA = sin((1 - t) * halfTheta) / sinHalfTheta;
			double ratioB = sin(t * halfTheta) / sinHalfTheta; 
			//calculate Quaternion.
			q.w = (qa.w * ratioA + qb.w * ratioB);
			q.x = (qa.x * ratioA + qb.x * ratioB);
			q.y = (qa.y * ratioA + qb.y * ratioB);
			q.z = (qa.z * ratioA + qb.z * ratioB);
		}
	}
}