surf画三维图问题

clear
clc
%syms t;
m=1;n=1;

for i=0:0.01:0.5;
for d=0.5:0.01:1;

D(m,n)=(1-sqrt(1-2*i))/2;
pb(m,n)=(2*D(m,n)*d - d + 1)^2/(4*(d + 1)) + (d*(2*D(m,n) + d - 1)^2)/(4*(d + 1));
n=n+1;
end
n=1;
m=m+1;
end
[d,i]=meshgrid(0.5:0.01:1,0:0.01:0.5);
surf(i,d,pb)
xlabel('i')
ylabel('d')
zlabel('pb')

如上述代码所示，绘出的三维图为什么会有一些多余的黑色线？而且都是首尾相连，不知道如何处理，求大神帮忙解决，非常感谢。