delta=0.45;K=100;lpad=20;iter=400;rho = [0 0 0 0 0.24123 0.75877];lambda = [0 0.38354 0.04237 0.57409];%rho = [0 0 0 0 0 1];%lambda = [0 0 1];mu_start = 2.5;mu_inc = 0.5;mu_term = 0.01;ticy = thresh_gauss_irreg(rho,lambda,K,delta,lpad,iter,mu_start,mu_inc,mu_term);toc