请编写函数fun()，其功能是计算并输出当x＜0.97时下列多项式的值，直到|Sn-Sn-1|＜0.000001为止。
Sn=1+0.5x+0.5(0.5-1)x2/2!+0.5(0.5-1)(0.5-2)x3/3!+…0.5(0.5-1)(0.5-2)…(0.5-n+1)xn/n!
例如：主函数从键盘给x输入0.21后，则输出为s=1.100000。
注意：部分源程序以存在文件test38_2.cpp中。
请勿改动主函数main和其他函数中的任何内容，仅在函数fun的花括号中填入所编写的若干语句。
文件test38_2.cpp的内容如下：
#include＜stdio.h＞
#include＜iostream.h＞
#include＜math.h＞
double fun(double x)
void main( )
double x, s;
cout＜＜"Input x:"＜＜endl;
cin＞＞x;
s=fun(x);
cout＜＜"s="＜＜s＜＜endl;