R/get.heat.R

get.heat <-
function(x, T0, k, t)
  {
	tim = t
        t1 = erf(x/(2*sqrt(k*tim)))
        T = T0*(0.5+0.5*c(t1))
    return(T)
  }

Try the geophys package in your browser

Any scripts or data that you put into this service are public.

geophys documentation built on May 1, 2019, 9:26 p.m.