В результате международного о сотрудничества между ИКИР ДВО РАН и Национальным университетом Узбекистана имени Мирзо Улугбека опубликована статья, в которой рассматривается неявная разностная схема для уравнения Дуффинга с производной дробного переменного порядка типа Римана–Лиувилля. Рассмотрены вопросы устойчивости и сходимости неявной разностной схемы. Приведены тестовые примеры для обоснования теоретических результатов. С помощью правила Рунге результаты неявной схемы сравниваются с результатами явной схемы. Построены фазовые траектории и осциллограммы для осциллятора Дуффинга с дробной производной переменного порядка типа Римана–Лиувилля, обнаружены хаотические режимы с использованием спектра максимальных показателей Ляпунова и сечений Пуанкаре. Для исследования вынужденных колебаний построены поверхности добротности, амплитудно-частотные и фазочастотные характеристики. Результаты исследования показали, что неявная разностная схема показывает более точные результаты, чем явная.