C言語を実行すると-infが出てきて困っています。
コンパイラはgccを使っています。
よろしくお願いします。
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define data 100//計算回数
double seiki(void)//正規乱数
{
double n;
double i;
double y;
double x1,x2;
double sigma;
double mean;
double Pi = 2*asin(1);
sigma = 1;
mean = 0;
{
x1 = (double)rand()/(RAND_MAX);
x2 = (double)rand()/(RAND_MAX);
y = sigma*sqrt(-2*log(x1))*sin(2*Pi*x2) + mean;
//printf("%f\n",y);
}
return (y);
}
int main(void)
{
double* price;
int i, j, k;
double a, b, c;
//メモリ確保
price=(double*)malloc(sizeof(double)*data+10);
for(i = 0;i<data;i++)
{
price[i+1] = price[i]+seiki();
printf("%lf\n",price[i+1]);
}
return 0;
}
No.3
- 回答日時:
#1 に加えて mean や sigma (って何?) も初期化されていないですね.
ついでにいうと, seiki って関数は間違ってませんか? 「0 以上 1 未満」の乱数を使うのがセオリーだったような.
No.1
- 回答日時:
>price[i+1] = price[i]+seiki();
iが0のとき、値が定まっていないprice[0]を参照しています。
malloc()で確保した領域の内容が特定の値(例えば0x00)で
埋まっているとは限りません。
>x1 = (double)rand()/(RAND_MAX);
>x2 = (double)rand()/(RAND_MAX);
rand()を使うのであれば、その前に一度だけ
乱数の種を初期化しておく必要があります。
その他、気になる点:
>price=(double*)malloc(sizeof(double)*data+10);
最後の+10の意味がわかりません。
定義しただけで使っていない変数が多々見受けられます。
余分なコードはそぎ落としましょう。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# c言語でユーザ関数を利用して複素数のべき乗と絶対値の数列を計算するプログラムが作りたいです。 3 2023/01/29 22:13
- C言語・C++・C# C言語のマクローリン展開ローラン展開のコードについて 3 2022/12/15 14:45
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# Cのdoubleの浮動小数点表示について 3 2023/04/17 13:14
- C言語・C++・C# LU分解法のピボッティングについて(C言語/gcc-9) 3 2022/07/11 23:10
- C言語・C++・C# (C言語)めちゃくちゃな値になってしまいます。 5 2022/08/13 11:55
- C言語・C++・C# C言語のエラーについて 2 2022/07/11 13:56
- C言語・C++・C# C 言語の Gauss Jordan 法について 2 2022/12/28 11:16
このQ&Aを見た人はこんなQ&Aも見ています
-
外出時に「待たせる妻」vs イライラする「待つ夫」は日本だけ?見習いたい海外事情
夫の家事参加に積極的なイメージのある海外でも、同様の事例はあるのか。結婚カウンセラーの佐竹悦子さんに伺ってみた。
-
C言語 配列の長さの上限
C言語・C++・C#
-
関数から配列を返すには?
C言語・C++・C#
-
C言語初心者の質問失礼します。
C言語・C++・C#
-
-
4
<math.h>があるのにsqrtが・・・
C言語・C++・C#
-
5
e^(x^2)の積分に関して
数学
-
6
C言語の2次元配列 容量が大きすぎる場合の対処方法
C言語・C++・C#
-
7
0除算して、落ちるプログラムと落ちないプログラム
C言語・C++・C#
-
8
分数の足し算をさせるプログラムが分かりません。どなたか分かりませんか?
C言語・C++・C#
-
9
C言語で-23乗を取り扱うには
C言語・C++・C#
-
10
c言語でのfscanfについて
C言語・C++・C#
-
11
C言語 exitの使い方
C言語・C++・C#
-
12
文字列から空白を取り除きたいのですが
C言語・C++・C#
-
13
構造体のメンバをfor文で回したい
C言語・C++・C#
-
14
シグマ公式・・・C言語
C言語・C++・C#
-
15
word 数式が一部表示されない
Word(ワード)
-
16
∫1/(x^2+1)^2 の不定積分がわかりません
数学
-
17
2分法で方程式の複数の解を自動的に求めるには?
C言語・C++・C#
-
18
相対誤差が小さいと判断する基準がわからないのでどのような値になったら小さいと判断してよいのか教えてほ
大学・短大
-
19
C言語の問題について
C言語・C++・C#
関連するカテゴリからQ&Aを探す
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
float型とdouble型の変数の違い...
-
C言語を実行すると-infが出てき...
-
doubleの変数にintとintの割り...
-
プログラムでの数字につく”f”の...
-
C 開放してるのにエラー(doubl...
-
doubleは常に%lfとするべきなのか
-
ニュートン法
-
相互相関関数
-
Cで3乗根を求める方法
-
floating point not loadedとは?
-
C言語で-23乗を取り扱うには
-
関数におけるif文とreturn文に...
-
C言語のプログラムで#include<m...
-
カウントアップタイマ
-
テイラー展開(C言語)
-
「割り算」 と 「分数の掛け算」
-
線形補間
-
物体が往復する動きを作りたい
-
EXE1→DLL→EXE2数値を受け渡す方法
-
-1.#IND00と出てしまうのですが...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
float型とdouble型の変数の違い...
-
doubleの変数にintとintの割り...
-
C言語を実行すると-infが出てき...
-
C 開放してるのにエラー(doubl...
-
至急です! マクロ定義で #defi...
-
c言語で、繰り返し文の中で、0....
-
関数におけるif文とreturn文に...
-
C言語 関数プロトタイプ宣言の...
-
C言語初心者 構造体 課題について
-
C言語の型による処理速度の違い
-
Cで3乗根を求める方法
-
C言語で-23乗を取り扱うには
-
2分法で方程式の複数の解を自...
-
doubleは常に%lfとするべきなのか
-
c言語のコンパイルエラー canno...
-
C言語で直角三角形の斜辺を求め...
-
C言語のプログラムで#include<m...
-
int とdoubleの比較
-
C++で外積
おすすめ情報