Programa pra Calcular Coeficiente Angular e Linear de Uma Determinada Reta

#include <iostream>

#include <cmath>

#include <vector>

using namespace std;

int main (){

setlocale(LC_ALL, "Portuguese");

double quadrado=0,soma_xi_quad=0, soma_yi_quad=0, soma_yi=0, soma_xi_yi=0, soma_xi=0, soma_xi_quadrado=0, soma_yi_quadrado=0, soma_sigma=0;

double a=0, b=0, d=0, sigma=0, sigma_p1=0, r=0, x_i, y_i, d_parte_1=0;

double a_parte_1=0, a_parte_2=0, a_parte_3=0;

double b_parte_1=0, b_parte_2=0, b_parte_3=0, b_parte_4=0;

double r_p1=0, r_p2=0,r_p3=0,r_p4=0;

double delta_a=0, delta_a_p1=0, delta_a_p2=0;

double delta_b=0, delta_b_p1=0, delta_b_p2=0, r_p=0,quadrado_p1=0;

int n, i;

cout<<"Este é um programa que calcula o coeficiente angular (A) e linear (B) de uma determinada curva."<<endl;

cout<<"Além dos parâmetros delta A, delta B e R^2."<<endl;

do{

cout<<"Insira o número de pontos presentes em sua curva:"<<endl;

cin>>n;

if(n<=2){

cout<<"O numero de pontos presentes na curva deve ser maior 2!!!"<<endl;

}

}while(n<=2);

vector<double>vec_xi;

vec_xi.reserve(n);

vector<double>vec_yi;

vec_yi.reserve(n);

for (i=0; i<n; i++){

cout<<"Insira o valor de X "<<i<<":"<<endl;

cin>>x_i;

vec_xi.push_back(x_i);

}

cout<<"---------------Atenção, agora você irá preencher os valores de Y do seu gráfico-------------------"<<endl<<endl;

for (i=0; i<n; i++){

cout<<"Insira o valor de Y "<<i<<":"<<endl;

cin>>y_i;

vec_yi.push_back(y_i);

}

for(i=0;i<n;i++){

soma_xi_quad = soma_xi_quad + (pow(vec_xi[i],2));

soma_xi = soma_xi + vec_xi[i];

soma_xi_yi = soma_xi_yi + (vec_xi[i]*vec_yi[i]);

soma_yi = soma_yi + vec_yi[i];

soma_yi_quad = soma_yi_quad + (pow(vec_yi[i],2));

}

soma_xi_quadrado = pow(soma_xi,2);

soma_yi_quadrado = pow(soma_yi,2);

d_parte_1 = (n*soma_xi_quad);

d = (d_parte_1- soma_xi_quadrado);

b_parte_1 = (soma_xi_quad*soma_yi);

b_parte_2 = (soma_xi*soma_xi_yi);

b_parte_3 = (b_parte_1-b_parte_2 );

b_parte_4