1:
2:
3:
4:
5:
6:
7:
8:

x=-5:.2:5;
y=sin(x.^2);
plot(x,y);
d=sign(diff(y));
ind=find(d~=[d(1) d(1:end-1)]);
hold on;
plot(x(ind),y(ind),'ro');
hold off;