xbasc()
N=64;
h=1/N;
x=0:h:1;
//--------------------------------------------------------

...

//--------------------------------------------------------
[L,I]=sort(a./b);              // scilab sortiert abwaerts
approx_EW=L(N-k)
exakter_EW=(k*%pi)^2
plot2d(x,sin(k*%pi*x),1)
plot2d(x,[0;X(:,I(N-k));0],2)
