A resposta do @Marcos Banikresposta do @Marcos Banik pode ser traduzida por exemplo para Python com um mínimo esforço:
>>> x = [1,2,3,4]
>>> y = [10,20,30,40]
>>> m = sum(a*b for (a,b) in zip(x,y)) - sum(x) * sum(y) / len(x)
>>> m /= sum(a**2 for a in x) - (sum(x)**2)/len(x)
>>> b = sum(y)/len(y) - m * sum(x)/len(x)
>>> m*5 + b
50
Entretanto nem toda linguagem imperativa possui o mesmo nível de expressividade, por isso vou apresentar uma versão mais elaborada (em JavaScript):
var x = [1,2,3,4];
var y = [10,20,30,40];
function produto(x, y) {
var ret = [];
for ( var i = 0 ; i < x.length ; i++ )
ret.push(x[i] * y[i]);
return ret;
}
function quadrados(x) {
var ret = [];
for ( var i = 0 ; i < x.length ; i++ )
ret.push(x[i] * x[i]);
return ret;
}
function somatorio(x) {
var ret = 0;
for ( var i = 0 ; i < x.length ; i++ )
ret += x[i];
return ret;
}
function media(x) { return somatorio(x) / x.length; }
var m = somatorio(produto(x,y)) - somatorio(x) * somatorio(y) / x.length;
m /= somatorio(quadrados(x)) - somatorio(x)*somatorio(x) / x.length;
var b = media(y) - m * media(x);
console.log(m*5 + b);