Heat Kernel for a One-Dimensional Rod
Calculate values of the heat kernel on a one-dimensional rod. The ends of the rod may be assumed to be insulated, or absorbing.
hotrod(len, xsource, xquery, sigma, ends=c("insulated", "absorbing"), nmax=20)
len |
Length of the rod. A single number or numeric vector. |
xsource |
Positions of the source points, from the left end of the rod
(in the same distance units as |
xquery |
Positions of the query points, from the left end of the rod
(in the same distance units as |
sigma |
Bandwidth for kernel. A single number or a numeric vector. |
ends |
Character string (partially matched) specifying whether the ends of the rod are assumed to be insulated or absorbing. |
nmax |
Number of terms in the infinite sum to use. A single integer or an integer vector. |
Computes the heat kernel as an infinite sum.
Number or numeric vector.
Greg McSwiggan and Adrian Baddeley Adrian.Baddeley@curtin.edu.au.
curve(hotrod(1, 0.1, x, 0.7)) # check it's a probability density f <- function(x) hotrod(1, 0.1, x, 0.7) integrate(f, 0, 1) ## absorbing ends curve(hotrod(1, 0.1, x, 0.7, ends="a"))
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.