function N=weierstrass(a,b)
% Den här matlabfunktionen beräknar funktionsserien 
% f(x) = summa 2^(-k/2)*sin((2^k)*x) (k löper från noll till oändligheten). 
% f(x) är kontinuerlig men ingenstans deriverbar!
tol=(b-a)/20;
N = ceil(2.*log((1+sqrt(2))./tol)./log(2)); % Här bestäms hur många termer
%% vi skall ta med i serien så att önskad noggrannhet erhålls.
x=linspace(a,b,999);xk=x;
r=1/sqrt(2);c=1;
y = sin(xk);% 0:e termen i funktionsserien.
for k=1:N;
	c=c*r;
   xk=rem(2*xk,2*pi); % Heltalsmultipler av 2 pi tas bort!
   y = y+c*sin(xk); % Term efter term läggs nu till föregående partialsumma.
end
figure(1);plot(x,y)
