How to use symsum command for a function of two variables and two indices?
8 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Sania Nizamani
am 17 Jul. 2022
Kommentiert: Sania Nizamani
am 18 Jul. 2022
Can anyone help to write a code in MATLAB for the expression given below:
I have tried with the following but it doesn't work: syms x y k1 k2 M_series=symsum(gamma(k1+1).*gamma(k2+1).*(x.^k1).*(y.^k2)./(factorial(k1).*factorial(k2).*gamma(k1+k2+1)),k1,0,Inf,k2,0,Inf);
3 Kommentare
Torsten
am 17 Jul. 2022
So since factorial(k1) = gamma(k1+1) and factorial(k2) = gamma(k2+1), these terms cancel out and you are left with
sum_i sum_j 1/(i+j)! * x^i * y^j
Is this correct ?
Akzeptierte Antwort
Torsten
am 17 Jul. 2022
syms x y i j
assume(i,'integer')
assume(j,'integer')
A = symsum(symsum(x^i*y^j/factorial(i+j),i,0,Inf),j,0,Inf)
x0 = 0:0.2:1;
y0 = 0:0.25:1;
[x0,y0] = meshgrid(x0,y0);
[n,m] = size(x0);
for i = 1:n
for j = 1:m
A0(i,j) = double(subs(A,[x,y],[x0(i,j),y0(i,j)]));
end
end
surf(x0,y0,A0)
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Calculus finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!