acc05dim1fin.m:
One Dimensional Test for Distributed Parameter Numerics
tolu=5.0e-06; toly1=5.0e-08; toly2=5.0e-07; toly3=5.0e-04;
Note:  toly2 and toly3 not used.
tolcory=5.0e-06; tolcorxi=5.0e-06;
Ncor=5; L=10
M=26; K=1600; mx=27; nt=1601,mu0=27; nu0=1601
T=5.0e+00; RX(1)=1.0e+00; RT=2.0e-02; RC=1.0e-02;
dt=3.1e-03; dx(1)=7.7e-02; s(3)=2.0e-01;
D=[4.2e-03,1.0e-15,2.2e-01];
a=[1.2e-02,8.6e-07,1.1e+01];
[k1,k2]=[1.0e+00,1.0e+00];
[wy1,wy3,wu3]=[1.0e-03,1.0e-04,1.0e-05];
[alp12,kap13,alp21,kap23]=[1.0e-04,1.0e-04,1.0e-04,1.0e-04];
q=[1.0e-01,0.0e+00,1.0e-01];
r=[1.0e-01,0.0e+00,0.0e+00];
s=[0.0e+00,0.0e+00,2.0e-01];
EPMR: max(D)*dt/dx(1)^2=1.1e-01;
Y(1,x,1)= 0.0e+00; 0.0e+00; 0.0e+00; 0.0e+00; 0.0e+00; 
Y(2,x,1)= 1.0e+00; 1.0e+00; 1.0e+00; 1.0e+00; 1.0e+00; 
Y(3,x,1)= 0.0e+00; 0.0e+00; 0.0e+00; 0.0e+00; 0.0e+00; 
mu3=[1,27,1601]; mxi3=[1,27,1601];

Ncormaxy=2; Ncormaxxi=1; Lmax=2;
DU3norm=1.2e-06; tolu*U3oldnorm=3.2e-06; U3oldnorm=6.4e-01;
DY1norm=7.3e-14; toly1*Y1oldnorm=4.8e-07; Y1oldnorm=9.6e+00;

DY2norm=2.2e-11; toly2*Y2oldnorm=8.0e-04; Y2oldnorm=1.6e+03;
DY3norm=2.4e-06; toly3*Y3oldnorm=2.4e-05; Y3oldnorm=4.7e-02;

Y(1,x,1)= 0.0e+00; 0.0e+00; 0.0e+00; 4.6e-118; 2.0e-02; 7.9e-160; 0.0e+00; 
Y(1,x,2)= -4.7e-217; 2.3e-271; 1.4e-216; 6.7e-40; 2.0e-02; 2.2e-62; -4.7e-217; 
Y(1,x,3)= -2.4e-100; 2.3e-135; 7.3e-100; 3.1e-19; 2.0e-02; 9.1e-23; -2.4e-100; 
Y(1,x,4)= -1.9e-41; 3.3e-57; 5.7e-41; 4.0e-18; 2.0e-02; 3.5e-21; -1.9e-41; 
Y(1,x,401)= -6.5e-12; 1.3e-12; 2.5e-11; 2.6e-06; 6.7e-03; 3.3e-07; -6.5e-12; 
Y(1,x,801)= -1.9e-09; 1.2e-09; 1.1e-08; 4.0e-05; 4.5e-03; 9.6e-06; -1.9e-09; 
Y(1,x,1201)= -1.9e-08; 3.9e-08; 2.1e-07; 1.3e-04; 3.7e-03; 4.5e-05; -1.9e-08; 
Y(1,x,1601)= 5.7e-09; 3.4e-07; 1.3e-06; 2.5e-04; 3.2e-03; 1.1e-04; 5.7e-09; 
100*y(1,mxm,nt)/y(1,mxm,1) =  16.21%;