MATLAB Answers

Help with a double series?

9 views (last 30 days)
DIANA CARTLON
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
a=input('Please enter the dimension, a\n');
b=input('Please enter the dimension, b\n');
h=input('Please enter the thickness, h\n');
E=input('Please enter the Modulus of Eleasticity, E\n');
v=input('Please enter the Poisson Ratio, v\n');
P=input('Please enter the Load, P\n');
xo=input('Please enter the Position, xo\n');
yo=input('Please enter the Position, yo\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
Walter Roberson on 28 Nov 2019
Those are infinite series. You should use symsum()

Sign in to comment.

Answers (1)

Jyothis Gireesh
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);

  0 Comments

Sign in to comment.