Лабы по МПО / Лабы / Лаб3 / СPP / MY_NAR_O

#include <stdlib.h>
#include <stdio.h>

#include "Sampler.h"
#define size 20

void main(void)
{
	int i,k;
	double a,b,x[size],y[size],y_calc[size];
	double sum_x,sum_y,sum_xy,sum_x2,sum_y2,xi,yi,sxy,sxx,syy;

for (k=0;k<80;k++)
{
SAMPLE;//15
	randomize();
SAMPLE;//17
	printf("Input Data: \n");
SAMPLE;//19
	for (i=0;i<size;i++)
	{
		x[i] = 10+i*10;
		y[i] = 10+i*10 + random(20) - 10;
	}
SAMPLE;//25
	sum_x = 0;
	sum_y = 0;
	sum_xy = 0;
	sum_x2 = 0;
	sum_y2 = 0;
	for (i=0; i<size; i++)
	{
		xi=x[i];
		yi=y[i];
		sum_x += xi;
		sum_y += yi;
		sum_xy += xi*yi;
		sum_x2 += xi*xi;
		sum_y2 += yi*yi;
	}
	sxx = sum_x2-sum_x*sum_x/size;
	sxy = sum_xy-sum_x*sum_y/size;
	syy = sum_y2-sum_y*sum_y/size;
	b=sxy/sxx;
	a=((sum_x2*sum_y-sum_x*sum_xy)/size)/sxx;
	for (i=0;i<size;i++)
		y_calc[i] = a + b * x[i];
SAMPLE;//48
	printf("a %3.1f,b %3.1f\n",a,b);
SAMPLE;//50
	for (i=0;i<size;i++)
		printf("y %3.1f\n",y_calc[i]);
SAMPLE;//53
}
}