Y1 = peso_laranja
Y2 = peso_toranja
X1 = diametro_laranja
X2 = diametro_toranja
n1 = np.size(peso_laranja)
n2 = np.size(peso_toranja)
a1 = (n1*np.sum(X1*Y1) - np.sum(X1)*np.sum(Y1))/(n1*np.sum(X1**2)-(np.sum(X1)**2))
a2 = (n2*np.sum(X2*Y2) - np.sum(X2)*np.sum(Y2))/(n2*np.sum(X2**2)-(np.sum(X2)**2))
b1 = np.mean(Y2) - a1*np.mean(X1)
b2 = np.mean(Y2) - a2*np.mean(X2)
y1 = a1*X1+b1
y2 = a2*X2+b2
plt.plot(diametro_laranja,peso_laranja)
plt.plot(X1,y1)
plt.plot(diametro_toranja, peso_toranja)
plt.plot(X2,y2)
plt.legend(['Peso Laranja', 'Crescimento esperado Laranja', 'Peso Toranja', 'Crescimento esperado Toranja'])