пример от ръководството:

Clear[A,b]
A=(({{5, 10}, {10, 30}}));b={8,20};
LinearSolve[A,b]

Clear[A,b]
A=(({{5, 10, 30}, {10, 30, 100}, {30, 100, 354}}));b={8,20, 70};
LinearSolve[A,b]

Пример2:

табулиране на функция:

Clear[xx,yy,f,x,a,b];
a=1.;b=2;
n=11; (*брой точки на таблицата*)
h=((b-a)/(n-1));
f[x_]:=x+Sin[(x^2)]
xx=Table[a+i*h,{i,0,n-1}]
yy=f[xx]

Table[(i^2),{i,1,5}]

търсим полином от първа степен  P1=a1*x + a0 => търсим a0 и a1:

броя на точките в кода е означен с n, а във формулата е с N

Clear[A,b,a]
A=(({{n, (∑_(i=1)^n)xx[[i]]}, {(∑_(i=1)^n)xx[[i]], (∑_(i=1)^n)((xx[[i]])^2)}}));
MatrixForm[A]
b={(∑_(i=1)^n)yy[[i]], (∑_(i=1)^n)xx[[i]]*yy[[i]]}
a=LinearSolve[A,b]
P1[x_]:=a[[2]]x+a[[1]]
P1[x]

Извод: Полиномът от първа степен е  P1*(x)=-0.697508x + 2.99685

графично представяне на резултатите:

построяваме табулираната функция (получените точки)

points=Table[{xx[[i]],yy[[i]]},{i,1,11}];
grpoints=ListPlot[points]

графика на полинома на линейна регресия:

grP1=Plot[P1[x],{x,1,2}, PlotStyle->{Black}]

Show[grpoints,grP1]

оценка на грешката:

errP1=(Sqrt[(∑_(i=1)^n)(((yy[[i]]-P1[xx[[i]]]))^2)])

P1[1.05]