% m2 --- Plot rootfcn(lambda) to find location of root
n = 10;
ll = linspace(0,13,n);
z = zeros(1,n);
for i = 1:n
  z(i) = rootfcn(ll(i));
end
plot(ll,z)
