Решение дифференциального уравнения методом Эйлера и методом Рунге-кутта

function [Y]=reshenie(f,a,b,a1,b1,ya,h) 
% f - функция правой части
% a и b - границы интервала
% a1, b1 - коэффициенты уравнения
% ya - начальное условие
% h - шаг
cc=1; 
H=0; 
while cc 
    % создание меню
    p=menu('Решение ДУ','Метод Эйлера','Метод Рунге-Кутта','Выход');
 
    % метод Эйлера
    if p==1 
        H=H+1; 
        fprintf(1,'\nЭксперимент %u\n', H) 
        disp('Метод Эйлера') 
        if H==0 
            p=0; 
        else p=1;
        end 
        [Y]=Eiler(f,a,b,a1,b1,ya,h,p) % вызов функции Eiler                               
        disp('Нажмите Enter') 
        pause 
    end 
 
    % Метод Рунге-Кутта
    if p==2 
        H=H+1; 
        fprintf(1,'\nЭксперимент %u\n', H) 
        disp('Метод Рунге-Кутта') 
        if H==0 
            p=0; 
        else p=1; 
        end 
        [Y]=RK(f,a,b,a1,b1,ya,h,p); % вызов функции RK               
        Y 
        disp('Нажмите Enter')
        pause 
    end 
 
    % Выход
    if p==3 
        return 
    end 
end