c=15
P=15
r_d=15
E=18*10^6
t=0.08
nu= 0.3
K= E*t^3/(12-12*nu^2)
[r theta]=meshgrid(0:0.5:15, 0:2*pi/30:2*pi)
w= ((P*r_d^4)/64*K)*(1-(r/r_d )^2 )^2+c
[x,y]=pol2cart(theta,r)
plot3(x,y,w)
end
r=0:0.2:5Q=linspace(0,2*pi,50)
[r,Q]=meshgrid(r,Q)
z=r.*exp(-r.^2).*cos(Q)
x=r.*cos(Q)
y=r.*sin(Q)
mesh(x,y,z)