d  (* data: h, t, c(0,t), s(t) *)

At   t = 0.01

j=3;  (*  c  *)
a={{d[[1,1]],d[[1,j]],"", "","",""}};
For[i=1,i<5,i++,
h=d[[i+1,1]];
ih=d[[i+1,j]];
i2h=d[[i,j]];
er=(ih-i2h)/3;
r=ih+er;
If[i>1,er2=(r-rold)/15; r2=r+er2;
a=Append[a,{h,ih,r,er/h^2,r2,er2/h^4}],
a=Append[a,{h,ih,r,er/h^2,"",""}]];
rold=r];
TableForm[a]
ermin=Min[a[[3,4]],a[[4,4]],a[[5,4]]];
ermax=Max[a[[3,4]],a[[4,4]],a[[5,4]]];
er=(ermax-ermin)d[[5,1]]^2;
Print["Rc from ",a[[5,3]]-er, "  to  ", a[[5,3]]+er,"  error = ",er]
ermin=Min[a[[3,6]],a[[4,6]],a[[5,6]]];
ermax=Max[a[[3,6]],a[[4,6]],a[[5,6]]];
er=(ermax-ermin)d[[5,1]]^4;
Print["R2c from ",a[[5,5]]-er, "  to  ", a[[5,5]]+er,"  error = ",er]

j=4;   (*  s  *)
a={{d[[1,1]],d[[1,j]],"", "","",""}};
For[i=1,i<5,i++,
h=d[[i+1,1]];
ih=d[[i+1,j]];
i2h=d[[i,j]];
er=(ih-i2h)/3;
r=ih+er;
If[i>1,er2=(r-rold)/15; r2=r+er2;
a=Append[a,{h,ih,r,er/h^2,r2,er2/h^4}],
a=Append[a,{h,ih,r,er/h^2,"",""}]];
rold=r];
TableForm[a]
ermin=Min[a[[3,4]],a[[4,4]],a[[5,4]]];
ermax=Max[a[[3,4]],a[[4,4]],a[[5,4]]];
er=(ermax-ermin)d[[5,1]]^2;
Print["Rs from ",a[[5,3]]-er, "  to  ", a[[5,3]]+er,"  error = ",er]
ermin=Min[a[[3,6]],a[[4,6]],a[[5,6]]];
ermax=Max[a[[3,6]],a[[4,6]],a[[5,6]]];
er=(ermax-ermin)d[[5,1]]^4;
Print["R2s from ",a[[5,5]]-er, "  to  ", a[[5,5]]+er,"  error = ",er]