- 締切済み
C言語 二分法
初投稿です。 お恥ずかしながらパソコンが苦手で、Cゲッが難しくてできません。 今回二分法です。 途中まではやったのですができません。 演習 0～1の乱数を12個発生させ,これらの平均をxi,yi とする。(s,tは乱数) xi=1/12(x1i+x2i+....x12i) yi=1/12(y1i+y2i+....y12i) このようなを1000個作り,(xi,yi)で散布図 を作りなさい。またx,yのそれぞれの平均を求 めよ。 この演習で #include<stdio.h> #include<stdlib.h> #include<math.h> #define eps 1.0e-5 float f(double x); void nibuin(void); int main () { int count; double a,b,m; count=0; printf("範囲の左の値を入力してください。\n"); scanf("%lf",&a); printf("範囲の右の値を入力してください。\n"); scanf("%lf",&a); if(count==1000){ printf("収束しませんでした。\n"); exit(1); } } while(!(fabs(a-b)<eps)); printf("解の値は%f\n収束するのに%d回かかりました。"m,count); } float f(double x) { reurn x*sin(x)+log(x); } まではできたのですが、 scanf("%lf",&a);とif(count==1000){の間に入る命令が打てません。 よろしくお願いします。
- masa55511
- お礼率0% (0/3)
- 数学・算数
- 回答数3
- ありがとう数0
- みんなの回答 （3）
- 専門家の回答
みんなの回答
- FT56F001
- ベストアンサー率59% (355/599)
二分法というタイトルと書きかけのプログラムは，雰囲気的に関係ありそうなので， 乱数の12個の平均を出す演習問題は入れ間違いで， 「x*sin(x)+log(x)=0を二分法で解くプログラムを作れ」 といった二分法の問題があるのかしら? 二分法で方程式f(x)=0を解くのなら， 1) f(a)*f(b)<0となる下限aと上限bを読み込む。 2) m=(a+b)/2に対してf(m)を計算する。 3) f(m)*f(a)>0なら下限aをmに置き換え，そうでなければ上限bをmに置き換える。 4) 2)～3)を，abs(a-b)がepsより大きく，かつ，カウントが1000以下である限り，繰り返す。 をコーディングすればよいわけです。
関連するQ&A
- 二分法のプログラムについて
下の用なプログラムを作ったのですがどうしても正しい答えを導くことができません。自分でもいろいろ調べてみましたがわかりません。誰かご教授宜しくお願いします。 #include<stdio.h> #include<stdlib.h> #define MAX 10 int n , count; double c[MAX+1]; double a,b,e; void nyuuryoku(void) { int i; printf("nの入力>"); scanf("%d",&n); if(n>MAX){printf("最大次数を超えている");exit(1);} else if(n<0){printf("nが負");exit(2);} else{for(i=0;i<=n;i++){printf("係数の値>");scanf("%lf",&c[i]);} }} double f(double x) {double y; int i; y = c[0]; for(i=1;i<=n;i++){ y=y*x+c[i];} return y; } void hani(void){ printf("aの値>");scanf("%lf",&a); printf("bの値>");scanf("%lf",&b); printf("eの値>");scanf("%lf",&e); if(e<=0){printf("eが0または負"); exit(3);} if(f(a)==0){printf("%f",f(a)); exit(4);} if(f(b)==0){printf("%f",f(b)); exit(5);} if(f(a)*f(b)>0){printf("初期値異常"); exit(6);}} double nibun(void) {double c; if(b>a){ while(b-a>e){ count++; c=(a+b)/2; if(f(c)==0){ return c;} if(f(a)*f(c)<0){b=c;} if(f(b)*f(c)<0){a=c;} } return a;} if(a>b){ while(a-b>e){ count++; c=(a+b)/2; if(f(c)==0){ return c;} if(f(b)*f(c)<0){a=c;} if(f(a)*f(c)<0){b=c;} } return a;} } void syutsuryoku(double x){ printf("x=%lf\n",x); printf("f(x)=%lf\n",f(x)); printf("繰り返し回数=%d\n",count); } int main(void){ double ans; count=0; nyuuryoku(); hani(); ans = nibun(); syutsuryoku(ans); }
- 締切済み
- その他（プログラミング・開発）
- C言語プログラム(二分法)について質問です。
C言語プログラム(二分法)について質問です。 以下作成したプログラムでは、aを入力すると, x=b - (log(a)) - (a+x)/2の解が求まります。（今回式は適当ですが。） 二分法で解を求めるというプログラムは作成できたのですが、 このプログラムで、aを１～１０まで変化させたときのxの値というようなループをプログラムでしたうえで、 a*xの値がもっとも大い点を求めるためにはどのようなプログラムを組めば（これを改良すれば）いいのでしょうか？ 具体的に行いたいことは、この二分法のプログラムをaとxの関数とし、a*xの最大点を求めたいのです。 質問が分かりにくいかもしれませんが、お願いします。 以下作成したプログラムです。 #include <stdio.h> #include <stdlib.h> #include <math.h> #define eps 1.0e-6 double b=0.3; double a=0; double f(double x); void nibun(void); int main() { nibun(); return 0; } void nibun(void) { int count; double x0,x1,m; printf("aの値\n"); scanf("%lf",&a); count=0; x0=-2*b; x1=b; do { count++; m=(x0+x1)/2.0; if(f(m)*f(x0)<0) x1=m; else x0=m; if(count==700) { printf("Error\n"); exit(1); } } while (!(fabs(x0-x1)<eps)); printf("解の値は %f\n",m); } double f(double x) { return(b - (log(a)) - (a+x)/2); }
- ベストアンサー
- C・C++・C#
- c言語
#include<stdio.h> #define N 5 int main(){ double x[] = {1, 3, 4, 5.5, 7}; double y[] = {3, 5, 6, 8, 11}; double a = 0.028; double eta = 0.001; double dEda; int ite, i ; for (ite = 1; ite <= 1000; ++ite){ dEda = 0; for(i = 0; i < N; i++){ dEda = dEda + (a * x[i] - y[i]) * x[i]; } a = if (ite == 1 || ite == 10 || ite == 100) printf("%4d: a = %.8lf\n", ite, a); } printf("%4d: a = %.8lf\n", 1000, a); return 0; } このプログラムの a = の部分にa-(dE/da)*((d^2)E/d(a^2))^-1を代入したいのですがうまくいきません。((d^2)E/d(a^2))の値を格納する変数を用意してちゃんと計算できるようにするにはどうすればよいでしょうか？ 補足 dE/da = Σ(ax-y)*x (d^2)E/d(a^2) = Σx^2 という計算になります。
- 締切済み
- C・C++・C#
- C言語のプログラミングがうまくいきません！
先日のプログラムを作り直しました。アドバイスのおかげでエラーは減りましたが、まだ完ぺきではないようです。恐らく、elseのところがいけないと思うのですが、どうしたらいいか分かりません。 ご指摘お願いします。 #include <stdio.h> #include <math.h> void main(void) { double x ,y ,z ,error ,menseki ; double ans; printf("x ?"); scanf("%lf", &x); printf("y ?"); scanf("%lf", &y); printf("z ?"); scanf("%lf", &z); if(x<y+z || y<x+z || z<x+y) { ans=(x+y+z)/2; menseki=sqrt(ans*(ans-x)*(ans-y)*(ans-z)); } if(x==y && y==z && z==x) { printf("正三角形です"); printf("面積は%lfです" , menseki); } else if(x==y || y==z || z==x) { printf("二等辺三角形です"); printf("面積は%lfです" , menseki); } else { printf("三角形です"); printf("面積は%lfです" , menseki); } else { printf("error"); } }
- ベストアンサー
- C・C++・C#
- Ｃ言語課題（数値計算）お願いします。
課題内容 ニュートン法を用いて√２の値を求められるプログラムを作る。 １０ケタまで表示する。 書いてみたソースコード #include <stdio.h> #include <stdlib.h> #include <math.h> #define max 1000 //最大繰り返し回数 #define eps 1.0e-10 //収束条件 double f(double x); double df(double x); void newton(void); int main() { newton(); return 0; } void newton(void) { int count; double a,newa; count=0; printf("初期値を入力してください。\n"); scanf("%lf",&a); for(;;) { count++; newa=a-f(a)/df(a); if(fabs(newa-a)<eps) break; a=newa; if(count==max) { printf("収束しませんでした。\n"); exit(1); } } printf("解の値は %f\n収束するのに %d 回かかりました。\n", newa,count); } double f(double x) { return x*x-2; } double df(double x) { return 2x; } コンパイルできません。 間違いとどうすればよいかを教えていただけるとありがたいです！！！
- 締切済み
- C・C++・C#
- C言語による数値計算を教えて下さい
「2次関数y = ax2 + bx + c の係数を最小2乗法で求めるための正規方程式を導出せよ．その上でキーボードからデータの組(xi, yi)を 5点入力してa,b,cを求めるプログラムを作成せよ(x1～x5,y1～y5)」 という問題なのですが、a,b,cを求める方程式を行列式の形で以下のように出せました。 |Σxi^4 , Σxi^3 , Σxi^2| |a| |Σxi^2*yi| |Σxi^3 , Σxi^2 , Σxi |*|b| = |Σxi*yi | |Σxi^2 , Σxi , 5 | |c| |Σyi | ですがこれをガウス消去法で解こうと思っているのですがどうにもプログラムができません。 どなたかC言語で数値計算できる方、助けてください。 以下のはサイエンス社の「C言語による数値計算入門」をもとに ガウス消去法で解くプログラムを自分なりに作ったのですが 当然エラーだらけです。 プログラムを作れる方、どうか教えてください。 #include <stdio.h> #include <math.h> /* ガウス消去法 */ double *simple_gauss( double a[3][3], double b[3] ); main(void) { double a[3][3]; /* マトリックスの宣言*/ double x[3]; /* 解ベクトルの宣言*/ double y[3]; /* 右辺ベクトルの宣言*/ double xx, yy; /* キーボード入力のための一時変数*/ int i, j; /* 繰り返しのための一時変数*/ for(i = 0; i < 2; i++) { /* 変数の初期化*/ x[i] = y[i] = 0; for(j = 0; j < 3; j++) a[i][j] = 0; } for(i=0; i<5; i++) { printf("\n(x, y) = "); scanf("%lf,%lf", &xx, &yy); a[1][1] += xx * xx * xx * xx; a[1][3] += xx * xx; a[2][2] = a[3][1]=a[1][3]; a[1][2] += xx * xx * xx; a[2][1]=a[1][2]; a[2][3] += xx; a[3][2]=a[2][3]; a[3][3] ++; y[2] += xx * yy; y[1] += xx * xx * yy; y[3] += yy; } y=*simple_gauss(a,y); return 0; } double *simple_gauss( double a[3][3], double y[3] ) { int i, j, k; double alpha, tmp; /* 前進消去 */ for( k = 0; k <= 1; k++) { for( i = k+1; i <= ; i++) { alpha = - a[i][k]/a[k][k]; for( j = k+1; j <= 2; j++) { a[i][j] = a[i][j] + alpha * a[k][j]; } y[i] = y[i] + alpha * y[k]; } } /* 後退代入 */ y[2] = y[2]/a[2][2]; for( k = 1; k >= 0; k--) { tmp = y[k]; for( j = k+1; j <= 2; j++) { tmp = tmp - a[k][j] * y[j]; } y[k] = tmp/a[k][k]; } return y; }
- ベストアンサー
- 数学・算数
- C言語の関数に関する質問ですが
C言語の初心者です。よろしくお願いいたします。 授業でこのような演習が出ました。 演習：実数x を入力したときの最大値を求めるプログラムを作れ． 実数x を入力すると，x; -x; x2; xの絶対値の平方根 の中で一番大きい値を答える プログラムを作れ（ファイル名はmax.c とする）． 表示は以下のようにする． Input x: -0.5 【Enter】 Answer is 0.707107. #include<stdio.h> #include<math.h> double max(double a, double b){ if( a > b) return a; else return b; } int main(void) { double x,y; printf(\"Input x: \"); scanf(\"%lf\",&x); y = max (x,-x); y = max (y,x*x); y = max (y,sqrt(fabs(x))); printf(\"Answer is %f.\\n\",y); } このように書けばうまく実行できますが、関数の中に関数を使えないでしょうか。うまく言えないですが、たとえば、以下のように書いてみましたが、うまく実行できません。どう直したらいいでしょうか、お忙しい中教えていただけたらうれしいです。 #include <stdio.h> #include <math.h> int max(double a,double b) { if (a<b) return b; else return a;} int main(void) { double x,result; printf(\"Input x:\"); scanf(\"%lf\",&x); result=max(max(x,-x),max(pow(x,２),sqrt(fabs(x)))); printf(\"%.2f\",result); return 0; } よろしくお願いいたします！！
- ベストアンサー
- C・C++・C#
- C言語の演習２
何度も似たような質問を出していますがよろしくお願いします。。。 問題は ２つの整数を入力してください。 整数Ａ：５４ 整数Ｂ：８４ Ａの値はＢの６４．２８５７１４％です。 というものです。あっているか確認お願いします。 #include<stdio.h> int main(void) { int vx, vy; puts("２つの整数を入力してください。"); printf("整数Ａ"); scanf("%lf",&vx"); printf("整数Ｂ"); scanf("%lf",&vy"); printf("Ａの値はｂの%f%です。/n",(double)(vx/vy)); return 0; } であっていますか？？お願いします。
- 締切済み
- C・C++・C#
- C言語のプログラムを見てください
ある100行の値がx列、y列の2列あるファイルを読み込んでそれを配列に入れ、yの最小値及びそれと同じ行にあるxの値を表示するプログラムを書きたいのですがy列の最小値を表示するプログラムを書き終えた所でコンパイルして実行してみると正しく値が表示されませんでした。それどころか実行するたびに値が変わってしまいます。どこがおかしいのかわからないため、ご指摘のほどよろしくお願いします。また、できれば同じ行にあるx列の値も表示させるプログラムを教えてください。 よろしくお願いします。 #include <stdio.h> #include <stdlib.h> #define N 100 int main(void) { int x[N],i; double y[N],min; FILE *fp; fp=fopen("book.dat","r"); if(fp==NULL){ puts("can't open file!"); exit(-1); } for(i=0;i<N;i++){ fscanf(fp,"%d %lf", &x[N],&y[N]); printf("x=%d\n y=%lf\n",x[N],y[N]); } min=y[0]; for(i=1;i<N;i++){ if(y[i]<min) min=y[i]; } fclose(fp); printf("最小値:%lf\n",min); return 0; }
- 締切済み
- C・C++・C#
- プログラムの作成
Cで、キ－ボ－ドから整数を1個(mとする)、実数を1個(xとする)それぞれ入力し、 S＝Σ(n＝0からmまで)1/2n+1×(x－1/x+1)^2n+1 の値を求めるプログラムを作ってみたのですが、あってるかどうかわかりますか？ #include<stdio.h> #include<math.h> main() { double y＝0; int i,m,x; printf(\" Input m: \") scanf(\"%lf\", &m); printf(\" Input x: \") scanf(\"%lf\", &x); for(i＝0; i <＝ m; i++) { y +＝ 1/(2*i+1)*pow((x－1)/(x+1),2*i+1) } printf(\" S ＝ %d＼n\", y) }
- ベストアンサー
- その他（プログラミング・開発）
