
計算がうまくいかないのですがどこを直せばよいでしょうか?
#include <stdio.h>
#include <math.h>
void input_cmplx(double *q)
{
printf("q_re = ");
scanf("%lf", q);
q++;
printf("q_im = ");
scanf("%lf", q);
}
void conj_cmplx(double *q, double *q_conj)
{
*q_conj = *q;
*(q_conj + 1) = - *(q + 1);
}
void product_cmplx(double *q1, double *q2,double *q_out)
{
*q_out = *q1 * *q2 - *(q1 + 1) * *(q2 + 1);
*(q_out + 1) = *q1 * *(q2 + 1) + *(q1 + 1) * *q2;
}
double abs_cmplx(double *q)
{
double b[2], c[2];
double *q1, *q2;
q1 = b;
q2 = c;
conj_cmplx(q, q1);
product_cmplx(q, q1, q2);
return pow(*b + *c, 0.5);
}
int main (void)
{
int N, i;
double *q1, *q2, *q3;
double c1[2], c2[2], c3[2];
double abs_cmplx(double *q);
q1 = c1;
q2 = c2;
q3 = c3;
input_cmplx(q1);
printf("項数N = ");
scanf("%d", &N);
for (i = 0; i < N; i++)
{
product_cmplx(q1, q2, q3);
printf("%+.3e %+.3ei", c3[0], c3[1]);
printf(" %+.3e\n",abs_cmplx(q1));
q1 = q3;
}
return 0;
}
A 回答 (3件)
- 最新から表示
- 回答順に表示
No.3
- 回答日時:
ちょっと見ただけでの疑問です
int main (void)
{
int N, i;
double *q1, *q2, *q3;
double c1[2], c2[2], c3[2];
double abs_cmplx(double *q); ← 書く必要ないのでは?
q1 = c1;
q2 = c2;
q3 = c3;
input_cmplx(q1); ← input_cmplx(c1); でよくない?
そもそも double *q1, *q2, *q3; もいらないのでは?
product_cmplx(q1, q2, q3); ← q2の実体C2はいつ値をセットしたの?
No.2
- 回答日時:
本当ならデバッガでいろいろな変数の値を追い掛けて「何がおかしいのか」を自分で調べるべきだ. この場合なら printf デバッグでもいけると思うけど.
なお C で「明示的に初期化しない自動変数の初期値は不定」だってことは理解してるよね?
No.1
- 回答日時:
q1 = q3;
の意味を理解できているのかな?
なお「計算がうまくいかない」ではなく
・どのようなデータを与えて
・どのような結果を期待したところ
・どのような結果になったのか
を具体的に書いてほしい.
complex で全てを解決すりゃいいんじゃね~かなぁ....
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# C言語のマクローリン展開ローラン展開のコードについて 3 2022/12/15 14:45
- C言語・C++・C# Cのdoubleの浮動小数点表示について 3 2023/04/17 13:14
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# C 言語の Gauss Jordan 法について 2 2022/12/28 11:16
- C言語・C++・C# LU分解法のピボッティングについて(C言語/gcc-9) 3 2022/07/11 23:10
- C言語・C++・C# C言語(構造体) 3 2022/07/05 20:08
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# LU分解法のピボット選択機能実装について(C言語・gcc-9) 1 2022/07/22 15:20
- C言語・C++・C# C言語のエラーについて 2 2022/07/11 13:56
関連するカテゴリからQ&Aを探す
今、見られている記事はコレ!
-
弁護士が語る「合法と違法を分けるオンラインカジノのシンプルな線引き」
「お金を賭けたら違法です」ーーこう答えたのは富士見坂法律事務所の井上義之弁護士。オンラインカジノが違法となるかどうかの基準は、このように非常にシンプルである。しかし2025年にはいって、違法賭博事件が相次...
-
釣りと密漁の違いは?知らなかったでは済まされない?事前にできることは?
知らなかったでは済まされないのが法律の世界であるが、全てを知ってから何かをするには少々手間がかかるし、最悪始めることすらできずに終わってしまうこともあり得る。教えてgooでも「釣りと密漁の境目はどこです...
-
カスハラとクレームの違いは?カスハラの法的責任は?企業がとるべき対応は?
東京都が、客からの迷惑行為などを称した「カスタマーハラスメント」、いわゆる「カスハラ」の防止を目的とした条例を、全国で初めて成立させた。条例に罰則はなく、2025年4月1日から施行される。 この動きは自治体...
-
なぜ批判コメントをするの?その心理と向き合い方をカウンセラーにきいた!
今や生活に必要不可欠となったインターネット。手軽に情報を得られるだけでなく、ネットを介したコミュニケーションも一般的となった。それと同時に顕在化しているのが、他者に対する辛らつな意見だ。ネットニュース...
-
大麻の使用罪がなかった理由や法改正での変更点、他国との違いを弁護士が解説
ドイツで2024年4月に大麻が合法化され、その2ヶ月後にサッカーEURO2024が行われた。その際、ドイツ警察は大会運営における治安維持の一つの方針として「アルコールを飲んでいるグループと、大麻を吸っているグループ...
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
C言語を実行すると-infが出てき...
-
C++に関する質問です
-
2次元配列からダブルポインタへ...
-
指数関数のプログラムを組んで...
-
至急です! マクロ定義で #defi...
-
C 開放してるのにエラー(doubl...
-
C言語★
-
変数の受け渡し double
-
二分探索で方程式の解を求める方法
-
C言語でのプログラム作成について
-
long doubleの表示方法
-
C言語のプログラム作成をお願い...
-
Cプログラムが円ではなく楕円...
-
sqrtを知らないと言われます
-
Cのプログラムに変換
-
定積分のプログラミングについて
-
γ関数のプログラム(初心者です)
-
2次元配列を関数の引数にする方...
-
jacobiのNewton法
-
変数宣言にいらないの入れると...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
float型とdouble型の変数の違い...
-
C言語を実行すると-infが出てき...
-
C 開放してるのにエラー(doubl...
-
c言語で、繰り返し文の中で、0....
-
doubleの変数にintとintの割り...
-
至急です! マクロ定義で #defi...
-
C言語の型による処理速度の違い
-
C言語 関数プロトタイプ宣言の...
-
2次方程式の解を求めるプログ...
-
関数におけるif文とreturn文に...
-
doubleは常に%lfとするべきなのか
-
int とdoubleの比較
-
C言語のプログラムで#include<m...
-
C言語で-23乗を取り扱うには
-
データ数の多い構造体配列
-
指数の表示
-
C言語のpow関数の不具合
-
c言語のプログラミングについて...
-
c言語のコンパイルエラー canno...
おすすめ情報
#include <complex.h>
を入れてみましたが
複素数を1+i,N=5とすると
+1.000e+00 +1.000e+00 i 1.414e+00
・・・
-4.000e+00 -4.000e+00 i 5.656e+00
となってほしいのですが
+6.951e-310 +6.951e-310i +1.732e+00
-6.951e-310 +6.951e-310i +1.732e+00
-0.000e+00 +0.000e+00i +0.000e+00
-0.000e+00 +0.000e+00i +0.000e+00
-0.000e+00 +0.000e+00i +0.000e+00
と表示されます。