function res = Ak(k)
% задаваме стойностите на х със съответната стъпка
x = 0:0.004:5 ; 
% записваме под интегралната функция у = x*sin^2(pi*x)cos(k*pi*x/4)
y = x.*sin(pi*x).*sin(pi*x).*cos(k*pi*x/4) ;
% решаваме интеграла чрез функцията trapz
res = trapz(x,y)*1/2 ;