
function [e,r]=testtrapzPC(f,x,t0,T,p=1,N0=20,M=5)
%EXAMPLE
%[e,r]=testtrapzPC(f=@(x,t) x,x=@(t) exp(t),t0=0,T=1,p=1,N0=20,M=5)
  N=N0;
  e=zeros(M,1);
  for k=1:M,
   [XX]=trapzPC(f,t0,T,x(t0),N,p);
   e(k)=norm(XX(end,:)'-x(T),'inf');
   N*=2;
  end 
  r=e(1:M-1)./e(2:M);
endfunction
