Optical pulse compression using dispersion decreasing fibres

We get 4 ODEs from non linear schrodinger equation which has pulse parameter such amplitude, pulse width, chirp and phase. I need generate a MATLAB code to solve these set of 4 ODEs and then that will be used to simulate the pulse compression using dispersion decreasing fibres . I need simulation part as well. Please help