clear all
open('leaky.mdl')
mu_range = logspace(0.1,100,100)';
figure(1)
ymax = zeros(length(mu_range),1);

for i = 1:length(mu_range)
    mu = mu_range(i);
    [t,x,y] = sim('leaky.mdl');
    ymax(i) = max(max(y));
    plot(t,y);
    hold on;
end