//アルキメデスの円周近似

#include<stdio.h>
#include<math.h>

int main(void)
{
int i;
long double a[100],b[100];
a[0]=4.0*sqrt(3.0);
b[0]=6.0;

for(i=0;i<99;i++)
{
a[i+1]=2*a[i]*b[i]/(a[i]+b[i]);
b[i+1]=sqrt(b[i]*a[i+1]);

printf("6X2^%d角形を用いた近似値：%Lf~%Lf\n",
i,b[i+1]/2,a[i+1]/2);
}
return 0;

}
桁数が足りなくなるのですがどう修正すればいいですか？