function [y,N,tol]=Wstrs(x,tol)
a=min(x);b=max(x);l=length(x);

if ~exist('tol')
   if a==b
      tol = 1E-4;
   else
      tol = (b-a)./(l.*100);
   end
end
xk=rem(2*x,2*pi);

%Nedan beräknas antalet termer (N) som skall tas med i serien 
%för att felet skall bli mindre än  tol.
N = ceil(2.*log((1+sqrt(2))./tol)./log(2));

y = zeros(size(x));
y= sin(xk)./sqrt(2);

for k=2:N;
   xk=2.*xk;
   xk=rem(xk,2*pi); 
   y = y+sin(xk)./sqrt(2.^k); 
end

