davidrpugh / davidrpugh/pyCollocation

Example Suggestions: Heat quation, Lorenz equations, Lotka-Volterra equations and Sturm-Liouville Equation

Open
#46 7 comments 0 reactions 0 assignees View on GitHub
Dominant language
Python
Stars
16
Forks
8
PR merge metrics
No merged PRs in 30d

Description

Hi,

I'm quite the fan of this Python module, and I wanted to suggest the following examples:
- The [Heat equation](https://en.wikipedia.org/wiki/Heat_equation), that is:

![](https://wikimedia.org/api/rest_v1/media/math/render/svg/6544d3fe02c45c7f36beb824d68073274d7b60b3)
- The [Lorenz equations](https://en.wikipedia.org/wiki/Lorenz_system):

![](https://wikimedia.org/api/rest_v1/media/math/render/svg/5f993e17e16f1c3ea4ad7031353c61164a226bb8)

Would be best to solve it with the original parameters: ![sigma](https://wikimedia.org/api/rest_v1/media/math/render/svg/f3e7e6795cde44bcd0913e18917d7de69575cca5), ![beta](https://wikimedia.org/api/rest_v1/media/math/render/svg/7639fbf77fddf40dc63b185e40a4676ac55d30c5) and ![rho](https://wikimedia.org/api/rest_v1/media/math/render/svg/e3a9e7187d23db35f1c1c1e80d3ce289d0d3092d). `x[0] = y[0] = z[0] = 1` on t ∈ [0,100].
- The [Lotka-Volterra equations](https://en.wikipedia.org/wiki/Lotka%E2%80%93Volterra_equations), that is:

![](https://wikimedia.org/api/rest_v1/media/math/render/svg/79752d662d4760abcc84c6f0bb94d708f17ff442)
- The [Sturm-Liouville Equation](https://en.wikipedia.org/wiki/Sturm%E2%80%93Liouville_theory), specifically the example:

![](https://wikimedia.org/api/rest_v1/media/math/render/svg/2e036052d7599dbf0a33db64d401039096a022c3)

where lambda is the eigenvalue of the equation. This problem is integrated over the range x ∈ [0,∞], although a range of `[0,100]` is good enough to get the first 10 or more eigenvalues. Plus it has the boundary conditions y(0) = y(∞) = 0. This equation has the exact solution:

![](https://wikimedia.org/api/rest_v1/media/math/render/svg/7d07e6b5511e18e50b805c790f3ff82ebb1e5bc8)

where Ai is the [Airy Ai function](https://en.wikipedia.org/wiki/Airy_function). The eigenvalues are the zeros of the Airy Ai function, with their negative sign replaced with a positive sign.

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.