%exm03132_2.m
t=linspace(0,3*pi,500);y=sin(t);  
%
z1=((t<pi)|(t>2*pi)).*y;										%<3>
w=(t>pi/3&t<2*pi/3)+(t>7*pi/3&t<8*pi/3);					%<4>
w_n=~w;																%<5>
z2=w*sin(pi/3)+w_n.*z1;											%<6>
shg,subplot(1,3,1),plot(t,y,':r'),ylabel('y')
subplot(1,3,2),plot(t,z1,':r'),axis([0 10 -1 1])
subplot(1,3,3),plot(t,z2,'-b'),axis([0 10 -1 1])  
disp('按任意鍵，繼續！'),pause
z=(y>=0).*y;														%<11>
a=sin(pi/3);
z=(y>=a)*a+(y<a).*z;												%<13>
clf
plot(t,y,':r');hold on;plot(t,z,'-b')
xlabel('t'),ylabel('z=f(t)'),title('逐段解析函數')
legend('y=sin(t)','z=f(t)',4),hold off 
disp('結束')