Version (1.06 MB) by Haihong
This is an n-body gravitational simulation, where Verlet algorithm and Octree method are applied.
Updated 1 Jan 2016

View License

1. You may run the scripts to see the examples.
2. You may also use your own case.
3. Your own case should be a struct array whose elements represent bodies’ initial states. Each element (body) contains 4 fields: mass (a positive scalar), position (a 1-by-3 vector, relative to the origin point), velocity (a 1-by-3 vector), and acceleration (a 1-by-3 vector).
You may check the examples in these scripts and in the folder TestCase.
Visualization-related settings are in the class Viewer. You may modify the code to see what happens :)
(a) marker size, (b) axis limit, (c)choice of whether fixed axis is applied, (d) choice of whether minor grid is displayed, (e) choice of whether stars’ trails are displayed, and (f) choice of whether star 1 is highlighted.
4. Note that this program cannot handle collisions - they just "pass through" each other (not true in reality). The script Script2starCollision shows what would happen in this situation.
5. Just a prototype, may has some bugs. Ctrl+C to terminate the program.

Cite As

Haihong (2024). GalaxyModel (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2015a
Compatible with any release
Platform Compatibility
Windows macOS Linux
Find more on Gravitation, Cosmology & Astrophysics in Help Center and MATLAB Answers

Community Treasure Hunt

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

Start Hunting!
Version Published Release Notes

Included a technical note.
Added a new tag.
Added a snapshot:)
Modified the description.
Included a technical note.

Changed the description.

Minor changes.