Solving Linear Ordinary Differential Equations (LODEs) by computing integrating factors for them

はじめる前に
新機能一覧
Maple ワークシートの作成
Mapleワークシートを共有
Maple ウィンドウのカスタマイズ
Maple システムのカスタマイズ
コネクティビティ
Mathematics
数学
- テンソル解析
- パッケージ
- ファンクションアドバイザ
- ベキ級数
- ベクトル解析
- 不活性関数
- 代数
- 因数分解と方程式解法
- 基本数学
- 基本的な注意事項
- 変分法
- 変換
- 幾何
- 微分代数方程式
- 微分幾何学
- 微分方程式
  - dsolve コマンド
    - numeric (数値)
    - アルゴリズム
      - 概要
      - 1 階 ODE
      - Abel
      - 因数分解
      - 完全に正確
      - 対称パターン
      - 教育
      - 積分因子 (1 次)
      - 積分因子 (2 次)
      - 線形
      - 超幾何学
    - 基本概要
    - 詳細概要
    - Abel
    - DESol
    - inttrans
    - Lie
    - ODESolStruc
    - series コマンド
    - システム
    - 初期条件
    - 区分関数
    - 参照
    - 対称パターン
    - 対話型
    - 形式的冪級数
    - 形式解
    - 微分因数分解
    - 教育
    - 積分因子 (1 次)
    - 積分因子 (2 次)
    - 設定
    - 超幾何学
  - DEtools パッケージ
  - PDEtools パッケージ
  - pdsolve コマンド
  - Rif パッケージ
  - Slode パッケージ
  - リー対称群による解法
  - 例題
  - 常微分方程式（ODE）の分類
  - 微分代数
  - 微分代数（diffalg）パッケージ
  - 微分形式
  - 偏微分方程式の解関数の確認（pdetest コマンド）
  - 初期条件
  - 常微分方程式の解関数の確認（odetest コマンド）
  - 常微分方程式用対話型 Maplet インターフェイス
- 微積分
- 数
- 数値計算
- 数学関数
- 数論
- 最適化
- 特殊関数
- 統計
- 線形代数
- 群論
- 評価
- 論理
- 過去のパッケージ
- 金融
- 離散数学
- piecewise examples
- グラフ電卓
- 区分関数
物理
プログラミング
グラフィックス
Student パッケージ
Science and Engineering
科学・エンジニアリング
アプリケーションと例題
リファレンス
システム
MapleSim
MapleSim ツールボックス
MapleSim ヘルプ
Tasks
Toolboxes
エラーメッセージガイド
マニュアル
数学アプリ

Description

Integrating factors for second and higher order linear ODEs

•	For linear ODEs (LODEs) of order 2 or greater, it is possible to calculate integrating factors by solving the adjoint of the LODE. This could be as difficult as the original problem, or much easier, depending on the example. This method is implemented in dsolve.

Examples

ode := diff(y(x), x, x) = (y(x)*x+y(x)*ln(x)*x^2+(diff(y(x), x))*ln(x)*x^2)*exp(x)/(ln(x)*x^2)-y(x)/(ln(x)*x^2)

ode := diff(y(x), `$`(x, 2)) = (y(x)*x+y(x)*ln(x)*x^2+(diff(y(x), x))*ln(x)*x^2)*exp(x)/(ln(x)*x^2)-y(x)/(ln(x)*x^2)

(1)

This ODE has the following adjoint

-(exp(x)*x-1)*y(x)/(ln(x)*x^2)+exp(x)*(diff(y(x), x))+diff(y(x), `$`(x, 2))

(2)

This adjoint equation is in turn solvable by dsolve

y(x) = ((Int(exp(-exp(x))/ln(x)^2, x))*_C1+_C2)*ln(x)

(3)

Now the solutions to the adjoint equation are integrating factors of the original LODE, so the two independent solutions implied in the general solution above

(4)

(Int(exp(-exp(x))/ln(x)^2, x)+2)*ln(x)

(5)

are integrating factors of ode. These integrating factors could also be found using the intfactor directly

Mu := ln(x), (Int(exp(-exp(x))/ln(x)^2, x))*ln(x)

(6)

Constructing solutions using integrating factors

How are these integrating factors transformed into a solution to the original problem? By using them to construct two first integrals; that is: two ODEs of lower order (in this case two first order ODEs). For that purpose it is provided the firint command which receives an exact ODE and returns a first integral. The idea is simple: an exact ODE is a total derivative - say dR/dx; firint returns the R + _C1:

-(ln(x)*x^2*exp(x)+x)*y(x)/x^2+ln(x)*(diff(y(x), x))+_C1 = 0

(7)

-((Int(exp(-exp(x))/ln(x)^2, x))*ln(x)^2*exp(x)*x^2+exp(-exp(x))*x^2+(Int(exp(-exp(x))/ln(x)^2, x))*ln(x)*x)*y(x)/(ln(x)*x^2)+(Int(exp(-exp(x))/ln(x)^2, x))*ln(x)*(diff(y(x), x))+_C1 = 0

(8)

Eliminating y' from these two first integrals (and replacing _C1 by _C2 in one of them) leads to the solution f(x,y(x),_C1,_C2) = 0 to this ode. So this process could be run interactively, as shown, or in one step