Поле наклона дифференциального уравнения и решение Коши

Question

Дано дифференциальное уравнение

$2y'=\frac{yx}{x^2+4}+\frac{x}{y}$ 

• Я должен нарисовать поле наклона в прямоугольнике с именем $P$, содержащий точку (2,1).
• Тогда в соответствующем интервале я должен найти решение Задача Коши для заданного дифференциального уравнения с начальным уравнение $y(x_{0})=y_{0}$ where $(x_{0},y_{0})$ определяется нажатие на прямоугольник $P$
• и в том же прямоугольнике $P$ график найденного приближения задачи Коши с данным исходным уравнением, приведенным выше.

Вот мое решение:

function Plotslope
x=-5:0.6:5;
y=-6:0.6:6;
delta=0.2; 
hold on
axis([-5,5,-6,6])
daspect([1,1,1])

for k=1:length(x)
        for m=1:length(y)
          eps=delta/(sqrt(1+ff(x(k),y(m))^2));
          plot([x(k)-eps, x(k)+eps],...
         [y(m)-eps*ff(x(k),y(m)),...
           y(m)+eps*ff(x(k),y(m))],'k');
        plot(x(k),y(m),'k.','LineWidth',0.2) 
        end
end

[x0,y0]=ginput(1);
plot(x0,y0,'bo')
[T,Y]=ode45(@ff,[x0,5],y0);
[T1,Y1]=ode45(@ff,[x0,-6],y0);
plot(T,Y,'r',T1,Y1,'r')

function z=ff(x,y) 
    z=(y*x)/(2*(x^2+4))+x/2;
end
end

У меня вопрос, как я могу написать решение для 3-й пули: «и в том же прямоугольнике $ P $ график найденного приближения задачи Коши с данным исходным уравнением, приведенным выше». или я это уже решил?