M = zeros(402,402)

M(202:402,1:201)=1;

g = sqrt(40/(200*pi));

M(1,202) = 1;
for i=1:200
    M(1,202) = M(1,202) - g^2;
end;

for i=1:200
    M(1+i,202) = g*10*i;
end;

for i=1:200
    M(1,202+i) = -(10*i);
    for j=1:200
        M(1,202+i) = M(1,202+i) - g*10*j;
    end;
end;

for i=1:200
    M(i+1,202+i) = -((10*i)^2);
end;

eigs(M)