Код MATLAB для диодной модели