function ecFR(N,I)
%Ionescu Mihai
%Grupa 1341
%Pentru
N=2, I=25

a=0;
b=pi;
alpha=2;
h=(b-a)/I;

    function f=f(x)
        f=sin(alpha*x)+(pi*x/10*alpha)*(cos(alpha*x)-1);
    end
    function k=k(x,y)
        k=x*y/50;
    end
        for i=1:I+1
            x(i)=(i-1)*h;
        end


        for i=1:I+1
            U(1,i)=f(x(i));
        end
        plot(x,U(1,:))

        for n=1:N
            for i=1:I+1
                for j=1:I+1
                    val(j)=k(x(i),x(j))*U(n,j);
                end
                U(n+1,i)=trapz(x,val)+f(x(i));

                plot(x,U(n+1,:))

            end
        end



end