You can find the convolution of any piecewise defined function
with another (off the domain of definition, they are assumed
to be zero). Here is
,
, and
, where
,
:
sage: x = PolynomialRing(QQ, 'x').gen() sage: f = Piecewise([[(0,1),1*x^0]]) sage: g = f.convolution(f) sage: h = f.convolution(g) sage: P = f.plot(); Q = g.plot(rgbcolor=(1,1,0)); R = h.plot(rgbcolor=(0,1,1))
show(P+Q+R)
.
See About this document... for information on suggesting changes.