Compartilho meu código:
1.Em partes:
#Toranjas
Xt = diametro_toranja
Yt = peso_toranja
nt = np.size(Xt)
at = (nt*np.sum(Xt*Yt) - np.sum(Xt)*np.sum(Yt))/(nt*np.sum(Xt**2)-np.sum(Xt)**2)
bt = np.mean(Yt) - at*np.mean(Xt)
xt = diametro_toranja
yt = at*Xt + bt
plt.plot (diametro_toranja, peso_toranja, color='green')
plt.plot(xt,yt, color='red')
#Laranjas
Xl = diametro_laranja
Yl = peso_laranja
nl = np.size(Xl)
al = (nl*np.sum(Xl*Yl) - np.sum(Xl)*np.sum(Yl))/(nl*np.sum(Xl**2)-np.sum(Xl)**2)
bl = np.mean(Yl) - al*np.mean(Xl)
xl = diametro_laranja
yl = al*Xl + bl
plt.plot (diametro_laranja, peso_laranja, color='orange')
plt.plot(xl,yl, color='b')
2.Todas as curvas em um único gráfico:
plt.plot (diametro_laranja, peso_laranja, color='orange')
plt.plot (diametro_toranja, peso_toranja, color='green')
plt.plot(xl,yl, color='b')
plt.plot(xt,yt, color='red')
plt.legend(['Laranja','Toranja'])
plt.xlabel('Diâmetro (cm)')
plt.ylabel('Peso (g)')
plt.title('Diâmetro x Peso das Toranjas e Laranjas')
Duvidas: nesta versão com todas as curvas na mesma área, como inserir as legendas?