Simulation of incoherently coupled multicomponent solitons on MATLAB

We show that incoherently coupled multicomponent beam evolution equations in the low-amplitude regime in biased photorefractive-photovoltaic crystals can exhibit the analytical solutions of multicomponent solitons, which are due to both the spatially nonuniform screening of the external electric field and the photovoltaic effect. These multicomponent solitons can be obtained for N self-trapped mutually incoherent optical beams, where N is an arbitrary whole number. When the photovoltaic effect is neglected, their nonlinear evolution equations change into that for multicomponent screening solitons, and these multicomponent solitons change into multicomponent screening solitons. When the external bias field is absent, their nonlinear evolution equations change into that for multicomponent photovoltaic solitons in closed and open circuits, and these multicomponent solitons change into multicomponent photovoltaic solitons in closed and open circuits.