clear all; %close all;

j1 = 1:15;
j2 = 1:14;
N1 = (length(j1)-1)./2;
N2 = length(j2)./2;
nv1 = -1*N1:-1*N1;
nv2 = -1*N2:-1*N2;
x1 = j1*(2*pi)./length(j1);
x2 = j2*(2*pi)./length(j2);

%phi1 = 2./(sqrt(3*pi./4)*pi.^(1/4)).*(1-(x1-pi).^2./((pi./4)^2)).*exp(-(x1-pi).^2./(2*(pi./4).^2));
%phi2 = 2./(sqrt(3*pi./4)*pi.^(1/4)).*(1-(x2-pi).^2./((pi./4)^2)).*exp(-(x2-pi).^2./(2*(pi./4).^2));
phi1 = -(x1-pi).^2 + pi.^2 - sin(2.*x1)-12*exp(-(x1-4).^2);
phi2 = -(x2-pi).^2 + pi.^2 - sin(2.*x2)-12*exp(-(x2-4).^2);
%phi1 = square(2*x1-1);
%phi2 = square(2*x2-1);

for ind = 1:length(nv1)
    a1(ind) = 1./length(j1).*sum(phi1.*exp(-i*nv1(ind).*x1));
end
for ind = 1:length(nv2)
    a2(ind) = 1./length(j2).*sum(phi2.*exp(-i*nv2(ind).*x2));
end
for ind = 1:length(j1)
    phir1(ind) = sum(a1.*exp(i.*nv1.*x1(ind)));
end
for ind = 1:length(j2)
    phir2(ind) = sum(a2.*exp(i*nv2.*x2(ind)));
end

figure, 
subplot(2,1,1)
plot([x1], [phi1], 'r'), hold on, plot([0 x1], [NaN 5*phir1],'b');
subplot(2,1,2)
plot([x2], [phi2], 'r'), hold on, plot([0 x2], [NaN 5*phir2],'b')