The piecewise function has a straightforward syntax. Every piece is specified by a Boolean condition followed by an expression. The next several Maple command lines make use of the following piecewise function:

This can be evaluated at arbitrary points. For example,

Because of the division by zero,  points such as x = 1 cannot be substituted. However, $\underset{x\to 1}{lim}f$  can be determined.

Such functions can be plotted to determine their behavior.

Besides evaluating limits, you can do operations such as computing derivatives, integrating, and solving differential equations with piecewise functions. Examples of solving DEs will be illustrated later. Using the same function, f(x), find its piecewise derivative.

Integrating this function gives:

Evaluate the integral with values for the lower and upper bounds.

where Si(x) is the Sine integral function.

Note: This works because discont is able to determine the potential discontinuities of piecewise functions. For example,

Consider a new function, g(x).

The piecewise function allows for common manipulations, such as simplifications.  The addition of the selector 'piecewise' indicates to simplify that it should only do simplifications as they apply to piecewise functions.  This is more efficient, in general.

Piecewise functions can be converted to a number of other functions and distributions.

Distributions can be converted back to piecewise functions.

The function piecewise lets us work with the CurveFitting[Spline] command. For example,

This spline can be graphed.

But now we can also integrate it.

Other piecewise functions can also be converted to piecewise and be properly manipulated. For example,

Another example of conversion, in this case with signum and abs.

Note: The discontinuity at 0 in the above function.

Another simple example:

It turns out to be a well-behaved, non-piecewise function.

We can also determine the differentiability class of a piecewise continuous function. For example, in the case of our previous spline function, newcubic, we have

This must be true for splines! However, when we check to see if it is $C^3$, we obtain

To determine the highest order of continuity and the problem points, enter:

Define f(x) as:

Derivatives can be found, and piecewise functions are returned.

The graph of the derivative can be drawn as well.

The straightforward calculation of the series of f around x=0 can be computed by using the series command.

This produces an answer with a superfluous order term. However, series can do better when using piecewise. Thus, it can be converted to a new function g(x).

We can then also plot this max function.

LaTeX output can also be generated.

Similarly, everything also works with min.

Converting to piecewise also works with _Envsignum0 := 1;

We can solve differential equations with piecewise functions in the coefficients. The type of equations that one can solve include all first-order methods using integration, Riccati, and higher-order methods including linear, Bernoulli, and Euler. Some examples are included in the sections that follow.

First-order Linear with Piecewise (Nontrivial) Coefficients

Consider the following DE: >  $\mathrm{eq}:=\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)\&plus;\mathrm{piecewise}\left(x<x^2-4\&comma;{\&ExponentialE;}^x\right)y\left(x\right)$ $\mathrm{eq}≔\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)+\left(\left\{\begin{array}{cc}{\&ExponentialE;}^x& x<x^2-4\\ 0& \mathrm{otherwise}\end{array}\right.\right)y\left(x\right)$ (3.1.1)   Now solve this DE. >  $\mathrm{dsolve}\left(\mathrm{eq}\&comma;y\left(x\right)\right)$ $y\left(x\right)=\left\{\begin{array}{cc}{\&ExponentialE;}^{-{\&ExponentialE;}^x}\mathrm{c\_\_1}& x<\frac{1}{2}-\frac{\sqrt{17}}{2}\\ {\&ExponentialE;}^{-{\&ExponentialE;}^{\frac{1}{2}-\frac{\sqrt{17}}{2}}}\mathrm{c\_\_1}& x<\frac{1}{2}+\frac{\sqrt{17}}{2}\\ {\&ExponentialE;}^{-{\&ExponentialE;}^{\frac{1}{2}-\frac{\sqrt{17}}{2}}-{\&ExponentialE;}^x+{\&ExponentialE;}^{\frac{1}{2}+\frac{\sqrt{17}}{2}}}\mathrm{c\_\_1}& \frac{1}{2}+\frac{\sqrt{17}}{2}\le x\end{array}\right.$ (3.1.2) Now an example of a linear equation with polynomial coefficients and 0-1 perturbation. >  $\mathrm{eq}:=\frac{{\&DifferentialD;}^2}{\&DifferentialD;x^2}y\left(x\right)\&plus;2x\left(\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)\right)\&plus;2y\left(x\right)\&equals;\mathrm{piecewise}\left(0<x\&comma;1\right)$ $\mathrm{eq}≔\frac{{\&DifferentialD;}^2}{\&DifferentialD;x^2}y\left(x\right)+2x\left(\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)\right)+2y\left(x\right)=\left\{\begin{array}{cc}1& 0<x\\ 0& \mathrm{otherwise}\end{array}\right.$ (3.1.3)   Solving the equation gives: >  $\mathrm{dsolve}\left(\mathrm{eq}\&comma;y\left(x\right)\right)$ $y\left(x\right)=\left\{\begin{array}{cc}{\&ExponentialE;}^{-x^2}\mathrm{erfi}\left(x\right)\mathrm{c\_\_1}+{\&ExponentialE;}^{-x^2}\mathrm{c\_\_2}& x<0\\ {\&ExponentialE;}^{-x^2}\mathrm{erfi}\left(x\right)\mathrm{c\_\_1}+{\&ExponentialE;}^{-x^2}\mathrm{c\_\_2}+\frac{{\&ExponentialE;}^{-x^2}{\&ExponentialE;}^{x^2}}{2}-\frac{{\&ExponentialE;}^{-x^2}}{2}& 0\le x\end{array}\right.$ (3.1.4) Now we attempt to solve a DE with y(x) inside the piecewise function. >  $\mathrm{eq}:=\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)\&equals;\mathrm{piecewise}\left(1<x\&comma;y\left(x\right)\&comma;1\right)$ $\mathrm{eq}≔\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)=\left\{\begin{array}{cc}y\left(x\right)& 1<x\\ 1& \mathrm{otherwise}\end{array}\right.$ (3.1.5) >  $\mathrm{dsolve}\left(\mathrm{eq}\&comma;y\left(x\right)\right)$ $y\left(x\right)=\left\{\begin{array}{cc}\mathrm{c\_\_1}+x& x<1\\ {\&ExponentialE;}^{x-1}\mathrm{c\_\_1}+{\&ExponentialE;}^{x-1}& 1\le x\end{array}\right.$ (3.1.6) And now solve a nonlinear differential equation with discontinuities depending on parameters. The assumption is needed because this function is discontinuous in 0 and a. >  $\mathrm{assume}\left(a\&gt;0\right)\&colon;$ >  $g:=\mathrm{convert}\left(-\mathrm{signum}\left(x\right) \cdot \left|a-\left|x\right|\right|\&comma;\&apos;\mathrm{piecewise}\&apos;\&comma;x\right)$ $g≔\left\{\begin{array}{cc}-\mathrm{a~}-x& x\le -\mathrm{a~}\\ \mathrm{a~}+x& x<0\\ -\mathrm{a~}+x& x<\mathrm{a~}\\ \mathrm{a~}-x& \mathrm{a~}\le x\end{array}\right.$ (3.1.7) >  $\mathrm{dsolve}\left(\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)\&equals;g \cdot y{\left(x\right)}^2\&comma;y\left(x\right)\right)$ $y\left(x\right)=\left\{\begin{array}{cc}\frac{2}{2x\mathrm{a~}+x^2+2\mathrm{c\_\_1}}& x<-\mathrm{a~}\\ \frac{2}{-2{\mathrm{a~}}^2-2x\mathrm{a~}-x^2+2\mathrm{c\_\_1}}& x<0\\ \frac{2}{-2{\mathrm{a~}}^2+2x\mathrm{a~}-x^2+2\mathrm{c\_\_1}}& x<\mathrm{a~}\\ \frac{2}{-2x\mathrm{a~}+x^2+2\mathrm{c\_\_1}}& \mathrm{a~}\le x\end{array}\right.$ (3.1.8)

Linear DE with Constant Coefficients and Piecewise Forcing Function

>  $\mathrm{eq}:=\frac{{\&DifferentialD;}^2}{\&DifferentialD;x^2}y\left(x\right)\&plus;y\left(x\right)\&equals;\mathrm{piecewise}\left(x<0\&comma;-x\&comma;0<x\&comma;x\right)$ $\mathrm{eq}≔\frac{{\&DifferentialD;}^2}{\&DifferentialD;x^2}y\left(x\right)+y\left(x\right)=\left\{\begin{array}{cc}-x& x<0\\ x& 0<x\end{array}\right.$ (3.2.1) The solution with given initial values: >  $\mathrm{dsolve}\left(\left\{y\left(0\right)\&equals;0\&comma;\mathrm{D}\left(y\right)\left(0\right)\&equals;1\&comma;\mathrm{eq}\right\}\&comma;y\left(x\right)\right)$ $y\left(x\right)=\left\{\begin{array}{cc}2\sin \left(x\right)-x& x<0\\ x& 0\le x\end{array}\right.$ (3.2.2) The general solution: >  $\mathrm{dsolve}\left(\mathrm{eq}\&comma;y\left(x\right)\right)$ $y\left(x\right)=\left\{\begin{array}{cc}\sin \left(x\right)\mathrm{c\_\_2}+\cos \left(x\right)\mathrm{c\_\_1}-x& x<0\\ \sin \left(x\right)\mathrm{c\_\_2}+\cos \left(x\right)\mathrm{c\_\_1}+x-2\sin \left(x\right)& 0\le x\end{array}\right.$ (3.2.3)   This also works in the case of output = basis. >  $\mathrm{dsolve}\left(\mathrm{eq}\&comma;y\left(x\right)\&comma;\mathrm{output}\&equals;\mathrm{basis}\right)$ $\left[\left[\cos \left(x\right)\&comma;\sin \left(x\right)\right]\&comma;-\left(\left\{\begin{array}{cc}x& x\le 0\\ 2\sin \left(x\right)-x& 0<x\end{array}\right.\right)\right]$ (3.2.4)

Nonlinear Riccati DE with Piecewise Continuous Coefficients

Consider the following differential equation. >  $\mathrm{eq}:=\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)\&equals;\mathrm{piecewise}\left(x\&gt;0\&comma;x\right)\cdot y{\left(x\right)}^2$ $\mathrm{eq}≔\frac{\&DifferentialD;}{\&DifferentialD;x}y\left(x\right)=\left(\left\{\begin{array}{cc}x& 0<x\\ 0& \mathrm{otherwise}\end{array}\right.\right){y\left(x\right)}^2$ (3.3.1) >  $\mathrm{eq\_sol}:=\mathrm{dsolve}\left(\mathrm{eq}\&comma;y\left(x\right)\right)$ $\mathrm{eq\_sol}≔y\left(x\right)=\left\{\begin{array}{cc}\frac{1}{\mathrm{c\_\_1}}& x<0\\ -\frac{2}{x^2-2\mathrm{c\_\_1}}& 0\le x\end{array}\right.$ (3.3.2)