testgjls.cpp

来自「eC++编译器源码」· C++ 代码 · 共 26 行

#include <GaussJ.h>
#include <LeastSquare.h>
#include <iostream.h>
#include <iomanip.h>

void main()
{PolyReturn pr; PolyArray pa; DataPoints d; LineReturn lr;
  //GaussJordan(pr, pa, 0, 0);
  d.noNums = 4;
  d.X[1] = 2.0;  d.Y[1] = 3.0;
  d.X[2] = 5.0;  d.Y[2] = 7.5;
  d.X[3] = 7.0;  d.Y[3] = 10.5;
  d.X[4] = 12.0;  d.Y[4] = 18.0;
  LineFit(d,0.0, lr);
  cout<<setw(10)<<lr.A<<lr.B<<lr.sigA<<lr.sigB<<lr.chi2<<lr.Q;
  d.X[1] = 2.0;  d.Y[1] = 5.0;
  d.X[2] = 5.0;  d.Y[2] = 26.0;
  d.X[3] = 7.0;  d.Y[3] = 50.0;
  d.X[4] = 12.0;  d.Y[4] = 145.0;
  pr.minFit = 3;
  pr.noTerms = 3;
  pr.listA[1]=1; pr.listA[2]=2; pr.listA[3]=3;
  PolyFit(d,pr);
  cout<<pr.minFit<<pr.coeff[1]<<pr.coeff[2];
};