28 May 2024 02:32:35 PM

demo2():
  C version
  Solve heat equation with exact solution
  u(x,t) = erf ( 0.5 * x / sqrt ( t ) ) )

  Final time T = 0.5
  Number of spatial nodes n = 10
  Number of time steps m = 32

  Explicit finite difference scheme:
  r = 0.390625
  Geomview script written to 'demo2_explicit.gv'
  Maple script written to 'demo2_explicit.mpl'
  Matlab script written to 'demo2_explicit.m'
  Absolute error = 0.132726

  Implicit finite difference scheme:
  r = 0.390625
  Geomview script written to 'demo2_implicit.gv'
  Maple script written to 'demo2_implicit.mpl'
  Matlab script written to 'demo2_implicit.m'
  Absolute error = 0.0484173

  Crank Nicolson linear equation solver:
  r = 0.390625
  Geomview script written to 'demo2_crank_nicolson.gv'
  Maple script written to 'demo2_crank_nicolson.mpl'
  Matlab script written to 'demo2_crank_nicolson.m'
  Absolute error = 0.0287708

  Seidman sweep scheme:
  r = 0.195312
  Geomview script written to 'demo2_seidman_sweep.gv'
  Maple script written to 'demo2_seidman_sweep.mpl'
  Matlab script written to 'demo2_seidman_sweep.m'
  Absolute error = 0.0210214

demo2():
  Normal end of execution.
28 May 2024 02:32:35 PM