clc;
clear all;
close all;
fs=5000;
dt=1/fs;
N=1*fs;
n=0:N-1;
t=dt*n;

%exp(-2*t) funkcja gasnaca
%sin funkcja oscylacyjna
do = exp(-2*t).*sin(2*pi*261.6*t);
re = exp(-2*t).*sin(2*pi*293.7*t);
mi = exp(-2*t).*sin(2*pi*329.6*t);
fa = exp(-2*t).*sin(2*pi*349.6*t);
sol = exp(-2*t).*sin(2*pi*391.9*t);
la = exp(-2*t).*sin(2*pi*440*t);
si = exp(-2*t).*sin(2*pi*493.9*t);

do = exp(-5*t).*sin(2*pi*261.6*t);
x = [do re mi fa sol la si la sol fa mi re do];
soundsc(x)
plot(x)