c3

//@version=5
indicator(title='2-Sine Wave Decomposition with Schumann Resonance', precision=12)

wave_height = input(100)
wave_duration1 = input(161)
wave_duration2 = input(217)
wave_phase1 = input(0)
wave_phase2 = input(0)

n = bar_index

f_sine_wave(_wave_height, _wave_duration, _wave_phase, _close) =>
    _pi = 3.14159265359
    _w = 2 * _pi / _wave_duration
    _sine_wave = _wave_height * math.sin(_w * n + _wave_phase)
    _sine_wave + _close

// Define function to generate Schumann resonance frequency
f_schumann_resonance() =>
    // Schumann resonance frequency varies from 7.83 Hz to 33.8 Hz
    _f = (33.8 - 7.83) * math.pow(math.sin(2 * math.pi * n / 80640), 2) + 7.83
    _f

close_price = input(close)

sine_wave1 = f_sine_wave(wave_height, wave_duration1, wave_phase1, close_price) + f_schumann_resonance()
sine_wave2 = f_sine_wave(wave_height, wave_duration2, wave_phase2, close_price) + f_schumann_resonance()

min_trigger = ta.crossover(sine_wave1, sine_wave2)
max_trigger = ta.crossunder(sine_wave1, sine_wave2)

plot(series=sine_wave1, title='Sine Wave 1', color=color.green)
plot(series=sine_wave2, title='Sine Wave 2', color=color.red)
plot(series=min_trigger ? sine_wave1 : na, title='Min Trigger', style=plot.style_cross, linewidth=3, color=color.black)
plot(series=max_trigger ? sine_wave1 : na, title='Max Trigger', style=plot.style_cross, linewidth=3, color=color.blue)