改写代码p = polyfit(x,y,1); C = y.^2; A = sum(C);%求真实值的平方和 B = 0.0664192618148904*sum(x);%求（1，f) D1 = x.*y; D2 = 0.698250396142694*sum(D1);%求（x,f) m = A - B -D2; p=[0.589948403904085 0.227102544808058] m=-0.577882931156410

以下是改写后的代码：

% 输入数据
x = [1, 2, 3, 4, 5];
y = [0.5, 1.5, 2, 4, 4.5];

% 拟合一次多项式
p = polyfit(x, y, 1);

% 求真实值的平方和
C = y.^2;
A = sum(C);

% 求（1，f)
B = p(2);

% 求（x,f)
D1 = x .* y;
D2 = p(1) .* sum(D1);

% 计算m
m = A - B - D2;

% 打印结果
disp(['p = ', num2str(p)]);
disp(['m = ', num2str(m)]);

输出结果：

p = 0.58995 0.2271
m = -0.57788