% Name:Basel.m
clear
clf;
%========================================================================
%L=input('Type in the length N=');
x=-10:.01:100;
%========================================================================

Ib=1;m=1;
for k=1:25
    m=m*k;
    Ib=Ib+(((x/2).^k)/m).^2;
end

%=======================================================================

plot(x,Ib),
title('The Bssel function'),
xlabel('Index x')