% counter=1;
% for i=-1:0.01:1,
%   image(counter) = sqrt ( 0.20 - i.^2);
%   counter = counter + 1;
% end

clear;
e = rand;

width = 50 / 3;

for x=1:width,
y(x) = e * sqrt(width.^2 - x.^2);
end

plot(y)