I've seen more than 5 videos on this topic. This is by far the best video that I've come across. Thanks for the wonderful explanation. Now I have crystal clarity :)
Good morning, I have a question to ask you: how can I adapt this system of equations in my case ?: d Is / dz = g_r * Ip * Is -alpha_s * Is; d Ip / dz = - (fp / fs) * g_r * Ip * Is -alpha_p * Ip; with g_r = 1 * 10 ^ -13 (coeff Raman), lambda_p = 1550nm, lambda_s = 1450 nm, alpha_s = alpha_p = 0.2 dB / km
Change ls to y(1), lp to y(2), d ls / dz to dy(1), and d lp / dz to dy(2). You didn't mention fp or fs, but you can define all of your constants in the functions as well.
Nothing stands out to me from your code, though I don't know the subject matter. Double check your initial conditions and the time you're integrating over.
I've responded as clearly and thoroughly as I can without going and doing a good amount of research on my own. If you have a specific question, I'll be happy to do my best to answer.