% example for pseudospectral computations

N=6
A=toeplitz([0 1/4 zeros(1,N-2)],[0 1 zeros(1,N-2)])
D=diag(2.^(1:N))
S=D*A*D^-1
N=64;
AA=toeplitz([0 1/4 zeros(1,N-2)],[0 1 zeros(1,N-2)]);
DD=diag(2.^(1:N));
SS=DD*AA*DD^-1;