x=0:0.1:6*pi;
y = sin(x);
Y = (y>0).*y;
plot(x,y, ':', x,Y, 'LineWidth', 2 )
