The PDEplot command produces the plot of the solution for a first order, linear or nonlinear partial differential equation (PDE), for given initial conditions.

Given initial data consisting of equations of a manifold in $\left(n\+1\right)$-space depending on $n-1$ parameters (for $n=2$, this is a curve) through which the required integral hypersurface must pass, the method of characteristics is used to generate the solution, resulting in a system of coupled first order ODEs to which numerical methods can be applied. The default is to use an internal Runge-Kutta fourth order (fixed step size) method. Other dsolve/numeric methods can be chosen by way of the method optional equation. Note, however, that use of these external methods tends to slow considerably the creation of plots.

The first argument, PDE, can be any first order partial differential equation. PDEs involving piecewise functions are handled using hardware arithmetic (evalhf) as of Maple V Release 4. All arguments after the first can be given in any order.

The inits parameter is a list of $n+1$ elements specifying the parametric form of a curve in $\left(n\+1\right)$-space through which the required integral surface must pass. The elements should depend on $n-1$ parameters. Using s and t as the parameters with $n=3$, for example, the elements of the list must be expressions such as $\left[\mathrm{f1}\left(s\,t\right)\,\mathrm{f2}\left(s\,t\right)\,\mathrm{f3}\left(s\,t\right)\,\mathrm{f4}\left(s\,t\right)\right]$, or equations of the form $\left[\mathrm{x2}=\mathrm{f1}\left(s\,t\right)\,u=\mathrm{f2}\left(s\,t\right)\,\mathrm{x1}=\mathrm{f3}\left(s\,t\right)\,\mathrm{x3}=\mathrm{f4}\left(s\,t\right)\right]$, where the $f\left(s\,t\right)$ are explicit expressions of s and t that can be evaluated numerically. The two possible inputs are equivalent but if equations are given the order of the elements is irrelevant.

The srange parameter is a list of ranges for the parameters used in the initial data. It should be given as $\[\mathrm{s1}=\mathrm{a1}\mathrm{..}\mathrm{b1},\mathrm{s2}=\mathrm{a2}\mathrm{..}\mathrm{b2},...\]$.

If $n=2$, by default the solution surface is displayed as an opaque, colored membrane, with the initial conditions highlighted in black. For $2<n$, by default a sequence of surfaces is shown in animation which together map out the required solution hypersurface. For more information, see plots[animate].

Remaining arguments are interpreted as options which are specified as equations of the form option = value. In particular, the style option allows one to plot the points as points only, or to interpolate them using line mode. See plot3d,option for more information.

This function is part of the PDEtools package, and so it can be used in the form PDEplot(..) only after executing the command with(PDEtools). However, it can always be accessed through the long form of the command by using PDEtools[PDEplot](..).

$\mathrm{with}\left(\mathrm{PDEtools}\right)\&colon;$

$\mathrm{pde1}≔\mathrm{diff}\left(u\left(x\&comma;y\right)\&comma;x\right)\mathrm{diff}\left(u\left(x\&comma;y\right)\&comma;y\right)-xy+u\left(x\&comma;y\right)=0$

$\mathrm{pde1}≔\left(\frac{\partial }{\partial x}u\left(x\&comma;y\right)\right)\left(\frac{\partial }{\partial y}u\left(x\&comma;y\right)\right)-xy+u\left(x\&comma;y\right)=0$

$\mathrm{PDEplot}\left(\mathrm{pde1}\&comma;\left[\cos \left(t\right)\&comma;\sin \left(t\right)\&comma;0\right]\&comma;t=-\mathrm{\pi }..\mathrm{\pi }\&comma;\mathrm{ic\_assumptions}=\left[\mathrm{diff}\left(u\left(x\&comma;y\right)\&comma;x\right)=-\cos \left(t\right)\right]\right)$

$\mathrm{PDEplot}\left(\mathrm{pde1}\&comma;\left[\cos \left(t\right)\&comma;\sin \left(t\right)\&comma;0\right]\&comma;t=-\mathrm{\pi }..\mathrm{\pi }\&comma;\mathrm{ic\_assumptions}=\left[\mathrm{diff}\left(u\left(x\&comma;y\right)\&comma;x\right)=\cos \left(t\right)\right]\right)$

$\mathrm{pde2}≔\left(y^2+{z\left(x\&comma;y\right)}^2+x^2\right)\mathrm{diff}\left(z\left(x\&comma;y\right)\&comma;x\right)-2xy\mathrm{diff}\left(z\left(x\&comma;y\right)\&comma;y\right)-2z\left(x\&comma;y\right)x=0$

$\mathrm{pde2}≔\left(y^2+{z\left(x\&comma;y\right)}^2+x^2\right)\left(\frac{\partial }{\partial x}z\left(x\&comma;y\right)\right)-2xy\left(\frac{\partial }{\partial y}z\left(x\&comma;y\right)\right)-2z\left(x\&comma;y\right)x=0$

$\mathrm{PDEplot}\left(\mathrm{pde2}\&comma;z\left(x\&comma;y\right)\&comma;\left[t\&comma;t\&comma;\frac{\sin \left(\frac{\mathrm{\pi }t}{0.1}\right)}{10}\right]\&comma;t=0..0.1\&comma;\mathrm{numchar}=40\&comma;\mathrm{orientation}=\left[-163\&comma;56\right]\&comma;\mathrm{basechar}=\mathrm{true}\&comma;\mathrm{numsteps}=\left[20\&comma;20\right]\&comma;\mathrm{stepsize}=0.15\&comma;\mathrm{initcolor}=\cos \left(t\right)t\&comma;\mathrm{style}=\mathrm{surfacecontour}\right)$

$\mathrm{pde3}≔\sin \left(\mathrm{diff}\left(u\left(x\&comma;y\&comma;z\right)\&comma;y\right)\right)=u\left(x\&comma;y\&comma;z\right)$

$\mathrm{pde3}≔\sin \left(\frac{\partial }{\partial y}u\left(x\&comma;y\&comma;z\right)\right)=u\left(x\&comma;y\&comma;z\right)$

$\mathrm{ics}≔\left[\cos \left(t\right)\sin \left(s\right)\&comma;\cos \left(s\right)\cos \left(t\right)\&comma;\cos \left(t\right)\&comma;\sin \left(t\right)\right],\left[t=0..\mathrm{\pi }\&comma;s=0..\mathrm{\pi }\right]$

$\mathrm{ics}≔\left[\cos \left(t\right)\sin \left(s\right)\&comma;\cos \left(s\right)\cos \left(t\right)\&comma;\cos \left(t\right)\&comma;\sin \left(t\right)\right],\left[t=0..\mathrm{\pi }\&comma;s=0..\mathrm{\pi }\right]$

$\mathrm{PDEplot}\left(\mathrm{pde3}\&comma;\mathrm{ics}\&comma;\mathrm{numsteps}=\left[-5\&comma;6\right]\&comma;\mathrm{stepsize}=0.1\&comma;\mathrm{axes}=\mathrm{none}\&comma;\mathrm{style}=\mathrm{surface}\&comma;\mathrm{numchar}=\left[16\&comma;16\right]\&comma;\mathrm{orientation}=\left[225\&comma;66\right]\&comma;\mathrm{lightmodel}='\mathrm{light2}'\right)$