MATLAB: Mesh plot for a series equation

mesh plot

Error using mesh (line 76),Z must be a matrix, not a scalar or vector. Can anyone please sort out for me the possible solution? Thanks in advance

L=2;
W=2;
x=0:0.2:2;
y=0:0.2:2;
T1=70;
T2=200;
n=1:0.2:1000;
[X,Y]=meshgrid(x,y);
for i=1:length(n);
    g_n=2/pi*sum(2/n(i)*sin(n(i)*pi*X/L).*sinh(n(i)*pi*Y/L)/sinh(n(i)*pi*W/L)));   
end 
Z=(T2-T1)*sum(g_n)+T1;
mesh(X,Y,Z)

Best Answer

See Roger Stafford’s Answer to: how can i write a matlab code for the summation below? please help.

Related Solutions

MATLAB: I have a problem with plotting a summation equation

Hi, you mean something like that ?

y=0:0.025:1;
x=0:0.05:0.7;
n=1:1:50;
[X,Y]=meshgrid(x,y);
for i=1:size(n,2)
    g_n=4*sum(1-cos(n(i)*pi)/(n(i)*pi)^3)*exp(-n(i)*pi*X).*sin(n(i)*pi*Y);
    g_xy(:,:,i)=g_n;
end
G=sum(g_xy,3);
mesh(X,Y,G)

MATLAB: How to perform symsum expansion with two known condition, in this case let say n and m and the n=m

You do not need the evalin(). You can use http://www.mathworks.com/help/symbolic/syms.html#inputarg_set

syms x y n m integer

You can call symsum() twice.

symsum( symsum(SomeExpression, n, 1, 5), m, 1, 5)

My tests show that under the assumption of integer for n, x, y, that m vanishes from the expression being summed, leaving

400*(-1)^(1+n) * pi * sinh(n * pi - (1/2)*pi)^2 * x * y / (2*n-1)

so the sum over m from 1 to 5 would just be 5 times the sum for n.

Best Answer

Related Solutions

MATLAB: I have a problem with plotting a summation equation

MATLAB: How to perform symsum expansion with two known condition, in this case let say n and m and the n=m

Related Question