Построить график функции с фишками с координатами

function PltWpt
x=-pi:pi/96:pi;
f=@(x)x.^3+x.^2.*cos(x)+x.*sin(x)+1;
y=f(x);
ha=axes;
line(x,y,'LineWidth',2,'ButtonDownFcn',@MyFun)
axis tight
set(ha,'ButtonDownFcn',@MyFun)
XL=get(gca,'XLim');
YL=get(gca,'YLim');
line(XL,[0,0],'Color','k','LineStyle',':','ButtonDownFcn',@MyFun)
line([0,0],YL,'Color','k','LineStyle',':','ButtonDownFcn',@MyFun)
hl=[];
ht=[];
hp=[];
  function MyFun(src,eventdata)
  delete([hl,ht,hp])
  CP=get(ha,'CurrentPoint');
  x0=CP(1,1);
  y0=CP(1,2);
  y1=f(x0);
  if sign(y0)~=sign(y1)
    y0=-y0;
  end
  if abs(y1)>abs(y0)
    y1=y0;
    x1=fzero(@(x)f(x)-y0,0);
  else
    x1=x0;
  end
  hl=line([0,x0,x0],[y0,y0,0],'Color','r','LineStyle',':');
  ht=text(x1,y1,['(',num2str(x1),';',num2str(y1),')'],'VerticalAlignment','top');
  hp=line([x1,x1],[y1,y1],'Marker','o','MarkerSize',4,...
    'MarkerEdgeColor','r','MarkerFaceColor','r');
  end
end

function PltWpt1
x=-pi:pi/96:pi; % вектор абсцисс
f=@(x)x.^3+x.^2.*cos(x)+x.*sin(x)+1; %указатель на анонимную функцию
y=f(x); % вектор ординат
ha=axes; %строим оси и получаем с них укзатель
line(x,y,'LineWidth',2,'ButtonDownFcn',@MyFun) %строим линию по ранее найденым
% векторам.
% Также учитываем что пользователь может нажать на этой линии. Если такое
% произойдёт, то выполнится функция MyFun
axis tight %устанавливаем диапазоны на осях равные диапазонам изменения 
% абсцисс и ординат нарисованой кривой
set(ha,'ButtonDownFcn',@MyFun) %Если пользователь нажмёт мышкой на белом полотне осей,
% то выполнится функция MyFun
XL=get(gca,'XLim'); % XL=[min(x),max(x)] c той лишь разницей
% что достал со свойства объекта оси
YL=get(gca,'YLim'); % YL=[min(y),max(y)]
line(XL,[0,0],'Color','k','LineStyle',':','ButtonDownFcn',@MyFun) % ось X
line([0,0],YL,'Color','k','LineStyle',':','ButtonDownFcn',@MyFun) % ось Y
box on % рамка вокруг осей
hlx=[];% Пустые массивы специально чтобы функции delete не было чего удалять
hly=[];% при первом щелчке мышкой
ht=[];
hp=[];
  function MyFun(src,eventdata)
  delete([hlx,hly,ht,hp]) % удаляем старые перпендикуляры, текст, точку
  CP=get(ha,'CurrentPoint'); % узнаём куда ж ткнул пользователь
  x0=CP(1,1);
  y0=CP(1,2);
  y1=f(x0);% расчитываем ординату кривой для выбраной пользователем абсциссы
  if sign(y0)==sign(y1) %если кривая по одну сторону оси х с выбраной пользователем точкой 
    hlx=line([0,x0],[y0,y0],'Color','r','LineStyle',':','ButtonDownFcn',@MyFun);%перепендикуляр к оси x
    if abs(y1)>abs(y0) %если кривая дальше от оси чем выбраная точка  
      hly=line([x0,x0],[0,y1],'Color','r','LineStyle',':','ButtonDownFcn',@MyFun);%то перепендикуляр опускаем с кривой
    else %если выбраная точка дальше от оси чем кривая
      hly=line([x0,x0],[0,y0],'Color','r','LineStyle',':','ButtonDownFcn',@MyFun);%то перепендикуляр опускаем с выбраной точки
    end
  else %если кривая и выбранная пользователем точка с разных сторон оси х 
    hlx=line([0,x0],[y0,y0],'Color','r','LineStyle',':','ButtonDownFcn',@MyFun);%перепендикуляр к оси x
    hly=line([x0,x0],[y0,y1],'Color','r','LineStyle',':','ButtonDownFcn',@MyFun);%перепендикуляр к оси y
  end
  ht=text(x0,y1,['(',num2str(x0),';',num2str(y1),')'],'VerticalAlignment','top','ButtonDownFcn',@MyFun);% подписываем точку на кривой
  hp=line([x0,x0],[y1,y1],'Marker','o','MarkerSize',4,...
    'MarkerEdgeColor','r','MarkerFaceColor','r','ButtonDownFcn',@MyFun); %собственно рисуем точку на кривой
  end
end

clear, clc
 
x = -pi:0.01:pi; % значения х
y = x.^3 + x.^2.*cos(x) + x.*sin(x) + 1; % значения функции
limt = [min(x) max(x) min(y) max(y)]; % пределы значений х и у
 
% рисуем:
hold on
plot(x,y,'b') % кривая
plot([min(x) max(x)],[0 0],'k',... % ось Х
    [0 0],[min(y) max(y)],'k') % ось Y
% сохраняем указатели(handle) объектов 
hp = plot(NaN,NaN,'or'); % точка тыка
hk = plot(NaN,NaN,'*r'); % точка на кривой
hpv = plot(NaN,NaN,'--r'); % вертикальный перпендикуляр
hph = plot(NaN,NaN,'--r'); % горизонтальный перпендикуляр
ht = text(NaN, NaN,'','BackgroundColor','g');
 
grid on % вкл. сетку
axis(limt) % устанавливаем границы окна
while 1 % бесконечный цикл
    [x0,y0] = ginput(1); % вводим одну точку
    set(hp,'XData',x0,'YData',y0);
    
    y1 = x0.^3 + x0.^2.*cos(x0) + x0.*sin(x0) + 1; % считаем координату на кривой
    set(hk,'XData',x0,'YData',y1);
    
    % координаты перпендикуляров:
    yu0 = min([0,y0,y1,min(y)]);
    yu1 = max([0,y0,y1,max(y)]);
    %yu1 = max(y).*(y0>0)+min(y).*(y0<0);
    xu = max(x).*(x0>0)+min(x).*(x0<0);
    set(hpv,'XData',[x0 x0],'YData',[yu0 yu1]);
    set(hph,'XData',[0 xu],'YData',[y1 y1]);
        
    % наносим подпись:
    str = ['(',num2str(x0,2),', ',num2str(y1,2),')'];
    set(ht,'Position',[x0+0.1, y1],'String',str)
end