I think there is something wrong but I cannot find it

H=[0.40229 0.97419 1.4829 2.2318 2.7139 3.336 3.9664];

I=[8.2823 43.043 98.033 219.98 323.67 487.7 687.84];

for n=1:5

pp = polyfit(H,I,n);

f = polyval(pp,H);

subplot(3,2,n);

plot(H,I,'.',H,f)

end

The following data represents the volume of ice needed to chill a given volume of hot liquid. (Use HotLiquid.txt and lce.txt from the DataFiles.zip folder. See instructions above for help.)+ Volume of Hot Liquid (cups)0.40 0.97 Volume of lce (cups)^ 1.48 3.34 3.97 8.2843.04+ 98.03- ...487.7687.84s^t a) Fit a 1st, 2nd, 3rd, and 4th degree polynomial to the dataset described by volume of liquid vs. volume of ice. Use subplot to draw all four. Plot your dataset as points, and your fitted equation as a solid line in the domain 1 s Vwater s 5. Use a legend to identify the data points and the fitted equation. + b) Which degree do you think best represents the data? Why? (Hint: how much ice is needed to chill 0 cups of hot liquid?) Answer in the appropriate box.* c) Use your chosen degree polynomial to find and print the volume of ice needed to chill 1.5 cups of hot liquid. 4 Self-Check:^ Part a). *- i500 Part c): Answer will vary depending on the degree chosen, but you can check the graph to make sure the answer is in the right ballpark.*