>> y=solve( ' D2y-2Dy-3*y=0 ' ， ' y(0)=0，y(1)=1 ' )

y=

       1/(exp(-1)-exp(3))*exp(-x)-1/(exp(-1)-exp(3))*exp(3*x)

>> y=simple(y)  %  this looks like a candidate for simplification

y=

       -(exp(-x)-exp(3*x))/(exp(3)-exp(-1))

>> pretty(y)  %   pretty it up

           exp(-x)-exp(3 x)

       - ---------------------

          exp(3) -exp(-1)

现在来绘制感兴趣的区域内的结果。

>> ezplot(y，[-6 2])

图22.3   符号函数y=-(exp(-x)-exp(3*x))/(exp(3)-exp(-1))  (-6≤x≤2)

微分方程组

函数dsolve也可同时处理若干个微分方程式，下面有两个线性一阶方程。

              =3f+4g            =-4f+3g