Equation is convergent but unable to extract the output

1 vue (au cours des 30 derniers jours)
Sun Heat
Sun Heat le 20 Mai 2022
Modifié(e) : Torsten le 20 Mai 2022
hello friends,
i run the below code and found that the equation is 'convergent', but i am unable to extract the output.
clc;close all;clear all;
syms m
k=360; h=0.562; theta=35;
alpha_m=((-1)*360*sqrt(1-((m^2)*(sind(theta)^2))));
term=(k*(1-(exp(-2*alpha_m*h)))/alpha_m);
y = vpa(symsum(term, m, 2, Inf));
if isinf(y)
fprintf('divergent\n');
else
fprintf('convergent\n');
end
I used
syms m
k=360; h=0.562; theta=35;
alpha_m=((-1)*360*sqrt(1-((m^2)*(sind(theta)^2))));
term=(k*(1-(exp(-2*alpha_m*h)))/alpha_m);
y = symsum(term, m, 2, Inf);
AA=double(y) % but it didn't work
Please help me. Thanks in advance.

Réponse acceptée

Torsten
Torsten le 20 Mai 2022
Modifié(e) : Torsten le 20 Mai 2022
The real part of y seems to converge, the imaginary part doesn't.
k=360; h=0.562; theta=35;
alpha_m = @(m)((-1)*360*sqrt(1-((m.^2)*(sind(theta)^2))));
term = @(m)(k*(1-(exp(-2*alpha_m(m)*h)))./alpha_m(m));
m = 2:10000;
y = cumsum(term(m))
plot(m,[real(y);imag(y)])

Plus de réponses (0)

Catégories

En savoir plus sur Symbolic Math Toolbox dans Help Center et File Exchange

Produits

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by