function [t,z,nz] =konvydin(P)
%Tekee konvoluutioytimen käyttäjän ampumista P pisteestä.
%
%[t,z,nz] =konvydin(7)
h=figure; hold on, grid on
axis([-1,1,1,10^5])
%xlabel('x')
%ylabel('log(y)')
eval(['title(''Ammu ' num2str(P) ' pistettä'')']);
for k=1:P
g(k,:)=ginput(1);
plot(g(k,1),g(k,2),'k*')
end

x=g(:,1);y=g(:,2);
[x,j]=sort(x);y=y(j);
[x,j,~] = unique(x);y=y(j);
xmin=(-pi+1000*x(1))/(1001);
xmax=(pi+1000*x(end))/(1001);
x=[-pi,xmin,x',xmax,pi]';y=[0,0,y',0,0]';
t=linspace(-pi,pi,1000);
z=interp1(x,y,t);
I=mean(z)*2*pi;
nz=z/I;
plot(t,z)
plot(t,nz,'r')

axesObjs = get(h, 'Children');
dataObjs = get(axesObjs, 'Children');
legend([dataObjs(1),dataObjs(2)],'normeerattu','murtoviiva')
end

