load
plot(the, elit,'*k')
hold on;
plot(the, eold,'or')
plot(the, enew,'ob')
xlim([0,90])
ylim([0, 2.0])