Create a circular 3D surf plot

i have the following code:

%% Assign constant
v=0.33;% poisson ratio
E=73.1*10^9; % young's modulus [N]
ti=0.1185;% thickness of bulkhead intact [m]
td=0.1185;% thickness of bulkhead Damage[m]
d=3.11; % diameter of bulkhead[m] 
ai=d/2; % radius of bulkhead[m]
r= linspace(-ai,ai) ;
z=linspace(0,ti);
q0=5; % pressure distribution [N/m^2]
%% compute stress for Bulkhead intact state
Mr=0.0625*q0*(3+v)*(ai^2-r.^2);
Sigmari = (Mr.*z)/(ti^3/12);

the output " Sigmari" reflects the stress at some point on a circular plate as a function of thr radius r and the thickness z. I want to create a plot that would the plate and the shear distribution on the plate with a color legend.

Thank you in advance.

Answer 1

[Walter Roberson]

%% Assign constant
v=0.33;% poisson ratio
E=73.1*10^9; % young's modulus [N]
ti=0.1185;% thickness of bulkhead intact [m]
td=0.1185;% thickness of bulkhead Damage[m]
d=3.11; % diameter of bulkhead[m] 
ai=d/2; % radius of bulkhead[m]
r = linspace(-ai,ai).' ;
z=linspace(0,ti);
q0=5; % pressure distribution [N/m^2]
%% compute stress for Bulkhead intact state
Mr=0.0625*q0*(3+v)*(ai^2-r.^2);
Sigmari = (Mr.*z)/(ti^3/12);

[x, y] = pol2cart(r, z);
pcolor(x, y, Sigmari);
axis equal