円周率を計算したいのですが、
#include<stdio.h>
#include <math.h>

int main(void)
{
int k;
double a,b,c,d,e,f;
a=0;b=0;c=0;d=0;e=0;f=0;


for(k=0;k<=1000000;k++)
{
f=1/pow(16,k);
b=4/(8*k+1);
c=2/(8*k+4);
d=1/(8*k+5);
e=1/(8*k+6);
a=a+f*(b-c-d-e);
}
printf("%f",a);
}
実行結果が４．００００となってしまいます。どこがおかしいかどなたか
教えていただけませんか？