The ODE Analyzer, a new interactive interface, has been created for numeric and symbolic solutions of ODE and ODE systems. The interface is easy to use, with many common options. In addition to computation of solutions and solution values, plots can also be generated. The interface is suitable as an educational tool; it can provide the corresponding Maple commands needed to compute the solution or plot. For information on launching the interface (and command-line options), see dsolve/interactive. For an illustrated description of the interface, see worksheet/interactive/dsolve.

First Order Nonlinear Equations

New solving methods for ODEs of Abel type.

New solving methods for ODEs that can be mapped into solvable Abel and Riccati type ODEs.

Second Order Linear Equations

New solving methods for ODEs that have non-rational coefficients or admit hypergeometric solutions.

New solving methods for rational coefficient ODEs that admit doubly-periodic non-Liouvillian solutions.

New integration methods for second order linear ODEs that explore non-standard forms of their Lie symmetries.

Miscellaneous

Update of many dsolve routines to return ODE solutions in terms of hypergeometric functions instead of integrals if possible.

Significant optimization of dsolve routines for the computation and use of differential invariants to solve nonlinear ODEs.

Significant optimization of dsolve routines by more quickly abandoning attempts to solve an ODE using algorithms that cannot solve it. For average problems, routines exit three to ten times faster.

Update of the dsolve function to use the LinearFunctionalSystems package for the computation of series solutions to linear ODE systems with polynomial coefficients.

New solving methods for PDE systems and optimization of the existing pdsolve routines for solving a single PDE.

Numerical Methods for Ordinary Differential Equations

More efficient Taylor series ODE solutions with high precision.

Interactive features added to the dsolve function output. You can query about and change the initial data, method used, and other parameters.

Ability to solve ODE systems containing user-defined functions.

Numerical Methods for Partial Differential Equations

Extension of all numeric PDE methods to solve nonlinear PDEs that are linear in the highest time derivative.

Addition of time-adaptive solution computation. For example, you can now compute a solution to within a specified error per unit time step.

Computation and plotting of error estimates for a solution as a function of time and space.

PDEtools[declare](y(x), prime=x);

$y\left(x\right)\mathrm{will\; now\; be\; displayed\; as}y$

$\mathrm{derivatives\; with\; respect\; to}x\mathrm{of\; functions\; of\; one\; variable\; will\; now\; be\; displayed\; with\; \text{'}}$

Solutions for more third order linear ODE families that depend on arbitrary functions:

diff(y(x),x,x,x) = g(x)*diff(g(x),x)*y(x) + 2*diff(g(x),x)*diff(y(x),x) + diff(g(x),x,x)*y(x) + g(x)^2*diff(y(x),x);

$\mathrm{y\text{'}\text{'}\text{'}}=g\left(x\right)\mathrm{g\text{'}}y+2\mathrm{g\text{'}}\mathrm{y\text{'}}+\mathrm{g\text{'}\text{'}}y+{g\left(x\right)}^2\mathrm{y\text{'}}$

dsolve((2),y(x));

$y=\mathrm{c\_\_1}{\ⅇ}^{\int g\left(x\right)\ⅆx}+\mathrm{c\_\_2}{\left(\int {\ⅇ}^{-\left(\int g\left(x\right)\ⅆx\right)}\ⅆx\right)}^2{\ⅇ}^{\int g\left(x\right)\ⅆx}+\mathrm{c\_\_3}{\ⅇ}^{\int g\left(x\right)\ⅆx}\left(\int {\ⅇ}^{-\left(\int g\left(x\right)\ⅆx\right)}\ⅆx\right)$

(1-I)*diff(y(x),x,x,x,x)-(3+5*I)*diff(y(x),x,x,x)+(-I*x+x-9+3*I)*diff(y(x),x,x)+(8*I-3*I*x-x)*diff(y(x),x)+(2+2*I-I*x-x)*y(x);

$\left(1-\mathrm{I}\right)\mathrm{y\text{'}\text{'}\text{'}\text{'}}-\left(3+5\mathrm{I}\right)\mathrm{y\text{'}\text{'}\text{'}}+\left(-\mathrm{I}x+x-9+3\mathrm{I}\right)\mathrm{y\text{'}\text{'}}+\left(8\mathrm{I}-3\mathrm{I}x-x\right)\mathrm{y\text{'}}+\left(2+2\mathrm{I}-\mathrm{I}x-x\right)y$

dsolve((4),y(x));

$y=\mathrm{c\_\_1}{\ⅇ}^{\frac{\left(-1+2\mathrm{I}+\mathrm{I}\sqrt{3}\right)x}{2}}+\mathrm{c\_\_2}{\ⅇ}^{-\frac{\left(\mathrm{I}\sqrt{3}+1-2\mathrm{I}\right)x}{2}}+\mathrm{c\_\_3}\left(-\left(\int \mathrm{AiryAi}\left(1-x\right){\ⅇ}^{-\frac{x\left(\mathrm{I}\sqrt{3}-1\right)}{2}}\ⅆx\right){\ⅇ}^{\frac{\left(-1+2\mathrm{I}+\mathrm{I}\sqrt{3}\right)x}{2}}+\left(\int \mathrm{AiryAi}\left(1-x\right){\ⅇ}^{\frac{x\left(1+\mathrm{I}\sqrt{3}\right)}{2}}\ⅆx\right){\ⅇ}^{-\frac{\left(\mathrm{I}\sqrt{3}+1-2\mathrm{I}\right)x}{2}}\right)+\mathrm{c\_\_4}\left(-\left(\int \mathrm{AiryBi}\left(1-x\right){\ⅇ}^{-\frac{x\left(\mathrm{I}\sqrt{3}-1\right)}{2}}\ⅆx\right){\ⅇ}^{\frac{\left(-1+2\mathrm{I}+\mathrm{I}\sqrt{3}\right)x}{2}}+\left(\int \mathrm{AiryBi}\left(1-x\right){\ⅇ}^{\frac{x\left(1+\mathrm{I}\sqrt{3}\right)}{2}}\ⅆx\right){\ⅇ}^{-\frac{\left(\mathrm{I}\sqrt{3}+1-2\mathrm{I}\right)x}{2}}\right)$

Solutions for second order linear ODEs in terms of elementary, hypergeometric, or special functions instead of uncomputable integrals:

-4*x*y(x) + diff(y(x),x) - 2*diff(y(x),x,x)*x + 2*diff(y(x),x,x)*x^3;

$-4yx+\mathrm{y\text{'}}-2\mathrm{y\text{'}\text{'}}x+2\mathrm{y\text{'}\text{'}}{x}^3$

dsolve((6));

diff(y(x),x,x) = cos(2*x)/sin(2*x)*diff(y(x),x)-2*y(x);

$\mathrm{y\text{'}\text{'}}=\frac{\cos \left(2x\right)\mathrm{y\text{'}}}{\sin \left(2x\right)}-2y$

dsolve((8));

$y=\mathrm{c\_\_1}{\sin \left(2x\right)}^{\raisebox{1ex}{$3$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}\mathrm{LegendreP}\left(\frac{1}{4}\,\frac{3}{4}\,\cos \left(2x\right)\right)+\mathrm{c\_\_2}{\sin \left(2x\right)}^{\raisebox{1ex}{$3$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}\mathrm{LegendreQ}\left(\frac{1}{4}\,\frac{3}{4}\,\cos \left(2x\right)\right)$

The solutions for Kamke's examples 43, 133, 263, 319, 353, 354, and 355 had uncomputable integrals in previous releases of Maple. In Maple 9, they are solved in terms of hypergeometric functions.

New elementary and hypergeometric solutions for linear ODEs families with non-rational coefficients:

diff(y(x),x,x) = 1/2*(2*x^3*y(x)*(x^2-1)^(1/2)-x^3*y(x)+diff(y(x),x)*x^4-3*diff(y(x),x)*x^2+2*diff(y(x),x))/x/(x^2-1)^2;

dsolve((10));

diff(y(x),x,x) = (-x^5*y(x)-3*x^3*y(x)+3*diff(y(x),x)*(x^2+3)^(1/2))/(x^2+3)^(3/2)/x;

$\mathrm{y\text{'}\text{'}}=\frac{-x^{5}y-3x^{3}y+3\mathrm{y\text{'}}\sqrt{x^2+3}}{{\left(x^2+3\right)}^{\raisebox{1ex}{$3$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}x}$

dsolve((12));

New doubly-periodic solutions for second order linear ODEs; in this case, the solution is computed directly from the ODE symmetries.

diff(y(x),x,x)-12*WeierstrassP(x,a,b)*y(x);

$\mathrm{y\text{'}\text{'}}-12\mathrm{WeierstrassP}\left(x\,a\,b\right)y$

dsolve((14));

$y=\mathrm{c\_\_1}\mathrm{WeierstrassPPrime}\left(x\,a\,b\right)+\mathrm{c\_\_2}\mathrm{WeierstrassPPrime}\left(x\,a\,b\right)\left(-\frac{2\left(6a{\mathrm{WeierstrassP}\left(x\,a\,b\right)}^2-9\mathrm{WeierstrassP}\left(x\,a\,b\right)b-a^2\right)\mathrm{WeierstrassPPrime}\left(x\,a\,b\right)}{4{\mathrm{WeierstrassP}\left(x\,a\,b\right)}^3-\mathrm{WeierstrassP}\left(x\,a\,b\right)a-b}-6a\mathrm{WeierstrassZeta}\left(x\,a\,b\right)+9bx\right)$

(WeierstrassPPrime(x,a,b) + WeierstrassP(x,a,b)^2)*diff(diff(y(x),x),x) + (WeierstrassP(x,a,b)^3 - WeierstrassP(x,a,b)*WeierstrassPPrime(x,a,b) - 6*WeierstrassP(x,a,b)^2 + 1/2*a)*diff(y(x),x) + (WeierstrassPPrime(x,a,b)^2 - WeierstrassP(x,a,b)^2*WeierstrassPPrime(x,a,b) - WeierstrassP(x,a,b)*(6*WeierstrassP(x,a,b)^2 - 1/2*a))*y(x) = 0;

$\left(\mathrm{WeierstrassPPrime}\left(x\,a\,b\right)+{\mathrm{WeierstrassP}\left(x\,a\,b\right)}^2\right)\mathrm{y\text{'}\text{'}}+\left({\mathrm{WeierstrassP}\left(x\,a\,b\right)}^3-\mathrm{WeierstrassP}\left(x\,a\,b\right)\mathrm{WeierstrassPPrime}\left(x\,a\,b\right)-6{\mathrm{WeierstrassP}\left(x\,a\,b\right)}^2+\frac{a}{2}\right)\mathrm{y\text{'}}+\left({\mathrm{WeierstrassPPrime}\left(x\,a\,b\right)}^2-{\mathrm{WeierstrassP}\left(x\,a\,b\right)}^2\mathrm{WeierstrassPPrime}\left(x\,a\,b\right)-\mathrm{WeierstrassP}\left(x\,a\,b\right)\left(6{\mathrm{WeierstrassP}\left(x\,a\,b\right)}^2-\frac{a}{2}\right)\right)y=0$

dsolve((16));

$y=\mathrm{c\_\_1}\mathrm{WeierstrassP}\left(x\,a\,b\right)+\mathrm{c\_\_2}{\ⅇ}^{\mathrm{WeierstrassZeta}\left(x\,a\,b\right)}$

New solutions for linear ODE problems that do not match any classification, but for which it is possible to compute symmetries of a non-standard form that can be used to determine a solution:

diff(diff(y(x),x),x) = 1/4*(-2*Psi(x)*Psi(2,x)+3*Psi(1,x)^2+Psi(x)^4)/Psi(x)^2*y(x);

$\mathrm{y\text{'}\text{'}}=\frac{\left(-2\mathrm{\Psi }\left(x\right)\mathrm{\Psi }\left(2\,x\right)+3{\mathrm{\Psi }\left(1\,x\right)}^2+{\mathrm{\Psi }\left(x\right)}^4\right)y}{4{\mathrm{\Psi }\left(x\right)}^2}$

dsolve((18));

$y=\frac{\mathrm{c\_\_1}{\ⅇ}^{-\frac{\mathrm{\Psi }\left(\mathrm{-1}\,x\right)}{2}}}{\sqrt{\mathrm{\Psi }\left(x\right)}}+\frac{\mathrm{c\_\_2}{\ⅇ}^{\frac{\mathrm{\Psi }\left(\mathrm{-1}\,x\right)}{2}}}{\sqrt{\mathrm{\Psi }\left(x\right)}}$

Simpler solutions in terms of elementary and special functions instead of hypergeometric functions:

diff(diff(y(x),x),x) = 1/4*(4*(x^n)^2*n^2-1+4*(x^n)^4*n^2+n^2)/x^2*y(x);

$\mathrm{y\text{'}\text{'}}=\frac{\left(4{\left(x^n\right)}^2{n}^2-1+4{\left(x^n\right)}^4{n}^2+n^2\right)y}{4x^2}$

dsolve((20));

$y=\mathrm{c\_\_1}{x}^{\frac{1}{2}-\frac{n}{2}}{\ⅇ}^{\frac{x^{2n}}{2}}+\mathrm{c\_\_2}{x}^{\frac{1}{2}-\frac{n}{2}}\mathrm{erf}\left(x^n\right){\ⅇ}^{\frac{x^{2n}}{2}}$

Simpler solution in terms of square roots instead of compositions of sin, cos, ln, and radicals:

-1/4*1/x/(x-1)*y(x)+1/2*(2*x-1)/x/(x-1)*diff(y(x),x)+diff(y(x),x,x) = 0;

$-\frac{y}{4x\left(-1+x\right)}+\frac{\left(2x-1\right)\mathrm{y\text{'}}}{2x\left(-1+x\right)}+\mathrm{y\text{'}\text{'}}=0$

dsolve((22));

$y=\mathrm{c\_\_1}\sqrt{x}+\mathrm{c\_\_2}\sqrt{-1+x}$

Quickly abandon attempt to solve an ODE for which no solution can be computed:

x^n*diff(y(x),x,x)+c*(a*x+b)^(n-4)*y(x)=0;

$x^n\mathrm{y\text{'}\text{'}}+c{\left(ax+b\right)}^{n-4}y=0$

dsolve((24));

$y=\mathrm{DESol}\left(\left\{\mathrm{\_Y\text{'}\text{'}}+\frac{c{\left(ax+b\right)}^{n-4}\mathrm{\_Y}\left(x\right)}{x^n}\right\}\,\left\{\mathrm{\_Y}\left(x\right)\right\}\right)$

Improved recognition of Elliptic function solutions:

diff(y(x),x)^2 = 1/(1+x^4*k^2+(-k^2-1)*x^2);

${\mathrm{y\text{'}}}^2=\frac{1}{1+x^4{k}^2+\left(-k^2-1\right)x^2}$

dsolve((26));

$y=\mathrm{EllipticF}\left(x\,k\right)+\mathrm{c\_\_1},y=-\mathrm{EllipticF}\left(x\,k\right)+\mathrm{c\_\_1}$

Kamke's example 80, in terms of an arbitrary function $f\left(x\right)$:

diff(y(x),x)+f(x)*sin(y(x))+(1-diff(f(x),x))*cos(y(x))-diff(f(x),x)-1, y(x);

$\mathrm{y\text{'}}+f\left(x\right)\sin \left(y\right)+\left(1-\mathrm{f\; \text{'}}\right)\cos \left(y\right)-\mathrm{f\; \text{'}}-1,y$

dsolve((28));

$y=2\arctan \left(\frac{-{\ⅇ}^{\int f\left(x\right)\ⅆx}+\left(\int {\ⅇ}^{\int f\left(x\right)\ⅆx}\ⅆx\right)f\left(x\right)+f\left(x\right)\mathrm{c\_\_1}}{\mathrm{c\_\_1}+\int {\ⅇ}^{\int f\left(x\right)\ⅆx}\ⅆx}\right)$

Implicit hypergeometric solution for Kamke's example 81 and a more general ODE class, of which example 81 is a special case:

diff(y(x),x) = -2*tan(x)*tan(y(x))+1;

$\mathrm{y\text{'}}=-2\tan \left(x\right)\tan \left(y\right)+1$

dsolve((30));

$\mathrm{c\_\_1}+\frac{\tan \left(x\right)}{{\left(\frac{\left(1+{\tan \left(y\right)}^2\right)\left(1+{\tan \left(x\right)}^2\right)}{{\left(\tan \left(x\right)\tan \left(y\right)-1\right)}^2}\right)}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}}+\frac{\left(\tan \left(y\right)+\tan \left(x\right)\right)\mathrm{hypergeom}\left(\left[\frac{1}{2}\,\frac{5}{4}\right]\,\left[\frac{3}{2}\right]\,-\frac{{\left(\tan \left(y\right)+\tan \left(x\right)\right)}^2}{{\left(\tan \left(x\right)\tan \left(y\right)-1\right)}^2}\right)}{2\left(\tan \left(x\right)\tan \left(y\right)-1\right)}=0$

diff(y(x),x) = -(P(y(x))^2+1)*(2*F(x)*P(y(x))-1)*diff(F(x),x)/(1+F(x)^2)/D(P)(y(x));

$\mathrm{y\text{'}}=-\frac{\left({P\left(y\right)}^2+1\right)\left(2F\left(x\right)P\left(y\right)-1\right)\mathrm{F\text{'}}}{\left(1+{F\left(x\right)}^2\right)\mathrm{D}\left(P\right)\left(y\right)}$

dsolve((32), y(x));

$\mathrm{c\_\_1}+\frac{F\left(x\right)}{{\left(\frac{\left({P\left(y\right)}^2+1\right)\left(1+{F\left(x\right)}^2\right)}{{\left(F\left(x\right)P\left(y\right)-1\right)}^2}\right)}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}}+\frac{\left(P\left(y\right)+F\left(x\right)\right)\mathrm{hypergeom}\left(\left[\frac{1}{2}\,\frac{5}{4}\right]\,\left[\frac{3}{2}\right]\,-\frac{{\left(P\left(y\right)+F\left(x\right)\right)}^2}{{\left(F\left(x\right)P\left(y\right)-1\right)}^2}\right)}{2\left(F\left(x\right)P\left(y\right)-1\right)}=0$

AIA 2p Abel type equations ("Abel Inverse Abel" 2-parameterized class, see Abel ODE papers in dsolve/references):

diff(y(x),x) = 2*y(x)*(2*x-3+b^2+y(x))/(y(x)*x+(x-1)*(-4+b^2+3*x));

dsolve((34));

$\mathrm{c\_\_1}+\frac{\left(b-1\right)\mathrm{LegendreP}\left(1\,b\,-\frac{-1+x}{\sqrt{-y}}\right)-\sqrt{-y}\mathrm{LegendreP}\left(0\,b\,-\frac{-1+x}{\sqrt{-y}}\right)}{\left(1-b\right)\mathrm{LegendreQ}\left(1\,b\,-\frac{-1+x}{\sqrt{-y}}\right)+\sqrt{-y}\mathrm{LegendreQ}\left(0\,b\,-\frac{-1+x}{\sqrt{-y}}\right)}=0$

Explicit solution free of integrals for Riccati type equations; in previous releases, the solution was implicit with integrals.

diff(y(x),x) = (x-exp(x))*y(x)^2+(exp(2*x)-x^2)/(-x+exp(x))*y(x)-exp(x)*(-x^2+x-1+x*exp(x))/(-x+exp(x))^2;

$\mathrm{y\text{'}}=\left(x-{\ⅇ}^x\right)y^2+\frac{\left({\ⅇ}^{2x}-x^2\right)y}{-x+{\ⅇ}^x}-\frac{{\ⅇ}^x\left(-x^2+x-1+x{\ⅇ}^x\right)}{{\left(-x+{\ⅇ}^x\right)}^2}$

dsolve((36),[sym_pat]); # the optional use of [sym_pat] methods speeds up this process

$y=\frac{\tan \left(\frac{\sqrt{4{\ⅇ}^x{x}^3-x^4-6{\ⅇ}^{2x}{x}^2+4{\ⅇ}^{3x}x-{\ⅇ}^{4x}}\left(-x^2+2{\ⅇ}^x-2\mathrm{c\_\_1}\right)}{4\left(x^2-2x{\ⅇ}^x+{\left({\ⅇ}^x\right)}^2\right)}\right)\sqrt{4{\ⅇ}^x{x}^3-x^4-6{\ⅇ}^{2x}{x}^2+4{\ⅇ}^{3x}x-{\ⅇ}^{4x}}+x^2-{\ⅇ}^{2x}}{2\left(-{\ⅇ}^{2x}+2x{\ⅇ}^x-x^2\right)}$

Explicit solution - free of RootOfs - with integrals of the arbitrary function $\mathrm{\sigma }\left(x\right)$; in previous releases, the solution was implicit or expressed using RootOf with integrals.

sigma(x) = -(y(x)-x*diff(y(x),x))*diff(y(x),x)/x/y(x)/diff(diff(y(x),x),x);

$\mathrm{\sigma }\left(x\right)=-\frac{\left(y-\mathrm{y\text{'}}x\right)\mathrm{y\text{'}}}{xy\mathrm{y\text{'}\text{'}}}$

dsolve((38));

$y=\mathrm{c\_\_2}{\ⅇ}^{\int \frac{{\ⅇ}^{-\left(\int \frac{1}{\mathrm{\sigma }\left(x\right)x}\ⅆx\right)}}{\int \frac{{\ⅇ}^{-\left(\int \frac{1}{\mathrm{\sigma }\left(x\right)x}\ⅆx\right)}\left(\mathrm{\sigma }\left(x\right)-1\right)}{\mathrm{\sigma }\left(x\right)}\ⅆx-\mathrm{c\_\_1}}\ⅆx}$

Simpler form of the solution; in previous releases, the solution was RootOf expressions with integrals.

diff(y(x),x,x)*(x-y(x))+2*diff(y(x),x)*(1+diff(y(x),x)) = 0;

$\mathrm{y\text{'}\text{'}}\left(-y+x\right)+2\mathrm{y\text{'}}\left(1+\mathrm{y\text{'}}\right)=0$

dsolve((40));

The dsolve function now uses LinearFunctionalSystems package routines for the computation of series solutions to linear ODE systems with polynomial coefficients.

PDEtools[declare](x(t), y(t), prime=t);

$x\left(t\right)\mathrm{will\; now\; be\; displayed\; as}x$

$y\left(t\right)\mathrm{will\; now\; be\; displayed\; as}y$

$\mathrm{derivatives\; with\; respect\; to}t\mathrm{of\; functions\; of\; one\; variable\; will\; now\; be\; displayed\; with\; \text{'}}$

{t*diff(x(t),t,t,t)-diff(y(t),t,t)-x(t)-t^3*y(t) = 0,
 t^2*diff(y(t),t,t)+diff(diff(x(t),t),t)-2*x(t)-t*y(t) = 0};

$\left\{t\mathrm{x\text{'}\text{'}\text{'}}-\mathrm{y\text{'}\text{'}}-x-t^{3}y=0\,t^2\mathrm{y\text{'}\text{'}}+\mathrm{x\text{'}\text{'}}-2x-ty=0\right\}$

dsolve((43), {x(t), y(t)}, 'series');

$\left\{x=x\left(0\right)+\mathrm{D}\left(x\right)\left(0\right)t+x\left(0\right)t^2+\left(\frac{\mathrm{D}\left(x\right)\left(0\right)}{3}+\frac{y\left(0\right)}{6}\right)t^3+\left(\frac{x\left(0\right)}{4}+\frac{\mathrm{D}\left(y\right)\left(0\right)}{12}\right)t^4+\left(-\frac{y\left(0\right)}{30}-\frac{x\left(0\right)}{40}-\frac{\mathrm{D}\left(x\right)\left(0\right)}{60}\right)t^5+O\left(t^6\right)\,y=y\left(0\right)+\mathrm{D}\left(y\right)\left(0\right)t-\frac{1}{2}x\left(0\right)t^2+\left(\frac{\mathrm{D}\left(x\right)\left(0\right)}{6}+\frac{y\left(0\right)}{6}\right)t^3+\left(\frac{5x\left(0\right)}{12}+\frac{\mathrm{D}\left(y\right)\left(0\right)}{6}\right)t^4+\left(-\frac{19y\left(0\right)}{120}-\frac{3x\left(0\right)}{40}-\frac{\mathrm{D}\left(x\right)\left(0\right)}{15}\right)t^5+O\left(t^6\right)\right\}$

Using the HINT option for PDE systems containing arbitrary functions:

{-diff(xi(x,y),y,y)+diff(xi(x,y),y)*h(y), diff(eta(x,y),x)*g(x)+diff(eta(x,y),x,x), diff(eta(x,y),y)*h(y)+eta(x,y)*diff(h(y),y)+2*diff(xi(x,y),y)*g(x)-2*diff(xi(x,y),x,y)+diff(eta(x,y),y,y), 2*diff(eta(x,y),x,y)+diff(xi(x,y),x)*g(x)-diff(xi(x,y),x,x)+2*diff(eta(x,y),x)*h(y)+xi(x,y)*diff(g(x),x)};

$\left\{{\mathrm{\eta }}_{x}g\left(x\right)+{\mathrm{\eta }}_{x,x}\,-{\mathrm{\xi }}_{y,y}+{\mathrm{\xi }}_{y}h\left(y\right)\,{\mathrm{\eta }}_{y}h\left(y\right)+\mathrm{\eta }\left(x\,y\right)h_y+2{\mathrm{\xi }}_{y}g\left(x\right)-2{\mathrm{\xi }}_{x,y}+{\mathrm{\eta }}_{y,y}\,2{\mathrm{\eta }}_{x,y}+{\mathrm{\xi }}_{x}g\left(x\right)-{\mathrm{\xi }}_{x,x}+2{\mathrm{\eta }}_{x}h\left(y\right)+\mathrm{\xi }\left(x\,y\right)g_x\right\}$

pdsolve( (45), [xi,eta], HINT=`+`);

$\left\{\mathrm{\eta }\left(x\,y\right)=\mathrm{c\_\_1}+\left(\mathrm{c\_\_6}+\int \left(\mathrm{c\_\_5}-\mathrm{c\_\_1}h\left(y\right)\right){\ⅇ}^{\int h\left(y\right)\ⅆy}\ⅆy\right){\ⅇ}^{-\left(\int h\left(y\right)\ⅆy\right)}\,\mathrm{\xi }\left(x\,y\right)=\left(\mathrm{c\_\_4}+\int \left(\mathrm{c\_\_3}+\mathrm{c\_\_2}g\left(x\right)\right){\ⅇ}^{\int -g\left(x\right)\ⅆx}\ⅆx\right){\ⅇ}^{-\left(\int -g\left(x\right)\ⅆx\right)}+\mathrm{c\_\_2}\right\}$

pdetest((46),(45));

$\left\{0\right\}$

Modifications to the taylorseries method to allow more efficient computation of high accuracy solutions.

Addition of interactive features to the procedure returns of the numeric ODE solvers allows you to query and change the initial data, query the last computed point (very useful for problems with singularities), and query the method used in the procedure, for example, rkf45. For more information, see dsolve/numeric/IVP.

Significant efficiency improvements to the numeric ODE boundary value problem solver, especially for difficult problems.

Ability to solve ODE systems containing user-defined functions. For more information, see the known option in dsolve/numeric.

Extension of all methods (schemes) to solve nonlinear PDEs that are linear in the highest time derivative.

Computation and plotting of error estimates for the solution as a function of time and space. For more information, see pdsolve/numeric/errorcontrol.

Addition of time-adaptive solution computation. For example, you can now compute a solution to within a specified error per unit time step. Note: Only the time step can be adapted. The solution must be computed with a number of spatial points sufficient to keep the spatial error within tolerance. This condition is monitored, and an error results if the spatial mesh is too coarse. For more information, see pdsolve/numeric/errorcontrol.