#Código para ajuste de reta utilizando Octave Online https://octave-online.net/#

pkg load optim 

x=[10.2;15;21;24.5;30;34.7;40.1;46;49.7;56];# Substitua com os valores de x 
y=[1.0;2.0;3.0;4.0;5.0;6.0;7.0;8.0;9.0;10.0];# Substitua com os valores de y 
sx=0.01*x # Substitua com o cálculo da incerteza de x ou diretamente com os valores
sy=0.01*y # Substitua com o cálculo da incerteza de y ou diretamente com os valores
F = [ones(size(x)(1), 1), x(:)];
[p, e_var, r, p_var, fit_var] = LinearRegression (F, y);
parameters_StdErr = [p, sqrt(p_var)]
b=parameters_StdErr(1,1);
sb=parameters_StdErr(1,2);
a=parameters_StdErr(2,1);
sa=parameters_StdErr(2,2);
leg=strcat("Y = (",num2str(a,2),"+-",num2str(sa,2),") X + (",num2str(b,2),"+-",num2str(sb,2),")")

#### Gráfico - Edite o título e os eixos ####
yFit = F * p;
g1=errorbar(x, y, sx,sx,sy,sy,'~>');
set(g1, 'LineStyle','none');
hold on;
plot(x, yFit, '--r')
ylabel('y(unidade)');
xlabel('x(unidade)');
title('Título');
legend ('Dados',leg,'location','northwest');