function integrals(a,b,n)
    h = (b-a)/n;
    score2 = 0;
    for i = a:b
        if i == a
            score = i^6 + 2*i^5 + 5*i^4 + 3*i^3 + 7*i^2 - i + 1;
        elseif i == b
            score1 = i^6 + 2*i^5 + 5*i^4 + 3*i^3 + 7*i^2 - i + 1;
        else
            score2 = score2 + i^6 + 2*i^5 + 5*i^4 + 3*i^3 + 7*i^2 - i + 1;
        end
    end
    finalScore = h * (((score+score1)/2) + score2);
    disp(finalScore);
end