sqrt( )= 精度 :1.0e-6 のとき sqrt( )= 精度 :1.0e-8 のとき sqrt( )= 精度 :1

Size: px

Start display at page:

Download "sqrt( )= 精度 :1.0e-6 のとき sqrt( )= 精度 :1.0e-8 のとき sqrt( )= 精度 :1"

なぎささわまつ
5 years ago
Views:

1 /******* 演習問題 1********************************** /***************en3-4_1.c****************/ doublecount; doublea,b; count=7.0; a=1.0; b=a+count/2.0; printf("%lf n",b); ***** / /*** ************ %ccen3-4_1.c *********************** /******* 演習問題 2********************************** /****************en3-4_2.c*********** #definemax1000 #defineaccuracy1.0e-1 doublea; doublea1,a2; doubletyu_ten; intcount; /*************** 数値入力と初期化 ***************/ printf("inputdata>>"); scanf("%lf",&a); a1=0.0; a2=a; /***************2 分法の繰り返し ****************/ for(count=0;counttyu_ten=(a1+a2)/2.0; if(tyu_ten*tyu_ten==a) if(tyu_ten*tyu_ten-a>0.0) a2=tyu_ten; else a1=tyu_ten; if(((a2-a1)/a1-accuracy)) printf("sqrt(%lf)=%lf n",a,tyu_ten); ************** 結果 (ACCURACY の値を変えたときの結果 ) 精度 :1.0e-1 のとき sqrt( )= 精度 :1.0e-2 のとき sqrt( )= 精度 :1.0e-4 のとき sqrt( )=

2 sqrt( )= 精度 :1.0e-6 のとき sqrt( )= 精度 :1.0e-8 のとき sqrt( )= 精度 :1.0e-10 のとき sqrt( )= をそれぞれデータとして入力してみたその結果精度が 1.0e-6 以上になると結果が同じになったそこで 2 よりもおおきい 5 の平方根をを精度が 1.0e-6 以上にしたときに調べたところ 5 の場合は精度が 1.0e-8 以上になると結果が同じになったまた標準関数を用いて求められる平方根は下の値のようになった上の結果と下の結果を見比べると標準関数の精度は 1.0e-6 から 1.0e-8 の間にあるのではないかと考えられたそのため再び精度が 1.0e-7 の時の 5 の平方根を求めたその結果 sqrt( )= となり標準関数の平方根の値と一致したため標準関数の精度は 1.0e-7 位ではないかと考えられる標準関数で平方根を求めた結果 sqrt( )= ************** /******* 演習問題 3********************************** en3-4_3.c****************/ #include<math.h> doublea,b,c; doubled,r,s,bunbo; doublex,x1,x2; printf("a=");scanf("%lf",&a); printf("b=");scanf("%lf",&b); printf("c=");scanf("%lf",&c); d=pow(b,2.0)-4.0*a*c; bunbo=2.0*a; if(a==0) x=-c/b; printf("x=%lf n",x); else if(d>0) x1=(-b+sqrt(d))/bunbo; x2=(-b-sqrt(d))/bunbo; printf("x1=%lf n",x1); printf("x2=%lf n",x2); elseif(d==0) x=(-b)/bunbo; printf("x=%lf n",x); else r=(-b)/bunbo; s=sqrt(-d)/bunbo;

3 printf("x1=%lf+%lfi n2",r,s); printf("x2=%lf-%lfi n",r,s); a=0 b=4 c=6 x= /******* 演習問題 4********************************** /***ensyu3-4.c***************/ doublenumber=1.2345e6; printf("%10.3lf n",number); printf("%10.3le n",number); %./ensyu e+06 /******* 演習問題 5********************************** /***************Ensyu3-5.c************/ doublematrix1[3][3],matrix2[3][3]; doublematrixsum[3][3],matrixmult[3][3]; intx; inty; for(x=0;x<3;x=x+1) for(y=0;y<3;y=y+1) printf("matrix[%d][%d]=",x,y); scanf("%lf",&matrix1[x][y]); matrixsum[x][y]=0; matrixmult[x][y]=0; for(x=0;x<3;x=x+1) for(y=0;y<3;y=y+1) printf("matrix[%d][%d]=",x,y); scanf("%lf",&matrix2[x][y]); /************Calculations***********/ for(x=0;x<3;x=x+1) for(y=0;y<3;y=y+1) matrixsum[x][y]=matrix1[x][y]+matrix2[x][y]; matrixmult[x][y]=matrix1[0][y]*matrix2[x][0]+matrix1[1][y]*matrix2[x][1]+matrix1[2][y]*matrix2[x][2]; /***********ResultDisplay*********/

4 printf("matrix1 n"); for(y=0;y<3;y=y+1) printf("%lf t%lf t%lf n",matrix1[y][0],matrix1[y][1],matrix1[y][2]); printf(" n n n"); printf("matrix2 n"); for(y=0;y<3;y++) printf("%lf t%lf t%lf n",matrix2[y][0],matrix2[y][1],matrix2[y][2]); printf("***************************************** n"); printf("matrixsum n"); printf("%lf t%lf t%lf n",matrixsum[0][0],matrixsum[0][1],matrixsum[0][2]); printf("%lf t%lf t%lf n",matrixsum[1][0],matrixsum[1][1],matrixsum[1][2]); printf("%lf t%lf t%lf n",matrixsum[2][0],matrixsum[2][1],matrixsum[2][2]); printf("******************************************** n"); printf("matrixmultiplication n"); printf("%lf t%lf t%lf n",matrixmult[0][0],matrixmult[0][1],matrixmult[0][2]); printf("%lf t%lf t%lf n",matrixmult[1][0],matrixmult[1][1],matrixmult[1][2]); printf("%lf t%lf t%lf n",matrixmult[2][0],matrixmult[2][1],matrixmult[2][2]); matrix[0][0]=1.1 matrix[0][1]=1.2 matrix[0][2]=1.3 matrix[1][0]=1.4 matrix[1][1]=1.5 matrix[1][2]=1.6 matrix[2][0]=1.7 matrix[2][1]=1.8 matrix[2][2]=1.9 matrix[0][0]=1.1 matrix[0][1]=1.2 matrix[0][2]=1.3 matrix[1][0]=1.4 matrix[1][1]=1.5 matrix[1][2]=1.6 matrix[2][0]=1.7 matrix[2][1]=1.8 matrix[2][2]=1.9 matrix matrix ***************************************** Matrixsum ******************************************** MatrixMultiplication /******* 演習問題 6************************* 6.(a) プログラム /*rep3-6.c*/ #include

5 #include intn,i; doublesum=0,s,ss=0,sss; printf("inputn>>");scanf("%d",&n); for(i=1;i<=n;i++) sum=sum+pow(i,-2); s=sum/n; for(i=1;i<=n;i++) ss=ss+pow(pow(i,-2)-s,2); ss=ss/n; sss=pow(ss,0.5); printf("sum=%lf ns=%lf nss=%lf nsss=%lf n",sum,s,ss,sss); Inputn>>18 sum= s= ss= sss= (b)1.645 付近で収束する 6.(c) プログラム /*rep3-6.c*/ #include #include #definemax inti,k; doublesum=0,s,ss=0,sss; printf("inputk>>");scanf("%d",&k); for(i=1;i<=max;i++) sum=sum+pow(i,-k); printf("sum=%lf n",sum); /*** 結果 Inputk>>2sum= Inputk>>3sum= Inputk>>4sum= Inputk>>5sum= Inputk>>6sum= Inputk>>7sum= Inputk>>8sum= Inputk>>9sum= Inputk>>10sum= いずれも (n= ) での値である k が大きくなると 1 に近付いていく ********* /******* 演習問題 7************************* /*ensyuu-3-7.c*/

6 #defineguu1 #definetyoki2 #definepaa3 inta,b; printf("inputnumber:guu(1),tyoki(2),paa(3)>>"); scanf("%d",&a); printf("inputnumber:guu(1),tyoki(2),paa(3)>>"); scanf("%d",&b); if(a!=b) switch(a) caseguu: if(b==tyoki)printf("winnerisa! n"); elseif(b==paa)printf("winnerisb! n"); casetyoki: if(b==paa)printf("winnerisa! n"); elseif(b==guu)printf("winnerisb! n"); casepaa: if(b==guu)printf("winnerisa! n"); elseif(b==tyoki)printf("winnerisb! n"); else printf("dorowgame! n"); return(0); %ccensyuu-3-7.c DorowGamed %ccensyuu-3-7.c DorowGame! InputNumber:GUU(1),TYOKI(2),PAA(3)>>2 WInnerisa! InputNumber:GUU(1),TYOKI(2),PAA(3)>>3 Winnerisb!

Microsoft PowerPoint - program.ppt [互換モード]

Microsoft PowerPoint - program.ppt [互換モード] プログラミング演習バージョン 1 担当教員 : 綴木馴プログラムの決まりについて学ぶおすすめする参考書ザ C 戸川隼人サイエンス社本日の予定 1. 授業の説明. 2. コンパイラーのインストール. プログラムの決まりについて学ぶ,P31 /* The most in C */ /* hello.c */ printf("hello,world n"); プログラムの決まり ( コメント )