Differential equation I can't solve

Click For Summary
SUMMARY

The discussion centers on solving the differential equation \(\frac{dx}{dt} = k \sqrt{t - x^2}\) with a positive constant \(k\). The original poster concludes that the equation lacks a closed-form solution and suggests using numerical methods, specifically Gnumeric, for approximation. They report that using Mathematica's NDSolve yields \(x = 0.60635\) at \(t = 1\) for \(k = 1\). Additionally, approximate solutions via Taylor series are provided for small and large \(t\), with a notable intersection at \(t = 2.354\).

PREREQUISITES
  • Understanding of differential equations, particularly first-order equations.
  • Familiarity with numerical methods for solving differential equations.
  • Knowledge of Taylor series expansions and their applications.
  • Experience with software tools like Mathematica and Gnumeric for numerical analysis.
NEXT STEPS
  • Explore numerical methods for differential equations using Gnumeric.
  • Learn about Mathematica's NDSolve function for solving differential equations.
  • Study Taylor series expansions in detail, focusing on their convergence properties.
  • Investigate the implications of varying the constant \(k\) in the differential equation.
USEFUL FOR

Mathematicians, physicists, and engineers who are dealing with complex differential equations and seeking numerical solutions. This discussion is particularly beneficial for those interested in numerical analysis and approximation techniques.

Irid
Messages
207
Reaction score
1
Hi,
I was solving my own-invented physical problem, and obtained a differential equation
[tex]\frac{dx}{dt} = k \sqrt{t-x^2}[/tex]
with k positive constant. I wish to solve it. I'm unable perform separation of variables here, and power-series method works poorly, too. I conclude that this equation does not have any closed-form or familiar series solution, so the best thing that I could do is use numerical methods [of an Excel type, but Gnumeric is my choice :)]. It turned out that with k=1, when t=1, x=0,6, exacty. This gives a clue that there might exist some sort of a solution (or it might not be any kind of a clue, I'm just being silly). How do I find it? How do I solve this equation?
 
Physics news on Phys.org
There's no exact solution. I'm assuming you're using a boundary condition of x=0 at t=0. Then, Mathematica's numerical routine NDSolve gives x=0.60635 at t=1.

You can solve it approximately by Taylor series for small t and large t. At small t,
[tex]x = {{\textstyle{2\over3}}t^{3/2}\left(1-{\textstyle{2\over21}}t^2<br /> +{\textstyle{10\over2079}}t^4 + \ldots\right)[/tex]
At large t,
[tex]x = t^{1/2}\left(1-{\textstyle{1\over8}}t^{-2}<br /> -{\textstyle{13\over128}}t^{-4} + \ldots\right)[/tex]
These cross at t=2.354, where they differ from the numerical solution by a worst-case 1%.

This is for k=1. But if [tex]x=f(t)[/tex] solves the equation for k=1, then [tex]x=k^{-1/2}f(kt)[/tex] solves it for general k.
 
Last edited:

Similar threads

Replies
2
Views
2K
Replies
13
Views
2K
Replies
7
Views
3K
  • · Replies 12 ·
Replies
12
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
2
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K