Help with a double series?

23 views (last 30 days)
DIANA CARTLON on 28 Nov 2019
Answered: Jyothis Gireesh on 30 Dec 2019
I am trying to write function for the equation, can someone tell me what i am doing wrong
clear,clc
E=input('Please enter the Modulus of Eleasticity, E\n');
v=input('Please enter the Poisson Ratio, v\n');
rec_plate_pt(a,b,P,xo,yo,x,y);
D=(E*h^3)/(12*(1-v^2))
c=(4*P)/(pi^4*a*b*D)
function w = rec_plate_pt(a,b,P,xo,yo,x,y)
s=0;
for m=1:Inf
for n=1:Inf
s=c*(sin((m*pi*xo)/a)*sin((n*pi*yo)/b)/((m^2/a^2+n^2/b^2).^2)*(sin((m*pi*x)/a))*(sin(n*pi*y)/b))
format long
fprintf('%d',w);
end
end
end

1 Comment

Walter Roberson on 28 Nov 2019
Those are infinite series. You should use symsum()

Jyothis Gireesh on 30 Dec 2019
As Walter mentioned in the previous comment the double summation may be implemented using a combination of symbolic variables and the symsum()” in MATLAB instead of defining a function.
One possible implementation is given below
syms Wc(x,y) m n;
Wc(x,y) = c*symsum(symsum(sin(m*pi*x0/a)*sin(n*pi*y0/b)*sin(m*pi*x/a)*sin(n*pi*y/b)/(m^2/a^2 + n^2/b^2)^2,n,1,Inf),m,1,Inf);