質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.37%
C

C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

Q&A

3回答

3121閲覧

モンテカルロ法を使ったc言語での円の面積の求め方

yukamu704

総合スコア0

C

C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

0グッド

0クリップ

投稿2021/06/18 09:40

前提・実現したいこと

ここに質問の内容を詳しく書いてください。
(例)PHP(CakePHP)で●●なシステムを作っています。
■■な機能を実装中に以下のエラーメッセージが発生しました。

発生している問題・エラーメッセージ

モンテカルロ法でc言語を使い、半径2の円の面積を求める問題で、半径1の円の面積は求めることができたのですが、半径2の面積がどうしても出てきません。どう直したらいいのか、ご教授お願いします

該当のソースコード

c言語

1#include <stdio.h> 2#include <stdlib.h> 3#include <math.h> 4#include <time.h> 5 6int main(void){ 7 8 double x,y; 9 int n = 0; 10 double p; 11 12 srand(time(NULL)); 13 14 for(int i = 1; i <= 1000000; i++){ 15 16 x = rand() / (RAND_MAX + 1.0) * 2; 17 y = rand() / (RAND_MAX + 1.0) * 2; 18 19 if((pow(x - 2.0, 2.0) + pow(y - 2.0, 2.0)) <= 2){ 20 21 n++; 22 23 } 24 25 } 26 27 p = 4.0 * n /1000000; 28 29 printf("求める面積は%f\n",p); 30 31 return 0; 32 33}

試したこと

pow関数の中の数字を変えたり、x, yの*の後を変えたりしてみました

補足情報(FW/ツールのバージョンなど)

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

maisumakun

2021/06/18 09:41

現状では、どのような問題があるのでしょうか?
guest

回答3

0

半径1のができているなら,その結果を4倍するだけで良いのでは.

正直に径の分だけ広げた世界でモンテカルロしなくてもいいよね.
乱数が円の内側に落ちる確率は一緒なので.


投稿2021/06/18 09:58

編集2021/06/18 10:00
fana

総合スコア11954

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

0

修正しました。

C

1#include <stdio.h> 2#include <stdlib.h> 3#include <math.h> 4#include <time.h> 5 6int main(void){ 7 8 double x,y; 9 int n = 0; 10 double p; 11 12 srand(time(NULL)); 13 14 for(int i = 1; i <= 1000000; i++){ 15 16 x = rand() / (RAND_MAX + 1.0) * 2; 17 y = rand() / (RAND_MAX + 1.0) * 2; 18 19 if((pow(x , 2.0) + pow(y , 2.0)) <= 4.0){ //修正 20 21 n++; 22 23 } 24 25 } 26 27 p = 4.0 * n /1000000 * 4.0 ; //修正 28 29 printf("求める面積は%f\n",p); 30 31 return 0; 32 33} 34

実行結果
求める面積は12.570192

投稿2021/06/18 10:14

tatsu99

総合スコア5493

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

0

x = rand() / (RAND_MAX + 1.0) * 2;

x = rand() / ((RAND_MAX + 1.0) * 2);って事になると思いますが
このへんの検証は行いましたか?

投稿2021/06/18 09:45

asm

総合スコア15149

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.37%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問