clear all;
close all;
relamdt=-4:0.1:4;
imlamdt=-4:0.1:4;
[RElamdt,IMlamdt]=meshgrid(relamdt,imlamdt);
axis square;
lamdt=RElamdt+i*IMlamdt;
sig=1+lamdt+lamdt.^2/2;
contour_levels=[1 2 3 4];
contour(RElamdt,IMlamdt,abs(sig),contour_levels,’linewidth’,4);
xlabel(‘Re \lambda \Delta t’,’Fontsize’,14);
ylabel(‘Im \lambda \Delta t’,’Fontsize’,14);
hold on;
plot([-4 4],[0 0],’k’);
plot([0 0],[-4 4],’k’);
title(‘contour |\sigma|’,’Fontsize’,14);