x = linspace(0,1 ,51)';
y = linspace(0,1,51)';
[X,Y] = meshgrid(x, y);
Z = exp(-2*pi^2*Y).*sin(2*pi*X);
surf(X, Y, Z);
ylabel('t');
xlabel('x');
zlabel('u(x,t)')
Lets have a look, so you can understand why that line. You wrote this:
Z=exp(-4.*(pi^2)*0.5*Y)*sin(2*pi*X);
^ element-wise multiply, but pointless.
But what is the point of element-wise multiplication with 4? 4 is a scalar number and so is pi, it serves no purpose whatsoever to use element-wise multiplication with them. So what parts of that calculation are not scalar? Answer: X and Y. Therefore any operation involving X and Y, or any other arrays derived from them will need to use element-wise operations. In this case this means the exp(..) and sin(..) terms will need to be multiplied using element-wise multiply. And thus I showed you (twice so far) this, which uses element-wise multiply for the terms that are arrays (and I did not use it on the scalar values like you did):
exp(-2*pi^2*Y).*sin(2*pi*X)
^ element-wise multiply, required.
or if you really want to have that 4 and 0.5:
exp(-4*pi^2*0.5*Y).*sin(2*pi*X)
You have now been told this multiple times:
Best Answer