r = linspace(0,100, 101)';
theta = linspace(0, 360, 61);
t = 1; % single t point for simplicity here
u = 1./r.*cos(2*pi*.05*r) .* ones(size(theta)); % generate some data to plot
xx = r.*cosd(theta);
yy = r.*sind(theta);
p = pcolor(xx, yy, u);
p.EdgeColor = 'none';