二次方程式の解を求めるプログラムで、実行時に"Segmentation fault"と表示されます。

Question

#前提・実現したいこと

　二次方程式f(x)=(ax - b)^2 + (cx - d)^2=0の解xを解の公式を用いて、求めるプログラムを記述しています。条件は以下の通りです。

・引数yの平方根を出力するmysqrt関数を定義すること
⇒関数頭部をdouble mysqrt(double y)とする
⇒結果に至るまで、計算途中結果と反復回数を表示すること
⇒反復回数が1000回を超えないこと

・<math.h>をインクルードしないこと(sqrt等使用しない)

・a=0の場合には、メッセージを表示し、正常終了すること
⇒<stdlib.h>で定義される標準ライブラリ関数exitを使用すること

・判別式の評価を行い、分岐処理により場合分けを行うこと
⇒D<0,D=0,D<0

##発生している問題・エラーメッセージ
虚数解の場合が表示されません。すべて解がb/aとして計算されてしまいます。

###該当のソースコード

//修正済み

#include <stdio.h>
#include <stdlib.h>


//解の公式の√(判別式)の部分の処理を行う関数

double mysqrt(double y){
	
	double a,b,c,d;
	
	int count=0; 				//反復回数
		
	printf(" %21.16g（反復回数%d 回）¥n", y, count);  
    
	while(count < 1000) {
    	count++;
    	y = -(a*d + b*c)*(a*d + b*c); //判別式
    
    	printf(" %21.16g（反復回数%d 回）¥n", y, count);  
    	
    }
    
    //値の更新
    y = -(a*d + b*c)*(a*d + b*c);
    
    return y;
}

int main(){
	
	double a=0,b=0,c=0,d=0;
	
	double num1=0,num2=0; 			//解
	double real=0, image=0;			//複素数を含む場合
	
	double y=0;		//判別式
	
	printf("係数を入力せよ\n");
	
	scanf("%lf",&a);
	scanf("%lf",&b);
	scanf("%lf",&c);
	scanf("%lf",&d);
	
	
	//入力データの確認
	if ((a*a+c*c)==0) {  
			printf("解の公式の分母は0以外の数にしてください¥n");  
			exit(EXIT_SUCCESS);  
	} 
	
	
	if(y==0){
		num1 = num2 = b/a;  //重解を持つ場合
		
		printf("解: x=%.2f, %.2f\n",num1,num2);
		
	}else{                      //虚数解をもつ場合
		real = (a*b+c*d)/(a*a+c*c);
		image = mysqrt(y)/(a*a+c*c);
		
		printf("解: %.2f+%.2fi, %.2f-%.2fi\n",real,image);
	}
	return 0;
}

試したこと

・判別式Dをyとし、その中身をある程度整理した形にした
・全ての変数を初期化
・解を求めるif文の箇所を訂正

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

Raspberry Pi 3A+を使用

Comments

[majiponi]

質問文は修正できます。まだ問題があるようでしたら、修正した結果としてのコードと、どのような問題があるかの記載をお願いいたします。

[退会済みユーザー]

修正しましたが、解がすべてb/a(重解の場合)として計算されてしまいます。

Answer 1

if (!(aa+cc)==0) {

ここ，おかしくないですか．
「0のときはダメ」ということをやりたいなら，それをそのまま素直に書けばよいのでは．

---

すべて解がb/aとして計算されてしまいます

その処理に入る条件は「yが0なら」となっています．
あなたのコードではyの値は0で初期化されていて，この判定の箇所までそのままですから，常に条件を満たします．

Comments

[退会済みユーザー]

ご回答ありがとうございます。ご指摘の箇所、訂正しました。

[fana]

一応回答に追記はしときましたが，

１点指摘されたら質問修正して別の問題が発覚して…ってのを「完成」まで繰り返すつもりなんでしょうか？
そういうやり方はこの場にふさわしくないと思いますし，
それ以前に，「自身で取り組む気が全く無いならば，やらなければいいのに」と感じてしまいます．

[退会済みユーザー]

分かりました。今後は気を付けます。

Answer 2

(ax - b)^2 + (cx - d)^2=0 を解くのよね？

(ax - b)^2 >= 0, (cx - d)^2 >= 0 だから
左辺が0になるのは　(ax - b)^2 = 0, (cx - d)^2 = 0 のとき。
つまり x = b/a かつ x = d/c
b/a != d/c のとき解はない。

Comments

[退会済みユーザー]

ご回答ありがとうございます。
ご指摘の通り、計算が誤っていました。修正します。

[episteme]

なんで低評価? なんか間違ったコト言ってる?

[退会済みユーザー]

私はまだ低評価を付ける条件を満たしていないので、そのような評価は付けられません。
episteme様の回答、とても参考になりました。ありがとうございます。

Answer 3

気づいた点だけ指摘、代入になっています。

    if (!(a=0)) {   // 誤
    if (!(a==0)) {  // 正

Comments

[退会済みユーザー]

ご回答ありがとうございます。
ご指摘の箇所、修正いたしました。

Answer 4

scanf("%lf",a);

scanfの引数はアドレスですよ

Comments

[退会済みユーザー]

ご回答ありがとうございます。
該当の部分、修正しました。"Segmentation fault"とは表示されなくなりましたが、aの数値が0以外にも関わらず、正常終了してしまいます。また、入力データの確認の部分の"a=0"を"a==0"と修正しました。

[maisumakun]

> aの数値が0以外にも関わらず、正常終了してしまいます

「2次方程式を解くプログラム」なのに、「2次方程式を入力して正常終了する」のが問題なのですか？

[majiponi]

質問。これ、問題文合ってますか? みた感じ、ad=bcのときに重解b/aを持ち、それ以外は虚数解になるんですが…。判別式も、文字通り問題文を解釈すると、y=-(ad+bc)^2になりますよ?

[退会済みユーザー]

その部分は、解の方程式の分母がもし0であれば正常終了するという記述になっています。"a"というのは、二次方程式a*x^2+b*x+cの解の公式における分母"2*a"を指します。
そこで今間違いに気づいたのですが、入力データの確認の部分は「aの数値」ではなく、「a^2+c^2の値」(f(x)の解の公式における分母）正しいです。説明不足で申し訳ありません。

[退会済みユーザー]

majiponi様
問題文は合っていますが、ご指摘の通り判別式の計算が間違っていました。解も重解と虚数解になりますね。訂正しておきます。