function Forum()

N = 5; h = 2*pi/N; x = h*(1:N);
 
D1 = toeplitz([-pi^2/(3*h^2)-1/6 ...
                 -.5*(-1).^(1:N-1)./sin(h*(1:N-1)/2).^2])

D2 = toeplitz([-pi^2/(3*h^2)-1/6 -.5*(-1).^(1:N-1)./sin(h*(1:N-1)/2).^2])

d=D1-D2