Tsunami is a natural disaster that have resulted in dreadful damages over time. Extensive researches have been conducted to scrutinize and counteract the natural hazard using three major research components which are: field monitoring, laboratory tests, and numerical methods. However, laboratory tests are high-priced and arduous. Numerical simulation overcomes these drawbacks and can be utilized in collaboration with laboratory tests. Recently, newly introduced meshless Lagrangian particle method called Smoothed Particle Hydrodynamics (SPH) has gained attention. In this paper, SPH method has been employed to simulate tsunami. A SPH code is developed from scratch. To validate the code, a traditional dam break simulation is conducted. Lastly, a tsunami model is simulated using the developed SPH code and compared with past experimental data. The results indicate that the code is in accordance with previous experimental data and numerical simulation. Whereby, there’s been a slight deviation arises in tsunami simulation. The velocity of the code is relatively less to that of the experimental data. Such inconsistencies could emerge due to a number of reasons, i.e. the choice of the SPH parameters and model simplification. Generally, the developed SPH code had a satisfactory performance to model tsunami and dam-break problem.