start = .25 ending = 2 steps = 10000 e = linspace(start,ending,steps); NOISE = 0; while (1+NOISE) < 1000; X = round(20.*rand(1,10)).*(.48+(rand*(NOISE/1000)*(-1)^(-round(rand)))); n1 = X(1)./e; n2 = X(2)./e; n3 = X(3)./e; n4 = X(4)./e; n5 = X(5)./e; n6 = X(6)./e; n7 = X(7)./e; n8 = X(8)./e; n9 = X(9)./e; n10 = X(10)./e; K = ((n1-round(n1)).^2+(n2-round(n2)).^2+(n3-round(n3)).^2+(n4-round(n4)).^2+(n5-round(n5)).^2+(n6-round(n6)).^2+(n7-round(n7)).^2+(n8-round(n8)).^2+(n9-round(n9)).^2+(n10-round(n10)).^2); [C,I]=min(K); V(NOISE+1) = .25 + I*((ending-start)/steps); NOISE = NOISE+1; end Nose = linspace(0,998,999); plot(Nose,V) hold y=.48.*ones(size(Nose)); plot(Nose,y,'red')