最小二乘算法 C 语言实现

//
//  main.m
//  test
//
//  Created by Jack on 15/11/20.
//  Copyright © 2015年 宇之楓鷙. All rights reserved.
//

#import <stdio.h>
#import <stdlib.h>
#import <malloc/malloc.h>
#import <math.h>


Smooth(double *x,double *y,double *a,int n,int m,double *dt1,double *dt2,double *dt3)//(x,y,a,n,m,dt1,dt2,dt3 )
//double *x; /*实型一维数组，输入参数，存放节点的xi值*/
//double *y; /*实型一维数组，输入参数，存放节点的yi值*/
//double *a; /*双精度实型一维数组，长度为m。返回m一1次拟合多项式的m个系数*/
//int n; /*整型变量，输入参数，给定数据点的个数*/
//int m; /*整型变量，输入参数，拟合多项式的项数*/
//double *dt1; /*实型变量，输出参数，拟合多项式与数据点偏差的平方和*/
//double *dt2; /*实型变量，输出参数，拟合多项式与数据点偏差的绝对值之和*/
//double *dt3; /*实型变量，输出参数，拟合多项式与数据点偏差的绝对值最大值*/
{
    int i ,j ,k ;
    double *s,*t,*b,z,d1,p,c,d2,g,q,dt;
    /*分别为s ,t ,b分配存贮空间*/
    s = (double *)calloc(n,sizeof(double));
    if(s == NULL)
    {
        printf("内存分配失败\n");