Hello! How do I solve ODE numerically. I tried, but I fall into a singularity when t is greater than 0.01. In the problem in question, when t = 0.03, my f(t) is approximately +0.003. I used several numerical methods, including CK45, but I was unable to get the solution for f(t) with any of them. I can't get out of the singularity!!!
Please help me!!!
Hello,
a singularity is also suspected in Wolfram Mathematica , see the attached file.
I suggest you follow MATLAB for numerical ODE/PDE, e.g., https://au.mathworks.com/help/matlab/math/solve-transistor-pde.html
https://au.mathworks.com/help/matlab/ref/pdeval.html
https://au.mathworks.com/help/matlab/math/partial-differential-equations.html#:~:text=The%20MATLAB%C2%AE%20PDE%20solver,change%20with%20respect%20to%20time.
Please check the equation. I have implemented a code to solve it and would gladly send it to you, but the behavior is so strange (blows up almost immediately) that I suspect there's an error. You are subtracting numbers that are almost zero and dividing by 0.00002 to get df/dt so that something appears to be wrong. Something must be happening around f=0.007; otherwise, the term (0.007-f) wouldn't be showing up in the deq.
- Badges
- Science topic
- Similar topics
- Mathematics