请编写函数fun，其功能是：计算并输出当x＜0.97时，如下多项式的值，直到|Sn－Sn－1|＜0.000001为止。

#include  <math.h>
#include  <stdio.h>
double fun(double x)
{
    double Sn=1.0,s1=0,t=1.0,b=1.0,r=1.0;
	int i;
	if(x<0.97)
	{
		for(i=1;;i++)
		{
			s1=Sn;
			b*=(0.5-i+1);
			t*=i;
			r*=x;
			Sn+=b/t*r;
			if(fabs(Sn-s1)<1e-6)
			  break;
		}
	} 
	return Sn;
}
void main()
{ int i;
  double x,s;
  FILE *out;
  printf("Input x: ");
  scanf("%lf",&x);
  s=fun(x);
  printf("s=%f\n ",s);