clear all; 



tauL = 2.5e-6;

tauP = 0.1e-6;

r = linspace(1,200,1e6);

dr = -r/(2*tauL); % Damping rate

wr = sqrt((r-1)/(tauL*tauP)-(r/(2*tauL)).^2); % Relaxation freq [s^-1]

sp = dr+j*wr;

sm = dr-j*wr;



plot(sp); hold on;

plot(sm,'r');

grid on;

xlabel('Re(s_{\pm})');

ylabel('Im(s_{\pm})');

title('6.638: Problem 5.1c, Root locus');