c言語の課題について

Question

前提・実現したいこと

学校の課題で1階常微分方程式 y’ = e^x − y の解をオイラー法、ホイン法、ルンゲ・クッタ法を用いて求めるという課題が出ました。プログラム自体はかけたのですが、下記の発生している問題を直したうえで、もう少し簡潔にできるところは直したいです。

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

最初xの値に対してその次の瞬間のxの値で求められるyの値が出てくるようになってしまいました。

該当のソースコード

#include <stdio.h>
#include <math.h>
#define X 0.0
#define Y 0.0
#define N 10.0

double func(double x,double y);
double euler(double x,double y,double h);
double heun(double x,double y,double h);
double runge_kutta(double x,double y,double h);

int main(void){
	double x,y,h = 1/N;
	
	printf("euler\n");
	x = X;
	y = Y;
	for( ;x<=1.0;){
		y = euler(x,y,h);
		x = x + h;
		printf("x = %f,y = %f\n",x,y);
	}
	
	printf("heun\n");
	x = X;
	y = Y;
	for( ;x<=1.0;){
		y = heun(x,y,h);
		x = x + h;
		printf("x = %f,y = %f\n",x,y);
	}
	
	printf("runge_kutta\n");
	x = X;
	y = Y;
	for( ;x<=1.0;){
		y = runge_kutta(x,y,h);
		x = x + h;
		printf("x = %f,y = %f\n",x,y);
	}
	return 0;
}

double func(double x,double y){
	
	return exp(x) - y;
}

double euler(double x,double y,double h){
	return y + h * func(x,y);
}

double heun(double x,double y,double h){
	double k_1,k_2;
	k_1 = h * func(x,y);
	k_2 = h * func(x+h,y+k_1);
	return y+(k_1+k_2)/2;
}

double runge_kutta(double x,double y,double h){
	double k_1,k_2,k_3,k_4;
	k_1 = h * func(x,y);
	k_2 = h * func(x+h/2,y+k_1/2);
	k_3 = h * func(x+h/2,y+k_2/2);
	k_4 = h * func(x+h,y+k_3);
	return y+(k_1+2*k_2+2*k_3+k_4)/6;
}

試したこと

for文の変化式をx += hにしておいたら問題が発生したのでyの値を求めた後にxの値を更新するようにしました。

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

ここにより詳細な情報を記載してください。

Comments

[y_waiwai]

どういう不具合なのかもちっと具体的に説明してください。
正常ならそれがどうなるのかも含めて

Answer 1

最初xの値に対してその次の瞬間のxの値で求められるyの値が出てくるようになってしまいました。

実行すると出てくる結果と期待していた結果を具体的に質問に書くと、やりたいことを明確に伝えやすいと思います。質問を編集して書いてください。

「最初xの値に対して」とは x == 0 のときのことだと思うので、 printf の位置を直してみました。
これは意図通りになっているでしょうか？

    printf("euler\n");
    x = X;
    y = Y;
    for( ;x<=1.0;){
        printf("x = %f,y = %f\n",x,y);
        y = euler(x,y,h);
        x = x + h;
    }
    printf("x = %f,y = %f\n",x,y);

Comments

[退会済みユーザー]

質問内容に言葉が足らなかった部分があってすみませんでした
考えていた通りになりました
ありがとうございます