% Shows use of simple sparse matrices
% You get print out after lines not ended with ;
% Everything after % is a comment
G=numgrid('S',5)
A=delsq(G);
n=size(A,1)
b=zeros(n,1);
b(7)=1;
x=A\b
X=v2g(x,G)
B=v2g(b,G)
surf(X)
figure(2)
mesh(B)