Sturm-Liouville eigenvalues calculation

17 visualizaciones (últimos 30 días)
Luigi Romano
Luigi Romano el 12 de Jun. de 2020
Respondida: Bjorn Gustavsson el 13 de Jun. de 2020
Hallo,
I am trying to solve a Sturm-Liouville eigenvalues problem. From the Robin BC, I get the following trascendental equation for the variable x:
tan(x*L/(2*S)) = 4*x*tau*S/((x*tau)^2+(tau*c*V)^2-4*S^2),
where L, S, tau, c and V are parameters. I should have ideally infinite solution (I would like to compute the first 100 solutions), but I cannot find them. I was trying with fsolve:
F= @(x) tan(x/(2*S)*l)-4*x*S*tau/((x*tau)^2+(tau*c*V)^2-4*S^2);
for i=1:1:1000
j = i-1;
a(i) = j*1000;
x0(i) = a(i);
mu(i) = fsolve(F,x0(i));
end
but I don't manage. Do you have any idea?
Thank you in advance
  1 comentario
darova
darova el 13 de Jun. de 2020
What is the problem? Can't you find a solution or what?

Iniciar sesión para comentar.

Respuestas (1)

Bjorn Gustavsson
Bjorn Gustavsson el 13 de Jun. de 2020
Have you checked the file-exchange contribution: discrete-orthogonal-polynomial. It solves this kind of problems numerically.

Categorías

Más información sobre Ordinary Differential Equations en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by