scilab 5_3 (The least squares method)

clc;clear;
//Функция, вычисляющая разность между экспериментальными
     //и теоретическими значениями.
     //Перед использованием необходимо определить
     //z=[x;y] - матрицу исходных данных - и
     //с - вектор начальных значений коэффициентов,
     //размерность вектора должна совпадать
     //с количеством искомых коэффициентов.
     function [zr]=G(c,z)
     zr=z(2)-c(1)-c(2)*z(1)-c(3)*z(1)^2-c(4)*z(1)^3
     endfunction
     //Исходные данные
     x=[-0.3 0.1 0.5 0.9 1.3 1.7 2.1];
     y=[-1.3 2.2 3.5 5.8 7.2 6.0 5.8];
     //Формирование матрицы исходных данных
     z=[x;y];
     //Вектор начальных приближений
     c=[0;0;0;0];
     //Решение задачи
     [a,err]=datafit(G,z,c)
          //Построение графика экспериментальных данных
      plot2d(x,y,-4);
      //Построение графика подобранной функции
      t=-0.3:0.01:2.1;
      Ptc=a(1)+a(2)*t+a(3)*t^2+a(4)*t^3;
      plot2d(t,Ptc);