Davie defined a Levy variant and the combination of single random variables to ensure that the diffusion matrix did not degenerate. The use of the method proposed by Davie, which is a combination of the Euler method and the exact combination, was investigated for applying the degenerate Levy diffusion approach to $(B_{ik}(Y))$. We use certain degenerate conditions of diffusion which contribute to order convergence. We also show MATLAB codes to apply the integrated solution to an SDE and observe a convergence behavior. We also evaluate the agreement between the theoretical values and the MATLAB numerical example.