跳至內容

Maple/常微分方程

維基教科書,自由的教學讀本

線性常微分方程

[編輯]

dsolve(diff(y(x), x, x) = y(x))

f := diff(y(x), x, x)+a^2*y(x);

Maple:

dsolve(f);

dsolve(diff(y(x), x, x, x, x) = y(x));

y(x) = _C1*exp(x)+_C2*exp(-x)+_C3*sin(x)+_C4*cos(x)


dsolve(diff(y(x), x, x, x, x, x, x, x) = diff(y(x), x, x))

y(x) = _C1+_C2*x+_C3*exp(x)-_C4*exp((-1/4-(1/4)*sqrt(5))*x)*sin((1/4)*sqrt(2)*sqrt(5-sqrt(5))*x)-_C5*exp((-1/4+(1/4)*sqrt(5))*x)*sin((1/4)*sqrt(2)*sqrt(5+sqrt(5))*x)+_C6*exp((-1/4-(1/4)*sqrt(5))*x)*cos((1/4)*sqrt(2)*sqrt(5-sqrt(5))*x)+_C7*exp((-1/4+(1/4)*sqrt(5))*x)*cos((1/4)*sqrt(2)*sqrt(5+sqrt(5))*x)


Maple 式:

dsolve(f);

y(x) = e/d+_C1*exp((1/6)*((36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(2/3)-12*c*a+4*b^2-2*b*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(1/3))*x/(a*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(1/3)))+_C2*exp(-(1/12*I)*(-I*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(2/3)+(12*I)*c*a-(4*I)*b^2-(4*I)*b*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(1/3)+sqrt(3)*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(2/3)+12*sqrt(3)*c*a-4*sqrt(3)*b^2)*x/(a*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(1/3)))+_C3*exp((1/12*I)*(I*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(2/3)-(12*I)*c*a+(4*I)*b^2+(4*I)*b*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(1/3)+sqrt(3)*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(2/3)+12*sqrt(3)*c*a-4*sqrt(3)*b^2)*x/(a*(36*c*b*a-108*d*a^2-8*b^3+12*sqrt(3)*sqrt(4*c^3*a-c^2*b^2-18*c*b*a*d+27*d^2*a^2+4*d*b^3)*a)^(1/3)))

分離變數型微分方程

[編輯]

  • Maple 可以直接得出解答,無需用戶作分離變數步驟。
  • f := diff(y(x), x) = (1+y(x)^2)/(1+x^2);
  • dsolve(f);


  • f := diff(y(x), x) = (1+y(x)^2)*(1+1/x)/(1+x^2)
  • dsolve(f);

歐拉型微分方程

[編輯]


  • f := diff(x(t), t) = (x(t)^2+t^2)/(t*x(t))
  • dsolve(f);
  • ,

  • dsolve(f);

伯努利微分方程

[編輯]

f := diff(x(t), t) = 5*x(t)/t+t*x(t)^(1/2)

  • dsolve(f);

Riccati 型微分方程

[編輯]

  • f := diff(x(t), t) = sin(t)*x(t)^2+2*sin(t)/cos(t)^2
  • dsolve(f);


  • f := diff(x(t), t) = t^5*x(t)^2+5*x(t)/t+6/t^2
  • dsolve(f);

非線性常微分方程

[編輯]

f := (1-x^2)*(diff(y(x), x, x))-2*x*(diff(y(x), x))+l*(l+1);

> dsolve(f);

y(x) = (1/2)*ln(x-1)*_C1+(1/2)*ln(x-1)*l^2+(1/2)*ln(x-1)*l-(1/2)*ln(x+1)*_C1+(1/2)*ln(x+1)*l^2+(1/2)*ln(x+1)*l+_C2


f := x^2*(diff(y(x), x, x))+x*(diff(y(x), x))+(x^2-v^2)*y(x)

dsolve(f);

y(x) = _C1*BesselJ(v, x)+_C2*BesselY(v, x)

f := (1-x^2)*(diff(y(x), x, x))-2*x*(diff(y(x), x))+l(l+1)*y(x);

dsolve(f);

y(x) = _C1*LegendreP((1/2)*sqrt(1+4*l(l+1))-1/2, x)+_C2*LegendreQ((1/2)*sqrt(1+4*l(l+1))-1/2, x)

f := diff(y(x), x, x)-x*(diff(y(x), x))+n*y(x);

dsolve(f);

y(x) = _C1*KummerM(1/2-(1/2)*n, 3/2, (1/2)*x^2)*x+_C2*KummerU(1/2-(1/2)*n, 3/2, (1/2)*x^2)*x

f := x*(diff(y(x), x, x, x))+x*(diff(y(x), x))+y(x)-1

dsolve(f);

y(x) = x*BesselJ(1, x)*_C3+x*BesselY(1, x)*_C2+(1/4)*π*x*(-BesselJ(1, x)*BesselY(0, x)+BesselY(1, x)*BesselJ(0, x))*(Pi*x*_C1*StruveH(-1, x)-2)



y(x) = exp(-(1/4)*x^4)*HeunB(-1/2, 0, -3/2, 1, -(1/2)*x^2)*_C2+exp(-(1/4)*x^4)*HeunB(1/2, 0, -3/2, 1, -(1/2)*x^2)*x*_C1+1

f := diff(y(x), x, x)+(1-x^2)*(diff(y(x), x))+x*y(x)-1

dsolve(f);

y(x) = HeunT(0, 6, -3^(1/3), (1/3)*3^(2/3)*x)*_C2+HeunT(0, -6, -3^(1/3), -(1/3)*3^(2/3)*x)*exp((1/3)*x*(-3+x^2))*_C1+x;


f := diff(y(x), x, x)-cos(x)*y(x)-1



y(x) = MathieuC(0, 2, (1/2)x| _C2 + MathieuS(0, 2, (1/2)x) _C1 - 2 | |


 /        /      1  \\//        /      1  \              /      
 |MathieuS|0, 2, - x|| |MathieuC|0, 2, - x| MathieuSPrime|0, 2, 
 \        \      2  // \        \      2  /              \      
                                                        \ 
 1  \           /      1  \              /      1  \\   | 
 - x| - MathieuS|0, 2, - x| MathieuCPrime|0, 2, - x|| dx| 
 2  /           \      2  /              \      2  //   | 
                                                        / 
                         /  /                                    
         /      1  \     | |  /        /      1  \\//        /   
 MathieuC|0, 2, - x| + 2 | |  |MathieuC|0, 2, - x|| |MathieuC|0, 
         \      2  /     | |  \        \      2  // \        \   
                         \/                                      
    1  \              /      1  \
 2, - x| MathieuSPrime|0, 2, - x|
    2  /              \      2  /
                                                    \          
            /      1  \              /      1  \\   |         /
  - MathieuS|0, 2, - x| MathieuCPrime|0, 2, - x|| dx| MathieuS|
            \      2  /              \      2  //   |         \
                                                    /          
       1  \
 0, 2, - x|
       2  /

>