Help with a double series?

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

[Walter Roberson]

Those are infinite series. You should use symsum()

Answer 1

[Jyothis Gireesh]

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