clc
clear all;
close all;

Fs=8000;
a=440;
q=2^(1/12);

c=(q^-9*a);
d=(q^-7*a);
e=(q^-5*a);
f=(q^-4*a);
g=(q^-2*a);
h=(q^2*a);
c2=(q^3*a);


bpm=60;
ct=60/bpm;
pul=2*ct;
os=1/2*ct;

tct=0:1/Fs:ct-1/Fs;
tpul=0:1/Fs:pul-1/Fs;
tos=0:1/Fs:os-1/Fs;

skokct = ones(1,8000);
krok=0;
for i=1:800
   krok=krok+1/800;
   skokct(i)=krok;
end
krok=1;
for i=7200:8000
   krok=krok-1/800;
   skokct(i)=krok;
end

skokos = ones(1,4000);
krok=0;
for i=1:400
   krok=krok+1/400;
   skokos(i)=krok;
end

krok=1;
for i=3600:4000
   krok=krok-1/400;
   skokos(i)=krok;
end

skokpul = ones(1,16000);
krok=0;
for i=1:1600
   krok=krok+1/1600;
   skokpul(i)=krok;
end

krok=1;
for i=14400:16000
   krok=krok-1/1600;
   skokpul(i)=krok;
end

Cos=cos(2*pi*c*tos).*skokos;
dos=cos(2*pi*d*tos).*skokos;
eos=cos(2*pi*e*tos).*skokos;
fos=cos(2*pi*f*tos).*skokos;
gos=cos(2*pi*g*tos).*skokos;
aos=cos(2*pi*a*tos).*skokos;
hos=cos(2*pi*h*tos).*skokos;
c2os=cos(2*pi*c2*tos).*skokos;

cpul=cos(2*pi*c*tpul).*skokpul;
dpul=cos(2*pi*d*tpul).*skokpul;
epul=cos(2*pi*e*tpul).*skokpul;
fpul=cos(2*pi*f*tpul).*skokpul;
gpul=cos(2*pi*g*tpul).*skokpul;
apul=cos(2*pi*a*tpul).*skokpul;
hpul=cos(2*pi*h*tpul).*skokpul;
c2pul=cos(2*pi*c2*tpul).*skokpul;

cct=cos(2*pi*c*tct).*skokct;
dct=cos(2*pi*d*tct).*skokct;
ect=cos(2*pi*e*tct).*skokct;
fct=cos(2*pi*f*tct).*skokct;
gct=cos(2*pi*g*tct).*skokct;
act=cos(2*pi*a*tct).*skokct;
hct=cos(2*pi*h*tct).*skokct;
c2ct=cos(2*pi*c2*tct).*skokct;

holubicka = [Cos dos eos fos gos aos hos gos c2ct c2ct hpul act act gpul fos aos fos dos eos gos eos Cos dct dct gpul fos aos fos dos eos gos eos Cos dct dct cpul];
sound(holubicka, Fs)