function  summation
maxN = input('Enter the maximum value of N required: ');
I(1)=f(1);
for N=2:maxN
    I(N)=I(N-1)+f(N);
end

disp(['Values of I_N'])
disp([1:N;I])