实验五参考答案

实验五：参考答案

1、编程实现，求Sn=a+aa+aaa+...+aa...a 之值，其中a 是一个数字，n 表示a 的位数。如：2+22+222+2222+22222（此时n=5），n,a 由键盘输入。（C 语言源程序文件名为51.c ，程序运行结果文件名为51.txt ，本程序若a 是2，n 是5，程序结果为24690）

【参考程序】

#include

void main()

{int a,n,i=1,sn=0,tn=0;

FILE *fp;

clrscr();

fp=fopen("e:\\5\\51.txt","w");

printf("a,n=:");

scanf("%d,%d",&a,&n);

while(i<=n)

{tn=tn+a;

sn=sn+tn;

a=a*10;

i++;

}

printf("a+aa+aaa+...=%d\n",sn);

fprintf(fp,"a+aa+aaa+...=%d\n",sn);

fprintf(fp, "12345678\n")/*输出学号*/

fclose(fp);

}

2、用矩形法求 。（C 语言源程序文件名为52.c ，程序运行结果文件名为52.txt ）

【参考答案】

#include

void main()

{float x,s,t,h;

FILE *fp;

s=0;h=0.001;

clrscr();

fp=fopen("e:\\5\\52.txt","w");

for(x=0;x<3;x+=h)

{t=(x*x+3)*h; s+=t;}

printf("s=%f\n",s);

fprintf(fp,"s=%f\n",s);

fprintf(fp, "12345678\n")/*输出学号*/

fclose(fp); ?

+30

)3*(dx x x

}

}

3、用牛顿迭代法求方程2x3-4x2+3x-6=0的根。（C语言源程序文件名为53.c，程序运行结果文件名为53.txt，程序运行结果2.0）

【参考程序】

#include/* newton4.c */

#include

main()

{float x1,x,f,f1;FILE *fp;

scanf("%f",&x);

fp=fopen("e:\\5\\53.txt","w");

do{

x1=x;

f=2.0*x*x*x-4.0*x*x+3*x-6;

f1=6.0*x*x-8.0*x+3;

x=x1-f/f1;}

while(fabs(x-x1)>=1e-5);

printf("2x^3-4x^2+3x-6=0,Its root is%f\n",x);

fprintf(fp,"2x^3-4x^2+3x-6=0,Its root is%f\n",x);

fprintf(fp, "12345678\n")/*输出学号*/

fclose(fp);

}

4、任意读入x的值（0

【参考程序】

#include

#include

main()

{float x, s=0, t;FILE *fp;

do

{printf("Input 0

scanf("%f ", &x);

}while(x<=0|x>=1);

fp=fopen("e:\\5\\54.txt","w");

t=x;

while(fabs(t)>=1e-6)

{s=s+t;

t=t*(-x); /*递推：利用前项求后项*/

}

printf("%f\n", s);

fprintf(fp,"%f\n", s);

fprintf(fp, "12345678\n")/*输出学号*/

fclose(fp);

}

5、编程输出形如下图的平行四边形。

【要求：行数、每行星号数（即列数）由键盘输入。】

【提示：使用双重循环完成此题。外循环控制行数、内循环控制每行的空格数和星号数，每行输出的空格数比前一行多一个。】

******

******

******

******

（C 语言源程序文件名为55.c ，程序运行结果文件名为55.txt ）

【参考程序】

#include

main()

{int i, j, k;

int H, L;

FILE *fp;

fp=fopen("e:\\5\\55.txt","w");

scanf("%d%d",&H,&L);

for(i=1; i<=H; i++) /*控制输出H 行*/

{for(k=1; k<=i; k++){printf(" "); fprintf(fp," ");} /*每行输出若干个空格，每行比前一行多输出一个空格*/

for(j=1; j<=L; j++){printf("*");fprintf(fp,"*");} /*每行输出L 个空格*/ printf("\n"); fprintf(fp,"\n");

}

fprintf(fp, "12345678\n")/*输出学号*/

fclose(fp);

}

四、附加题：

6、计算输出下列式子的和（x 和n 由键盘读入，其中0=5）：

s=x+x 2/(1×2)+x 3/(2×3)+x 5/(3×5)+……+x fn /(f n-1×f n ) 即n n n

f f f x x x x x s ????-+++++=15

3

2 (533221)

其中：f i （i=0,1,2,3,……,n ）是Fibonacci 数列中的第i 项。

Fibonacci 数列是：1,1,2,3,5,8,13,21,34,……。

Fibonacci 数列的生成规则是： 当i>1时，f i =f i-1+f i-2

（C 语言源程序文件名为56.c ，程序运行结果文件名为56.txt. 当n=10，x=0.5时，运行结果： s=0.648015）

【参考程序】

#include

#include

main()

{float x, s, t, f1=1, f2=1, f3;

int n, i;

FILE *fp;

fp=fopen("e:\\5\\56.txt","w");

do

{printf("Input 0=5\n");

scanf("%f,%d",&x,&n);

}while(x<=0||x>=1||n<5); /*此do...while循环确保输入正确*/ s=x;

for(i=2; i<=n; i++)

{f3=f1+f2;

t=pow(x, f3)/(f2*f3);

s=s+t;

f1=f2; f2=f3;

}

printf("s=%f\n", s);

fprintf(fp,"s=%f\n", s);

fprintf(fp, "12345678\n");/*输出学号*/

fclose(fp);

}

程序也可以写成

#include

#include

main()

{float x, s, t, f1=1, f2=1, f3,a1,a2,a;

int n, i;

FILE *fp;

fp=fopen("e:\\5\\56.txt","w");

do

{printf("Input 0=5\n");

scanf("%f,%d",&x,&n);

}while(x<=0||x>=1||n<5); /*此do...while循环确保输入正确*/ s=x; a1=x;a2=x;

for(i=2; i<=n; i++)

{f3=f1+f2; a=a1*a2;

t=a/(f2*f3);

s=s+t;

f1=f2; f2=f3;

a1=a2;a2=a;

}

printf("s=%f\n", s);

fprintf(fp,"s=%f\n", s);

fprintf(fp, "12345678\n");/*输出学号*/

fclose(fp);

}

7、古代印度国王要褒奖他的聪明能干的宰相达依尔（国际象棋的发明者），问他要什么，达依尔回答：“陛下只要在国际象棋棋盘的第一个格子中放一颗麦子，第二个格子中放两颗麦子，以后每个格子中的麦子数都是前格的两倍，如此放满64格，我就心满意足了”。国王想这不难办到，但一袋麦子很快用完了，全印度的麦子都不够，请编程序计算所需麦子的体积(设1立方米的麦子约为1.4×108颗)。（C语言源程序文件名为57.c，程序运行结果文件名为57.txt）

【参考程序】

#include

void main()

{double s=1;

int i;

FILE *fp;

clrscr();

fp=fopen("e:\\5\\57.txt","w");

for(i=2;i<=64;i++)

s+=2*s;

printf("%lf\n",s/1.4e8);

fprintf(fp,"%lf\n",s/1.4e8);

fprintf(fp, "12345678\n")/*输出学号*/

fclose(fp);

}