رسم مقادیر یک ماتریس دو بعدی به صورت سیمی (mesh) و شامل دیواره های کناری، با دستور meshz ، در متلب (MATLAB)

سلام دوستان عزیز.من هر چی این کد متلب رو اجرا میکنم در خط  mesh(xx,yy,w)خطا میده.کسی متونه تصحیح بکنه؟clc
clear
close all
syms x y n m
format short g
P=5000; % Concentrated Load -->N
t=5; % Plate Thickness --> mm
E=207000; % Modulus of Elasticity --> Mpa
nu=0.3; % Possion's Ratio
% Plate dimension --> mm
a=500;
b=500;
D=E*t^3/(12*(1-nu^2));
z=t/2;
w=(4*P/(pi^4*a*b*D))*symsum(symsum((-1)^(m+n)*sin((2*m+1)*pi*x/a)*...
sin((2*n+1)*pi*y/b)/(((2*m+1)/a)^2+((2*n+1)/b)^2)^2,n,0,10),m,0,10);
sigmax=-E*z/(1-nu^2)*(diff(w,x,2)+nu*diff(w,y,2));
sigmay=-E*z/(1-nu^2)*(nu*diff(w,x,2)+diff(w,y,2));
xx=0:10:500;
yy=0:10:500;
[xx,yy]=meshgrid(xx,yy);
w=subs(w,{x,y},{xx,yy});
max(max(w))
figure(1)
mesh(xx,yy,w)
figure(2)
sigma_x=subs(sigmax,{x,y},{xx,yy});
mesh(xx,yy,sigma_x)
figure(3)
sigma_y=subs(sigmay,{x,y},{xx,yy});
mesh(xx,yy,sigma_y)
wmax=max(max(w));
sigma_x_max=max(max(sigma_x));
sigma_y_max=max(max(sigma_y));
disp(wmax)