blob: c21e12a770858c14f23020f1971a448af2b4b935 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
//Example 5.8
clc
clear
x = 1:4;
y = [1.65 2.7 4.5 7.35];
X = x;
Y = log10(y);
n = length(Y);
M1 = [sum(X.^2) sum(X); sum(X) n];
M2 = [sum(X.*Y); sum(Y)];
M = M1\M2;
B = M(1);
A = M(2);
a = 10^A;
b = B/log10(%e);
disp(round(a), "a =")
disp(round(b*10^4)/10^4, "b =")
|
