プログラマーの気分になれるタイピング

(％．３ｅ％．３ｅ％．３ｅ￥ｎ”，)

%.3e %.3e %.3e\n",

(ｔ，ｘ，ｘ１，ｘ２，ｘ３，ｅ１，ｅ２，ｅ３）；)

t, x, x1, x2, x3, e1, e2, e3);

(ｘ＝ｅｘｐ（－（ｔ＋ｈ）＊（ｔ＋ｈ）／２．０）＋１；)

x=exp(-(t+h)*(t+h)/2.0)+1;

(／＊おいらーほう＊／)

/*オイラー法*/

(ｘ１＝（ｔ－ｔ＊ｘ１）＊ｈ＋ｘ１；)

x1=(t-t*x1)*h+x1;

(ｅ１＝ｆａｂｓ（ｘ１－ｘ）；)

e1=fabs(x1-x);

(／＊ほいんほう＊／)

/*ホイン法*/

(ｋ１＝ｔ－ｔ＊ｘ２；)

k1=t-t*x2;

(ｘｎ＝ｋ１＊ｈ＋ｘ２；)

xn=k1*h+x2;

(ｋ２＝ｔ＋ｈ－（ｔ＋ｈ）＊ｘｎ；)

k2=t+h-(t+h)*xn;

(ｘ２＝（ｋ１＋ｋ２）＊ｈ／２．０＋ｘ２；)

x2=(k1+k2)*h/2.0+x2;

(ｅ２＝ｆａｂｓ（ｘ２－ｘ）；)

e2=fabs(x2-x);

(／＊るーげんくったほう＊／)

/*ルーゲンクッタ法*/

(ｋ１＝ｔ－ｔ＊ｘ３；)

k1=t-t*x3;

(ｘｎ＝ｋ１＊ｈ／２．０＋ｘ３；)

xn=k1*h/2.0+x3;

(ｋ２＝ｔ＋ｈ／２．０－（ｔ＋ｈ／２．０）＊ｘｎ；)

k2=t+h/2.0-(t+h/2.0)*xn;

(ｘｎ＝ｋ２＊ｈ／２．０＋ｘ３；)

xn=k2*h/2.0+x3;

(ｋ３＝ｔ＋ｈ／２．０－（ｔ＋ｈ／２．０）＊ｘｎ；)

k3=t+h/2.0-(t+h/2.0)*xn;

(ｘｎ＝ｋ３＊ｈ＋ｘ３；)

xn=k3*h+x3;

(ｋ４＝ｔ＋ｈ－（ｔ＋ｈ）＊ｘｎ；)

k4=t+h-(t+h)*xn;

(ｘ３＝（ｋ１＋２＊ｋ２＋２＊ｋ３＋ｋ４）／６．０＊ｈ＋ｘ３；)

x3=(k1+2*k2+2*k3+k4)/6.0*h+x3;

(ｅ３＝ｆａｂｓ（ｘ３－ｘ）；)

e3=fabs(x3-x);

(｝)

}

(ｒｅｔｕｒｎ（０）；)

return(0);

(｝)

}