/* as*(t^3)+bs*(t^2) etc */

volatile double spline,as,bs,cs,ds,tmp;

tmp = 0.001; // остальные тоже с вменяемыми значениями

spline = ds + tmp * (cs + tmp * (bs + as*tmp));

//---------------------------------------------------------------------------
#include <windows.h>
#include <winbase.h>
#include <stdio.h>
#include <sys\timeb.h>
#pragma hdrstop

//---------------------------------------------------------------------------
static unsigned long long base_m;
static double base_d=0;
static long long ddm=0;

#pragma argsused
double prec_time(void)
{
LARGE_INTEGER counter_b;
LARGE_INTEGER counter_e;
LARGE_INTEGER freq;
unsigned long long m;
struct timeb ct;
double d;

QueryPerformanceFrequency(&freq);
QueryPerformanceCounter(&counter_b);
ftime(&ct);
QueryPerformanceCounter(&counter_e);
m=(((unsigned long long)counter_b.QuadPart>>1)+((unsigned long long)counter_e.QuadPart>>1));
//m - новая поправка
d=ct.time*1000.0+ct.millitm; //Собственно время в миллисекундах
if (base_d)
{
long long dm;
double dd=d-base_d;
base_m+=(dd)*((double)freq.QuadPart/1000.0);
dm=m-base_m;
ddm+=dm;
if (dd>0)
{
printf("%d\n",ddm);
if (ddm<0)
base_m-=freq.QuadPart/300000;
else
base_m+=freq.QuadPart/300000;
ddm=0;
}
base_d=d;
return d+((double)dm/((double)freq.QuadPart/1000.0));
}
else
{
base_m=m;
base_d=d;
return d;
}
}


int main(int argc, char* argv[])
{
int i;
struct timeb ct;
double pt=prec_time();

for(i=0;i<10000;i++)
{
double t1;
double t2;
ftime(&ct);
t1=ct.time*1000.0+ct.millitm;
t2=prec_time();
printf("%f - %f = %f, interval = %f\n",t1,t2,t2-t1,t2-pt);
pt=t2;
Sleep(2);
}
return 0;
}