function dydt = odefunc(t, y, strfunc)

% Функция правой части дифференциального уравнения
% определяется строкой strfunc, которая задает функцию переменных t, y 
% Создается из строки при помощи функций inline и feval 

% Если нет переменных t, y, то их формально добавляем, 
% чтобы функция считалась зависящей от t, y

strfunction = strcat(strfunc, '+ 0*t + 0*y') ;

F = inline(strfunction) ;
dydt = feval(F, t, y) ;
