cmna (version 1.0.5)

wave: Wave Equation using

Description

solve heat equation via forward-time central-space method

Usage

wave(u, alpha, xdelta, tdelta, n)

Arguments

u

the initial values of u

alpha

the thermal diffusivity coefficient

xdelta

the change in x at each step in u

tdelta

the time step

n

the number of steps to take

Value

a matrix of u values at each time step

Details

The heat solves the heat equation using the forward-time central-space method in one-dimension.

Examples

# NOT RUN {
speed <- 2
x0 <- 0
xdelta <- .05
x <- seq(x0, 1, xdelta)
m <- length(x)
u <- sin(x * pi * 2)
u[11:21] <- 0
tdelta <- .02
n <- 40
z <- wave(u, speed, xdelta, tdelta, n)
# }