Regresie Trigo

#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#include <math.h>
#define NRMAX   10

#define pi 3.141592

void regtrig (int m, float x[NRMAX], float y[NRMAX], float omega, float *a, float *b)
{
    float sy,scos2x,scosx,sycosx;
    int i;

    sy=0;   scos2x=0;   scosx=0;    sycosx=0;
    for (i=0;i<m;i++)
        {
            sy=sy+y[i];
            scos2x=scos2x+(cos(omega*x[i])*cos(omega*x[i]));
            scosx=scosx+cos(omega*x[i]);
            sycosx=sycosx+(y[i]*cos(omega*x[i]));
        }
        *a=(sy*scos2x-scosx*sycosx)/(m*scos2x-scosx*scosx);
        *b=(m*sycosx-scosx*sy)/(m*scos2x-scosx*scosx);
}

int main (void)
{
    int i,m;
    float a,b,omega;
    float   X[NRMAX] = {0,pi/2, pi, 3*pi/2, 2*pi},
            Y[NRMAX] = {2, 1, 0, 1, 2};

    printf("REGRESIA TRIGONOMETRICA");
    printf("\nIntroduceti numarul experientelor:");
    scanf("%d",&m);
    printf("omega=");
    scanf("%f",&omega);
/*
    printf("\nIntroduceti abscisele si ordonatele functiei:");
    for (i=0;i<m;i++)
        {
            printf("\nX[%d]=",i);
            scanf("%f",&X[i]);
            printf("Y[%d]=",i);
            scanf("%f",&Y[i]);
        }
*/
    regtrig(m,X,Y,omega,&a,&b);
    printf("Rezultatul regresiei trigonometrice este: a=%f; b=%f",a,b);

    getch();
    return 1;
}