step = 0.1;
N = 2551;
x = 0;
y = 0;

for i = 1:N
    y = y + 0.5*x + (x*x)*cos(floor(x/4)-32);
    x = x + step;
end

fprintf('%f \n', y)