
function dy = odefunc3(t, yv, func)

% Вектор-функция правых частей системы дифференциальных уравнений
% определяется массивом строк func: F1(t, x, y, z), F2(t, x, y, z), F3(t, x, y, z), 
% Создается из строки при помощи функций inline и feval 

% Если нет переменных t, x, y, z, то правильнее их добавить формально:
strfunction1 = strcat(func(1,:), '+ 0*t + 0*x + 0*y + 0*z') ;
strfunction2 = strcat(func(2,:), '+ 0*t + 0*x + 0*y + 0*z') ;
strfunction3 = strcat(func(3,:), '+ 0*t + 0*x + 0*y + 0*z') ;

F1 = inline(strfunction1) ;
F2 = inline(strfunction2) ;
F3 = inline(strfunction3) ;

dy(1,:) = feval(F1, t, yv(1), yv(2), yv(3)) ;
dy(2,:) = feval(F2, t, yv(1), yv(2), yv(3)) ;
dy(3,:) = feval(F3, t, yv(1), yv(2), yv(3)) ;





