c=[0,0];
t=linspace(0,2*pi,1000);
xx=[];
yy=[];
for radius=0:5:45
a=ones(1,1000)*radius;
[x,y] = pol2cart(t,a);
xx=[xx;x];
yy=[yy;y];
end

 plot(xx(2,:),yy(2,:))

v=10:10:100;
v=v';
v=v*ones(1,1000);
pcolor(xx,yy,v)

shading interp
axis square