A cosmological boundary value problem
显示 更早的评论
Dear Community,
Recently I have encountered the following system of odes:
where κ and 
known constants with the boundary conditions 
, 
. I would like to construct a proper code by using bvp4c or bvp5c but I am not sure how to handle the boundary conditions.
known constants with the boundary conditions , . I would like to construct a proper code by using bvp4c or bvp5c but I am not sure how to handle the boundary conditions. I would be very grateful for any help you could give me. Have a nice day!
6 个评论
Jan
2021-6-7
dSdt = S / (3 * G) must fail at G(0) = 0. Then the derivative is not finite and the integration will reply infinite results. This is not a probölemk of Matlab, but a mathematical one.
James Tursa
2021-6-7
编辑：James Tursa
2021-6-7
Good point. In fact all three equations have 0/0 at the boundary condition.
Szilárd Zsóka
2021-6-7
Both of you are right. Numerically I set 0 to 0.01. The problem has an asymptotic power series solution. As 
then 
and 
. Our goal is to find asymptotic solution for the system which cannot be approximated by power series. I have several initial guesses but I am not very good at numerical calculations.
then and . Our goal is to find asymptotic solution for the system which cannot be approximated by power series. I have several initial guesses but I am not very good at numerical calculations.
James Tursa
2021-6-7
Seems like you still have a mathematical problem because asymptotically the rhs has 1/sqrt(t). Sure you can arbitrarily start the problem at 0.01, but does this really have meaning if the actual math doesn't match and tends to infinity? Where are these equations coming from?
Szilárd Zsóka
2021-6-8
Szilárd Zsóka
2021-6-8
By considering the power series solutions and the fact that G can be used as new independent variable I have reformulated the system into the following form:
with 
and 
boundary conditions. Can it solved with https://www.mathworks.com/help/matlab/math/solve-bvp-with-singular-term.html? What do you think?
and boundary conditions. Can it solved with https://www.mathworks.com/help/matlab/math/solve-bvp-with-singular-term.html? What do you think? 回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
2021-6-7
2021-6-8
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
