% [x,y,z]=ndgrid(linspace(-3,3,120));
% f=((-(x.^2) .* (z.^3) -(9/80).*(y.^2).*(z.^3)) + ((x.^2) + (9/4).* (y.^2) + (z.^2)-1).^3);
% P=patch(isosurface(x,y,z,f,0));
% set(P,'FaceColor','r');
%
% daspect([1 1 1]);
% view(3)
% colormap('flag')
% camlight('left')
% camlight('headlight')
% lighting phong
% xlabel('X')
% ylabel('y')
% zlabel('z')
% title('fighting')
%
% n=100;
% x=linspace(-3,3,n); %nirv
% y=linspace(-3,3,n);
% z=linspace(-3,3,n);
% [X,Y,Z]=ndgrid(x,y,z);
% F=((-(X.^2) .* (Z.^3) -(9/80).*(Y.^2).*(Z.^3)) + ((X.^2) + (9/4).* (Y.^2) + (Z.^2)-1).^3);
%
% isosurface(F,0)
% lighting phong
% caxis
% axis equal
% colormap('flag'); %颜色
% view([60 30]);
% a=['我','爱','你'];
% for i=1:5
% pause(0.5)
% if i>2
% text(70,25+4*i,65,[num2str(a(i))],'Fontsize',20,'color','m')
% else
% text(70,35+3.5*i,70,[num2str(a(i))],'Fontsize',20,'color','m')
% end
% end
% 最炫民族风
% fs = 44100; % sample rate
% dt = 1/fs;
% T16 = 0.125;
% t16 = [0:dt:T16];
% [temp k] = size(t16);
% t4 = linspace(0,4*T16,4*k);
% t8 = linspace(0,2*T16,2*k);
% [temp i] = size(t4);
% [temp j] = size(t8);
% % Modification functions
% mod4=(t4.^4).*exp(-30*(t4.^0.5));
% mod4=mod4*(1/max(mod4));
% mod8=(t8.^4).*exp(-50*(t8.^0.5));
% mod8=mod8*(1/max(mod8));
% mod16=(t16.^4).*exp(-90*(t16.^0.5));
% mod16=mod16*(1/max(mod16));
% f0 = 2*146.8; % reference frequency
% ScaleTable = [2/3 3/4 5/6 15/16 ...
% 1 9/8 5/4 4/3 3/2 5/3 9/5 15/8 ...
% 2 9/4 5/2 8/3 3 10/3 15/4 4 ...
% 1/2 9/16 5/8];
% % 1/4 notes
% do0f = mod4.*cos(2*pi*ScaleTable(21)*f0*t4);
% re0f = mod4.*cos(2*pi*ScaleTable(22)*f0*t4);
% mi0f = mod4.*cos(2*pi*ScaleTable(23)*f0*t4);
% fa0f = mod4.*cos(2*pi*ScaleTable(1)*f0*t4);
% so0f = mod4.*cos(2*pi*ScaleTable(2)*f0*t4);
% la0f = mod4.*cos(2*pi*ScaleTable(3)*f0*t4);
% ti0f = mod4.*cos(2*pi*ScaleTable(4)*f0*t4);
% do1f = mod4.*cos(2*pi*ScaleTable(5)*f0*t4);
% re1f = mod4.*cos(2*pi*ScaleTable(6)*f0*t4);
% mi1f = mod4.*cos(2*pi*ScaleTable(7)*f0*t4);
% fa1f = mod4.*cos(2*pi*ScaleTable(8)*f0*t4);
% so1f = mod4.*cos(2*pi*ScaleTable(9)*f0*t4);
% la1f = mod4.*cos(2*pi*ScaleTable(10)*f0*t4);
% tb1f = mod4.*cos(2*pi*ScaleTable(11)*f0*t4);
% ti1f = mod4.*cos(2*pi*ScaleTable(12)*f0*t4);
% do2f = mod4.*cos(2*pi*ScaleTable(13)*f0*t4);
% re2f = mod4.*cos(2*pi*ScaleTable(14)*f0*t4);
% mi2f = mod4.*cos(2*pi*ScaleTable(15)*f0*t4);
% fa2f = mod4.*cos(2*pi*ScaleTable(16)*f0*t4);
% so2f = mod4.*cos(2*pi*ScaleTable(17)*f0*t4);
% la2f = mod4.*cos(2*pi*ScaleTable(18)*f0*t4);
% ti2f = mod4.*cos(2*pi*ScaleTable(19)*f0*t4);
% do3f = mod4.*cos(2*pi*ScaleTable(20)*f0*t4);
% blkf = zeros(1,i);
% % 1/8 notes
% do0e = mod8.*cos(2*pi*ScaleTable(21)*f0*t8);
% re0e = mod8.*cos(2*pi*ScaleTable(22)*f0*t8);
% mi0e = mod8.*cos(2*pi*ScaleTable(23)*f0*t8);
% fa0e = mod8.*cos(2*pi*ScaleTable(1)*f0*t8);
% so0e = mod8.*cos(2*pi*ScaleTable(2)*f0*t8);
% la0e = mod8.*cos(2*pi*ScaleTable(3)*f0*t8);
% ti0e = mod8.*cos(2*pi*ScaleTable(4)*f0*t8);
% do1e = mod8.*cos(2*pi*ScaleTable(5)*f0*t8);
% re1e = mod8.*cos(2*pi*ScaleTable(6)*f0*t8);
% mi1e = mod8.*cos(2*pi*ScaleTable(7)*f0*t8);
% fa1e = mod8.*cos(2*pi*ScaleTable(8)*f0*t8);
% so1e = mod8.*cos(2*pi*ScaleTable(9)*f0*t8);
% la1e = mod8.*cos(2*pi*ScaleTable(10)*f0*t8);
% tb1e = mod8.*cos(2*pi*ScaleTable(11)*f0*t8);
% ti1e = mod8.*cos(2*pi*ScaleTable(12)*f0*t8);
% do2e = mod8.*cos(2*pi*ScaleTable(13)*f0*t8);
% re2e = mod8.*cos(2*pi*ScaleTable(14)*f0*t8);
% mi2e = mod8.*cos(2*pi*ScaleTable(15)*f0*t8);
% fa2e = mod8.*cos(2*pi*ScaleTable(16)*f0*t8);
% so2e = mod8.*cos(2*pi*ScaleTable(17)*f0*t8);
% la2e = mod8.*cos(2*pi*ScaleTable(18)*f0*t8);
% ti2e = mod8.*cos(2*pi*ScaleTable(19)*f0*t8);
% do3e = mod8.*cos(2*pi*ScaleTable(20)*f0*t8);
% blke = zeros(1,j);
% % 1/16 notes
% do0s = mod16.*cos(2*pi*ScaleTable(21)*f0*t16);
% re0s = mod16.*cos(2*pi*ScaleTable(22)*f0*t16);
% mi0s = mod16.*cos(2*pi*ScaleTable(23)*f0*t16);
% fa0s = mod16.*cos(2*pi*ScaleTable(1)*f0*t16);
% so0s = mod16.*cos(2*pi*ScaleTable(2)*f0*t16);
% la0s = mod16.*cos(2*pi*ScaleTable(3)*f0*t16);
% ti0s = mod16.*cos(2*pi*ScaleTable(4)*f0*t16);
% do1s = mod16.*cos(2*pi*ScaleTable(5)*f0*t16);
% re1s = mod16.*cos(2*pi*ScaleTable(6)*f0*t16);
% mi1s = mod16.*cos(2*pi*ScaleTable(7)*f0*t16);
% fa1s = mod16.*cos(2*pi*ScaleTable(8)*f0*t16);
% so1s = mod16.*cos(2*pi*ScaleTable(9)*f0*t16);
% la1s = mod16.*cos(2*pi*ScaleTable(10)*f0*t16);
% tb1s = mod16.*cos(2*pi*ScaleTable(11)*f0*t16);
% ti1s = mod16.*cos(2*pi*ScaleTable(12)*f0*t16);
% do2s = mod16.*cos(2*pi*ScaleTable(13)*f0*t16);
% re2s = mod16.*cos(2*pi*ScaleTable(14)*f0*t16);
% mi2s = mod16.*cos(2*pi*ScaleTable(15)*f0*t16);
% fa2s = mod16.*cos(2*pi*ScaleTable(16)*f0*t16);
% so2s = mod16.*cos(2*pi*ScaleTable(17)*f0*t16);
% la2s = mod16.*cos(2*pi*ScaleTable(18)*f0*t16);
% ti2s = mod16.*cos(2*pi*ScaleTable(19)*f0*t16);
% do3s = mod16.*cos(2*pi*ScaleTable(20)*f0*t16);
% blks = zeros(1,k);
% % Melody by Schau_mal
% part0 = [mi1f la0e la0e do1f mi1f ...
% re1e re1s mi1s re1e do1e re1e do1e la0f ...
% mi1f la0e la0e do1f mi1f ...
% so1e re1s mi1s re1e do1e re1e do1e ti0e so0e ...
% mi1f la0e la0e do1f mi1f ...
% re1e re1s mi1s re1e do1e re1e do1e la0e so0e ...
% mi1f la0e la0e do1f mi1f ...
% so1e mi1e blkf blkf blkf ...
% ];
% part1 = [la0f la0e so0e la0f la0e do1e ...
% do1f re1e do1e la0f la0f ...
% do1f do1e so0e do1e re1e mi1e so1e ...
% so1e mi1e re1f mi1f mi1f ...
% la1e la1e la1e so1e mi1e mi1f do1e ...
% la0e la0e la0e mi1e re1s mi1s re1e re1f ...
% mi1e mi1e so1e mi1e re1e mi1e re1e do1e ...
% la0f so0f la0f la0f ...
% ];
% part2 = [mi1e mi1e so1e mi1e mi1e so1e so1e la1e ...
% do2e la1e so1f la1s do2s la1e la1f ...
% la0f la0e so0e la0f do1f ...
% re1e mi1s re1s do1e re1e mi1f mi1f ...
% la0e la1e la1e so1e re1e mi1s re1s do1e re1e ...
% mi1f mi1f blke blke blkf ...
% do1e la0e la0e do1e re1f so0e so0e ...
% mi1e so1e mi1e re1e do1f do1f ...
% la0e do1e re1e mi1e re1e do1e so0e mi0e ...
% la0f la0f blke blke blkf ...
% ];
% part3 = [la0f la0e so0e la0f do1f ...
% re1e mi1s re1s do1e re1e mi1f mi1f ...
% la0e la1e la1e so1e re1e mi1s re1s do1e re1e ...
% mi1f mi1f blke blke blkf ...
% do1e la0e la0e do1e re1f so0e so0e ...
% mi1e so1e mi1e re1e do1f do1e do1e ...
% la0e do1e re1e mi1e so1e mi1e mi1e so1e ...
% la1f la1f la1f la1f ...
% ];
% part4 = [la1e la1s la1s la1e la1e la1e la1s so1s mi1e re1e ...
% re1e re1s re1s mi1e mi1s so1s mi1e mi1s re1s do1e do1s la0s ...
% la0f la0e so0e la0f la0e do1e ...
% re1e mi1s re1s do1e re1e mi1f mi1f ...
% la1e so1e mi1e re1e so1e mi1e re1e do1e ...
% do1f do1f la0s do1s re1s mi1s re1s do1s la0s do1s
% ];
% part5 = [do2e do2s do2s la1e la1s la1s so1e so1s so1s mi1e mi1s mi1s ...
% re1e mi1s re1s do1e la0s so0s la0s so0s do1s re1s mi1s so1s la1s re2s ...
% do2f do2f blks blks blks blks do1e re1e ...
% mi1f mi1f mi1f so1e mi1e ...
% la1f la1f la1e do1e so1e mi1e ...
% re1f re1e re1s re1s re1e re1e do1e re1e ...
% mi1f mi1e mi1s mi1s mi1e re1s do1s ti0e do1s re1s ...
% mi1f mi1f mi1f so1e mi1e ...
% do2f la1f la1f la1e do1e ...
% re1f so1f so1f la1f ...
% ti1f ti1f ti1f ti1f ...
% ];
% part6 = [blkf blkf mi1e so1e mi1e so1e ...
% mi1f la0e la0s la0s do1f la0e mi1s la0s ...
% do1e do1s do1s re1e do1s re1s mi1f mi1f ...
% mi1f la0e la0s la0s so1f re1e re1s re1s ...
% mi1f mi1f mi1s re1s do1s la0s mi0s re0s mi0s so0s ...
% do1f la0e la0s la0s re1f so0e so0s so0s ...
% mi0f so0e so0s so0s do1f do1f ...
% la0f do1e do1s la0s mi1e mi1s mi1s re1e re1s mi1s ...
% ];
% % Combination, v1 is complete version, v2 is simple version.
% v1 = [part0 part1 part1 part2 part3 part4 part0 part1 part1 part2 part3 part5 part3 part6 part3];
% v2 = [part0 part1 part1 part2 part3 part5 part3 part6 part3];
% % Let's rock ^_^
% s = v1;
% s = s/max(s);
% sound(s,fs);
% 五星红旗
function varargout=saxplaxliu(varargin)
x=[0 0 40 40 0];
y=[0 13