% Sitzung S. 216/217

u = inline('heat([0,0],[1,1],t,0,0,[5,0,0,0],h)','t','h');
t = 0.424; order = 2; tol = 5e-14;
nmin = 4; [val,err,ampl] = richardson(tol,order,nmin,u,t);
val, err = max(err,ampl*eps)

nmin = 8; [val,err,ampl] = richardson(tol,order,nmin,u,t);
val, err = max(err,ampl*eps)
