gistfile1.txt

# -*- coding: utf-8 -*-
from pylab import *

N = 100
start = 0
end = 1

A = rand()
B = rand()

x = linspace(start,end,N)
y = A*x + B
y += randn(N)/10

p = polyfit(x,y,1)

figure()

title('Regressao linear')
plot(x,y,'o',label='A=%.2f, B=%.2f' % (A,B))
plot(x,polyval(p,x),'-',label='A=%.2f, B=%.2f' % tuple(p))
legend(loc='best')

show()