R=30;
[x,y]=meshgrid(-100:100);
z=128*double(x.^2+y.^2<R^2);
z=z.*y/R;
z=z+128;

imshow(z/256)