x=0:pi/100:4*pi;
y1=sin(x);
y2=cos(x);
k=find(abs(y1-y2)<0.001);
x1=x(k);
y3=sin(x1);
plot(x,y1,'g--',x,y2,'b',x1,y3,'rp')
x=0:pi/100:4*pi;
y1=sin(x);
y2=cos(x);
k=find(abs(y1-y2)<0.001);
x1=x(k);
y3=sin(x1);
plot(x,y1,'g--',x,y2,'b',x1,y3,'rp')