ブラウン運動を非整数階微積分するプログラムを作成したのですが、上手く動きません。どこがわるいのでしょうか？？
実行結果は数値を出し、エクセルでグラフ化しています。（window.hが使えないため）
#include <stdio.h>　　#include <stdlib.h>
#include <time.h>　　#include <math.h>
#define KMAX 100　　#define XMAX 300
#define PI 3.14159265358979323846
double gar()
{double r = 0.0;
int i;
for(i = 0; i < 10; i++){r += ((double)rand() / RAND_MAX) - 0.5;}
return r;}
int main(void)
{int x, k;　　　　　float d;　scanf("%f",&d);
double ek = PI * 2.0 * 0.7 / XMAX;　double a = 5.0 – 2*d;
double f[XMAX];　　double g[KMAX];
srand((unsigned)time(NULL));
for(k = 1; k < KMAX; k++){
g[k]=pow(( ek * rk ),-a/2);}
for(x = 0; x < XMAX; x++){　　f[x] = 0.0;
for(k = 1; k < KMAX; k++){　　double kx = ek * rk * x;
f[x] = f[x]+g[k]*(gar() * cos(kx) - gai() * sin(kx));}
printf("%g\n",f[x]);}
}