Вычислите значение функции y=x^3-x^2+16x-43 для x изменяющегося в диапазоне от-4 до 4 включительно с шагом 0.5. нужно написать на паскале

Program task;
var i,a:real;
begin     
i:=-4;     
while (i<=4) do begin           
a:=(i*i*i)-(i*i)+16*i-43;           
writeln(a:0:2);           
i+=0.5;     
end;     
readln;
end.