Equation is convergent but unable to extract the output
    4 views (last 30 days)
  
       Show older comments
    
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. 
0 Comments
Accepted Answer
  Torsten
      
      
 on 20 May 2022
        
      Edited: Torsten
      
      
 on 20 May 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)])
0 Comments
More Answers (0)
See Also
Categories
				Find more on Symbolic Math Toolbox in Help Center and File Exchange
			
	Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
