// Bat dau tinh sai so
		for (int v = start; v <= end; v++)
			saiso = saiso + pow((dau[v][1] - a*dau[v][0] - b),2);
		saiso = sqrt(1.00 / (double)D * saiso);