function  [FY]=INTEGR(N,F)
%
%	***  Integration Subroutine  ***
%
%       Real F(0:10,0:10), FY, D1     定義變量
%       Integer N, M
%
        
        SPLINE;            %求出擬合系數(shù)
        FY = 0.0;%FY = 0.0;                      輸出變量先置0
        for I = 1:N  %Do 35 I = 0,N-1  
            D1 = 0.0;       %D1 = 0.0
            for M = 5:-1:2  %Do 7 M = 4,1,-1
                D1 = (D1 + F(M,I)/(M-1))*F(6,I); %D1 = (D1 + F(M,I)/Float(M))*F(5,I)
            end      %7     Continue
            FY = FY + D1;%      FY = FY + D1
        end%35    Continue
%      Return
       return