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')