%%% KOD GŁÓWNY WYWOŁYWANY W WIERSZU POLECEŃ POWINNO WSZYSTKO DZIAŁAĆ PO ZAŁADOWANIU MPLIKÓW TZN ZAPISANIU NA DYSKU
T=1;
t=-T:0.01:2*T;
w0=2*pi/T;
m=30;
k=1:30;
%F1=((((-1).^k)- 1)./((pi.^2)*(k.^2)));
%F2=(2*sin((pi*k))./(2) + 0.5*pi.*k.*(1 - ((-1).^k))./((pi^2)*(k.^2)));
F= ((((-1).^k)- 1)./((pi.^2)*(k.^2)))-j*(2*sin((pi*k))./(2) + 0.5*pi.*k.*(1 - ((-1).^k))./((pi^2)*(k.^2)));

F0=0;
Fs=F0;
for k=1:30
    Fk=F(k);
    Fs=Fs+2*abs(Fk)*cos(w0*t*k+angle(Fk));
end

subplot(2,1,1);
plot(t,Fs,'linewidth',2);
grid on;
xlabel('t [ms]');
ylabel('f(t)');
title('funkcja nr. 17 z zad 1');

F0 = 0;
for k=1:30
    Fk = F;
end

F(2) = -1i/4;

Fs = F0;
fpr = 0;
for k=1:30
    Fs = Fs + 2*abs(F(k))*cos(w0*k*t + angle(F(k)));
    fpr = fpr + abs(F(k))^2;
end
hold on;
plot(t, Fs);

%G = [F0 2*F];
%I=0:k;
%M=abs(G);
%fi = angle(G) * 180/pi;
%x=[I' M' fi']


hold off;


Fsk=sqrt(quad('f17r2',0,1))
fsk_parserv = sqrt(F0*2 + 2*fpr)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%Mplik 1; 'mf1.m'
function y=mf1(t)
y=(4*t).^2

%%% Mplik2; 'mf2.m'
function y=mf2(t)
y=1^2;

%%% Mplik3;'mf3.m'

function y=mf3(t);
y=(-4*(t-1/3)).^2;

%%% Mplik4;'mf4.m'

fuction y= mf4(t)
y=-1^2;

%%% Mplik5 'calka.m'

Fsk=(quad(mf1,0,1/4)+quad(mf2,1/4,1/2) + quad(mf3,1/2,3/4)+quad(mf4,3/4,1))

%%% Mplik6'f17r2.m'

function y= f17r2(t);
y=f17r(t).*f17r(t);

%%%Mplik7'f17r.m'

function y=f17r(t);

y=4*t.*rectpuls(t-0,1/4)+1.*rectpuls(t-1/4,1/2)-4*(t-1/2).*rectpuls(t-1/2,3/4)-1.*rectpuls(t-3/4,1);