• ベストアンサー

微分方程式をオイラー法でときたい

高階常微分方程式 y” =f(x、y) y(0)=yo y’(0)=y’o この式と初期値でオイラー法を使って解きたいのですが... オイラー法を二回使えばよい、一階の連立方程式に なおせばよい。という意味がいまいちつかめません。 教えていただけると助かります。お願いします。

質問者が選んだベストアンサー

  • ベストアンサー
  • Tacosan
  • ベストアンサー率23% (3656/15482)
回答No.2

一階の連立方程式に~というのは, 多分 y' = z とおくと (z, y)' = (f(x, y), z), (z(0), y(0)) = (y'_0, y_0) と書けるということでしょう.

aichikyuhaku
質問者

お礼

なるほど! dy/dx=f(x,y) , y(x_0)=y_0 を解くオイラー法の公式 y_k+1=y_k+f(x_k,y_k)Δx があったとすると z'=f(x,y) z(0)=y'_0は z_k+1=z_k+f(x_k,z_k)Δx y'=z y(0)=y_0は y_k+1=y_k+zΔx とすればよいということですか?

その他の回答 (1)

  • saru11
  • ベストアンサー率20% (1/5)
回答No.1

回答と問題の両方の情報をもう少し示していただければ、わかりそうですが・・・。

aichikyuhaku
質問者

お礼

すみませんm(_ _;)m  問題文はこれだけなのです。 二回使えばよいとかは、検索等して調べたのです...

関連するQ&A

  • 常微分方程式の問題で。

    常微分方程式の問題で オイラー法、ルンゲ・クッタ法、ルンゲ・クッタ・ギル法、ミルン法、アダムス法より2つの方法を用いて解け。 y''+y'/x+y=0 初期条件y'/x=-1/2 0<=x<=4 h=0.4 という問題なのですが、基本的に解き方がわかりません。 二階微分なので連立ルンゲ・クッタで解くのかな?っと思ったのですが 結局それからどうすればいいのかわかりません。 2つの方法とありますが、これはどの方法でも解けるものなのでしょうか? 本当はC言語のプログラムを組まなければいけないのですが、 この理論がわからないのでおねがいします。

  • 高階連立常微分方程式の数値計算

    4次のルンゲクッタ法を用いた数値計算を勉強しています. 1階連立常微分方程式と高階常微分方程式は理解でき,プログラムも作成することができました. 次に高階の連立常微分方程式を解こうと思ったら,頭が混乱してしまいました. 4次のルンゲクッタ法を用いて高解連立常微分方程式を解く考え方を教えて頂ければ嬉しいです. また何か良い参考書があれば教えて頂きたいと思います. よろしくお願いします.

  • オイラー法、2次ルンゲクッタ法、4次ルンゲクッタ法のC言語プログラムに

    オイラー法、2次ルンゲクッタ法、4次ルンゲクッタ法のC言語プログラムについて教えてください! 課題なのですが、まったくわからず困ってます>< 1 常微分方程式 dy/dx=f(x,y),y(0)=1 の数値解をオイラー法を用いて計算するプログラムを作為せよ。ただし、f(x,y)=3-6x^2-4x+2xyとする。 2 α=1,β=1,γ=1/2,σ=1/2 の場合の2次ルンゲクッタ法を考える。1と同じ常微分方程式(f(x,y)も同じ)を考え、その数値解を求めるプログラムを作成せよ。また、オイラー法と2次ルンゲクッタ法の実行結果を示して、2つの近似精度を比較せよ。 3 1と同じ常微分方程式(f(x,y)も同じ)を考え、その数値解を4次ルンゲクッタ法を使って求めるプログラムを作成せよ。また、オイラー法、2次ルンゲクッタ法、4次ルンゲクッタ法の実行結果を示して、3つの近似精度を比較せよ。 以上の3つです。 休んでいた自分が悪いのですが、ネットで調べてもよくわからなくて… わかる方、よろしくおねがいします…

  • 常微分方程式、4次のルンゲクッタ法

    (d^2x/dt^2)-2(dy/dt)=f(x) (d^2y/dt^2)+2(dx/dt)=g(y) この連立常微分方程式を4次のルンゲクッタ法で解くためにはどうすればいいのでしょうか?

  • 変分法のオイラーの微分方程式について

    変分法の基本について質問があります。 (1) オイラーの微分方程式は、 d/dx(∂f/∂y')-∂f/∂y=0 となっていますが、これを実際に関数 I=integral(a→b) f(x,y,y')dx にあてはめて解く際に、 第一項では y'だけ考えて y は定数のように扱い、逆に第二項では y だけ考えて y'を定数のように扱うものと思います。 本来yとy'は互いに結びついているはずなのに、なぜこのように独立した変数であるかのように扱えるのでしょうか? (2) また同じくこの第一項を d/dx F_y' と書くとすると、計算するときには F_xy'+F_yy'y' + F_y'y' y'' と全微分のように3つの別れるのはどうしてなのでしょうか? よろしくお願いします。

  • 2階の微分方程式

    こんにちは。現在、微分方程式に取り組んでいます。 おそらくとても単純な所で詰まってしまって、困っています。よろしくお願いします。 式(1) G=dx/dt=py-ax+i 式(2) F=dy/dt=qx-by+j   を使って (1)G=F=0の時、x0、y0を求める (2)n=x-x0,m=y-y0とし、n(t),m(t)それぞれに対する2階微分方程式を求めよ という問題です。 (1)は連立方程式を解いて x=(bi+pj)/(ab-pq)  y=(aj+qi)/(ab-pq) までは出せたのですが、これらをそれぞれx0,y0と考えてしまってよいのでしょうか? 「x,yの2階の微分方程式にする」ようなヒントがあったのですが、そのヒントでかえって混乱しています。 (2)はx0、y0で詰まってしまったので。。。止まっています。 すみませんが、よろしくお願いします。

  • 2階常微分方程式

    x*y''-y'=0 この微分方程式を解きたいのですが解法がわかりません・・・ y"の係数がx^2であればオイラーの微分方程式としてy=x^αを代入して解を得ることができますが この形ですとα=0となってしまい解を求めることができません。 このような定数係数以外の時の2階常微分方程式の解の求め方について説明をお願いします。

  • 微分方程式について

    この連立微分方程式の解き方がわかりません。どのようにして解けばいいのでしょうか。  x''=Ax/r^2+y'B  y''=Ay/r^2-x'B (A,Bは定数 x'',y''はそれぞれx,yの二階微分、x,yはそれぞれx,yの一階微分 r=(x^2+y^2)^1/2) どなたかよろしくお願いします。

  • 常微分方程式の数値計算

    実験で常微分方程式の数値計算をしました y'(x)=2xy,y(0)=1.0(解はy=ex2)←eのx二乗の微分方程式できざみ幅h=0.1、区間〔0,2〕の条件で オイラー法と改良オイラー法とルンゲクッタ法を使って 表計算ソフトを使って求めよとあるのですが、 表計算で近似値のy座標とか Δt後のxの変化分K1とか求めるのですが、全然意味がわからなくて、 何かいい方法はないでしょうか?

  • 常微分方程式の問題

    微分方程式x'=x^2-t-1で初期値がx(0)=1の問題で、 オイラー法とルンゲクッタ法は出せたのですが、 それらと比較するために出す解析解がわかりません。 どなたかよろしくお願いします。