【MATLAB】请问这个作图语句哪里有问题呢,为什么出现空图像?

【MATLAB】请问这个作图语句哪里有问题呢,为什么出现空图像?
a=1;b=1;u0=1;
pi=3.14159265;
da=a/100;db=b/100;
[x,y]=meshgrid(0:da:a,0:db:b);
u=0*x;N=20;
for I=0:N
n=2*I+1;
c=n*pi/a;d=n*sinh(c*b);
u=u+sinh(c*x).*sinh(c*y)/d;
end
u=u*4*u0/pi;
contour(x,y,u,20);
新手刚开始学MATLAB,希望大神解答!