my differential equation general solution form requires very high or low values, exceeding double precision capability, while I know its solution is finite

1 view (last 30 days)
hosein Javan
hosein Javan on 17 Aug 2020
Edited: hosein Javan on 19 Aug 2020
as a brief explanation, I'm trying to solve a first order bessel equation which is part of a more complicated PDE problem, whose solution is as following form:
now we have a bunch of these functions and not one, so there are for example 10 integral coefficients like c1(1), c1(2),... and c2(1),c2(2),....
now when I apply the boundary conditions, a set of linear equations are made by which I can calculate c1, c2,... . in the form "A*C=B".
now as x is near zero, I1(x) goes to zero and K1(x) goes to inf, and as x goes up, I1(x) goes to inf and K1(x) goes to zero. given that y(x) is finite and something less than 0.07, c1 and c2 go to inf or zero depending on their factor. and here the problem begins. working with very high and low values decrease accuracy and might lead to inf or zero.
any idea to overcome this issue?
  10 Comments
hosein Javan
hosein Javan on 19 Aug 2020
Hello again Mr. Goodmanson. thanks for your patience. yes the boundary condition is somehow like that. in fact it states that and must be continuous at each y1,...,y5. where "k" and "m" are constants for each "n" and has been defined before. and the regions are distinguished by y1,...,y5 boundaries. p is the frequency across only "x" and not "y", so for each region it is the same.
Mr. Goodmanson what about this idea: the input of bessel function dimension is radian which is "dimensionless", however I was wondering if we change that unit, what change should I do to the equations?

Sign in to comment.

Answers (0)

Community Treasure Hunt

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

Start Hunting!

Translated by