2021-10-30

计算方法上机实践四 4.2

syms x y n;
n=5;
x=[1 2 3 4 5];
y=[4 4.5 6 8 8.5];
f1=sym(x.^2);
f2=sym(x.*y);
a = (sum(f1)*sum(y)-sum(x)*sum(f2))/(sum(f1)*n-(sum(x)^2))
b = (n*sum(f2)-sum(x)*sum(y))/(n*sum(f1)-sum(x)^2)
plot(x,a+b*x)