Created by **Nils Berglund** and optimized by **Marco Mancini** C code for videos on YouTube Channel https://www.youtube.com/c/NilsBerglund Below are parameter values used for different simulations, as well as initial conditions used in function animation. Some simulations use variants of the published code. The list is going to be updated gradually. #### 31 May 2026 - A four-sided parabolic resonator with Neumann boundary conditions ### **Program:** `wave_billiard.c` **Initial condition in function `animation()`:** ``` /* add oscillating waves */ for (source = 0; source < N_SOURCES; source++) { wave_source_x[source] = 0.0; wave_source_y[source] = 0.0; source_amp[source] = INITIAL_AMP; } source_periods[0] = OSCILLATING_SOURCE_PERIOD; for (source = 0; source < N_SOURCES; source++) { dperiod = source_periods[source]; phase = i - (int)(dperiod*(double)((int)((double)i/dperiod))); if ((ADD_OSCILLATING_SOURCE)&&(phase == 1)&&(i