#include <stdlib.h>
#include <math.h>
#include "real.h"
#define frac(x) x-floor(x)

static int j_quasi=0;

void quasi(int n,real *a,real *b, real *x) {

int i; static double *z;

if(j_quasi==0){
	z=(double*)malloc(n*sizeof(double));
	for(i=0;i<n;i++) z[i]=frac(cbrt((double)(i+1.234567890987654321)));
	}

j_quasi++; for(i=0;i<n;i++) x[i]=a[i]+frac(j_quasi*z[i])*(b[i]-a[i]);

}