The functional partial differential equation (FPDE) for cell division, $$ \begin{align*} &\frac{\partial}{\partial t}n(x,t) +\frac{\partial}{\partial x}(g(x,t)n(x,t))\\ &\quad = -(b(x,t)+\mu(x,t))n(x,t)+b(\alpha x,t)\alpha n(\alpha x,t)+b(\beta x,t)\beta n(\beta x,t), \end{align*} $$![](//static-cambridge-org.ezproxyberklee.flo.org/binary/version/id/urn:cambridge.org:id:binary:20210420022557267-0751:S1446181121000055:S1446181121000055_eqnu1.png)
is not amenable to analytical solution techniques, despite being closely related to the first-order partial differential equation (PDE) $$ \begin{align*} \frac{\partial}{\partial t}n(x,t) +\frac{\partial}{\partial x}(g(x,t)n(x,t)) = -(b(x,t)+\mu(x,t))n(x,t)+F(x,t), \end{align*} $$![](//static-cambridge-org.ezproxyberklee.flo.org/binary/version/id/urn:cambridge.org:id:binary:20210420022557267-0751:S1446181121000055:S1446181121000055_eqnu2.png)
which, with known
$F(x,t)$, can be solved by the method of characteristics. The difficulty is due to the advanced functional terms
$n(\alpha x,t)$ and
$n(\beta x,t)$, where
$\beta \ge 2 \ge \alpha \ge 1$, which arise because cells of size x are created when cells of size
$\alpha x$ and
$\beta x$ divide.
The nonnegative function,
$n(x,t)$, denotes the density of cells at time t with respect to cell size x. The functions
$g(x,t)$,
$b(x,t)$ and
$\mu (x,t)$ are, respectively, the growth rate, splitting rate and death rate of cells of size x. The total number of cells,
$\int _{0}^{\infty }n(x,t)\,dx$, coincides with the
$L^1$ norm of n. The goal of this paper is to find estimates in
$L^1$ (and, with some restrictions,
$L^p$ for
$p>1$) for a sequence of approximate solutions to the FPDE that are generated by solving the first-order PDE. Our goal is to provide a framework for the analysis and computation of such FPDEs, and we give examples of such computations at the end of the paper.