clc
clear
base=1000;
gamma=0.00001;
power=1;
lamda_min=5;
iter=1:20000:10000000;
lamda=max(base*(gamma*iter+1).^-1,lamda_min);
plot(iter,lamda)
clc
clear
base=1000;
gamma=0.00001;
power=1;
lamda_min=5;
iter=1:20000:10000000;
lamda=max(base*(gamma*iter+1).^-1,lamda_min);
plot(iter,lamda)