程序設(shè)計實驗報告（第二次）

程序設(shè)計實驗報告（第二次）

班級：機自七班

姓名：李鵬飛

學號：20110401726

求500以內(nèi)的所有親密對數(shù)。

⑴源程序代碼。

#include<stdio.h>

intfacsum(intm)

(

intsum=l,f=2;

for(f=2;f<=m/2;f++)

(

if(m%f==0)

sum+=f;

)

returnsum;

)

voidmain()

(

intm=3,n,k;

for(m=3;m<=500;m++)

(

n=facsum(m);

k=facsum(n);

if(m==k&&m<=n)

{printf(,,%d,%d\n,,,m,n);}

)

(2)運行結(jié)果。

■3°C:\Users\asus\Desktop\Debug\e2_l.exe'

6,6

28,28

220,284

496,496

Pressanykeytocontinue

(3)結(jié)果說明。

自定義一個函數(shù)，并調(diào)用該函數(shù)計算并顯示出500以內(nèi)的所有親

密對數(shù)。

二、利用復化梯形公式計算定積分。

<1>求第一個函數(shù)的定積分。

(1)源程序代碼。

mainl.c

#include<stdio.h>

#include<math.h>

#includeHsab.h"

doublef(doublex)

(

doubleresult;

rcsult=x*x*cxp(x);

returnresult;

main()

(

doublea,b,result;

intn;

printf(Hpleaseinputthevalueofa,b,n\nn);

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

result=sab(a,b,n);

printf("theresultisresult);

return0;

)

sab.c

#include<stdio.h>

doublef(doublex);

doublesab(doublea,doubleb,intn)

(

doubleh,result,x1,x2,x3=0,t;

intk;

h=(b-a)/n;

xl=f(a);

x2=f(b);

for(k=1;k<=n-1;k++)

(

t=a+k*h;

x3=x3+f(t);

result=h*(x1+x2)/2+h*x3;

)

returnresult;

)

Sab.h

doublef(doublex);

doublesab(doublea,doubleb,intn);

⑵運行結(jié)果。

■3"C:\USERS\ASUS\DESKTOP\2\3\Debug\3.exe'

pleaseinputtheualueofa,b,n

018

theresultis0.728890

Pressanykeytocontinue

<2>求第二個函數(shù)的定積分

(1)源程序代碼。

main2.c

#include<stdio.h>

#includeHsab.h"

doublef(doublex)

(

doubleresult;

result=l/(25+x*x);

returnresult;

)

main()

(

doublea,b,result;

intn;

printf("pleaseinputthevalueofa,b,n\n");

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

result=sab(a,b,n);

printf("theresultis%lf\n",result);

return0;

sab.c與sab.h同上。

(2)運行結(jié)果。

■3'C:\USERS\ASUS\DESKTOP\2\22\Debug\22.exe'

pleaseinputtheualueofa,b,n

kl110

theresultis0.078938

Pressanykeytocontinue

(3)結(jié)果說明。

此程序運用到了工程管理的思想。

三、計算Ackerman函數(shù)。

(1)源程序代碼。

#include<stdio.h>

intack(intn,intx,inty)

(

inta;

if(n==0)

a=x+1;

if(n==l&&y==0)

a=x;

if(n==2&&y==0)

a=0;

if(n==3&&y==0)

a=l;

if(n>=4&&y==0)

a=2;

if(n!=O&&y!=O)

a=ack(n-1,ack(n,x,y-1),x);

returna;

)

main()

{

intn,x,y,result;

printf("pleaseinputthevalueofn,x,y:\nM);

scanf(,'%d%d%d',,&n,&x,&y);

if(n>=0&&x>=0&&y>=0)

(

result二ack(n,x,y);

printf("ack(%d,%d,%d)=%d\n",n,x,y,result);

)

else

{printf(H輸入錯誤，請重新輸入\n")；

scanf(H%d%d%d",&n,&x,&y);

result=ack(n,x,y);

printf(,'ack(%d,%d,%d)=%d\n,,,n,x,y,result);

(2)運行結(jié)果。

"C:\Users\asus\Desktop\Debug\e2_3.exe'

hleaseinputtheualueofn,x『y：