I need to plot the light-current characteristic of laser (Pout vs I) by solving using ode45.
dy(1)= ((ni*I)/(q*V))-(N/tau)-(v_g*a*((N-N_tr)/(1+epsilon*N_p))*N_p);
dy(2)= (T*v_g*a*((N-N_tr)/(1+epsilon*N_p))*N_p) + (B_sp*T*R_sp) - (N_p/tau_p);
Pout = N_p*Vp*h*v*alpha_m*v_g;
how can I do that?