Vericidad de las resistencias

En este programa se ingresan los valores de una resistencia o varias resistencias y de aquí el programa calcula un promedio y rango, tanto individual como general y posteriormente imprime los resultados en listas y al final en dos graficas.
En esta ocasión se presentan dos formas de hacer el programa, una en la cual solo se usan listas y matrices y en el siguiente el ciclo for con matrices.

Primer programa
Usando matrices y listas

from pylab import *
print("1º resistencia")
n1a=int(input("Dime la primer medicion de la resistencia: "))
n1b=int(input("Dime la segunda medicion de la resistencia: "))
n1c=int(input("Dime la tercera medicion de la resistencia: "))
p1=((n1a+n1b+n1c)/3)
a1=np.array([n1a, n1b, n1c])
r1=((a1.max())-(a1.min()))
print()
print("2º resistencia")
n2a=int(input("Dime la primer medicion de la resistencia: "))
n2b=int(input("Dime la segunda medicion de la resistencia: "))
n2c=int(input("Dime la tercera medicion de la resistencia: "))
p2=((n2a+n2b+n2c)/3)
a2=np.array([n2a, n2b, n2c])
r2=((a2.max())-(a2.min()))
print()
print("3º resistencia")
n3a=int(input("Dime la primer medicion de la resistencia: "))
n3b=int(input("Dime la segunda medicion de la resistencia: "))
n3c=int(input("Dime la tercera medicion de la resistencia: "))
p3=((n3a+n3b+n3c)/3)
a3=np.array([n3a, n3b, n3c])
r3=((a3.max())-(a3.min()))
print()
print("4º resistencia")
n4a=int(input("Dime la primer medicion de la resistencia: "))
n4b=int(input("Dime la segunda medicion de la resistencia: "))
n4c=int(input("Dime la tercera medicion de la resistencia: "))
p4=((n4a+n4b+n4c)/3)
a4=np.array([n4a, n4b, n4c])
r4=((a4.max())-(a4.min()))
print()
print("5º resistencia")
n5a=int(input("Dime la primer medicion de la resistencia: "))
n5b=int(input("Dime la segunda medicion de la resistencia: "))
n5c=int(input("Dime la tercera medicion de la resistencia: "))
p5=((n5a+n5b+n5c)/3)
a5=np.array([n5a, n5b, n5c])
r5=((a5.max())-(a5.min()))
print()
print("6º resistencia")
n6a=int(input("Dime la primer medicion de la resistencia: "))
n6b=int(input("Dime la segunda medicion de la resistencia: "))
n6c=int(input("Dime la tercera medicion de la resistencia: "))
p6=((n6a+n6b+n6c)/3)
a6=np.array([n6a, n6b, n6c])
r6=((a6.max())-(a6.min()))
print()
print("7º resistencia")
n7a=int(input("Dime la primer medicion de la resistencia: "))
n7b=int(input("Dime la segunda medicion de la resistencia: "))
n7c=int(input("Dime la tercera medicion de la resistencia: "))
p7=((n7a+n7b+n7c)/3)
a7=np.array([n7a, n7b, n7c])
r7=((a7.max())-(a7.min()))
print()
print("8º resistencia")
n8a=int(input("Dime la primer medicion de la resistencia: "))
n8b=int(input("Dime la segunda medicion de la resistencia: "))
n8c=int(input("Dime la tercera medicion de la resistencia: "))
p8=((n8a+n8b+n8c)/3)
a8=np.array([n8a, n8b, n8c])
r8=((a8.max())-(a8.min()))
print()
print("9º resistencia")
n9a=int(input("Dime la primer medicion de la resistencia: "))
n9b=int(input("Dime la segunda medicion de la resistencia: "))
n9c=int(input("Dime la tercera medicion de la resistencia: "))
p9=((n9a+n9b+n9c)/3)
a9=np.array([n9a, n9b, n9c])
r9=((a9.max())-(a9.min()))
print()
print("10º resistencia")
n10a=int(input("Dime la primer medicion de la resistencia: "))
n10b=int(input("Dime la segunda medicion de la resistencia: "))
n10c=int(input("Dime la tercera medicion de la resistencia: "))
p10=((n10a+n10b+n10c)/3)
a10=np.array([n10a, n10b, n10c])
r10=((a10.max())-(a10.min()))
print()
print("11º resistencia")
n11a=int(input("Dime la primer medicion de la resistencia: "))
n11b=int(input("Dime la segunda medicion de la resistencia: "))
n11c=int(input("Dime la tercera medicion de la resistencia: "))
p11=((n11a+n11b+n11c)/3)
a11=np.array([n11a, n11b, n11c])
r11=((a11.max())-(a11.min()))
print()
print("12º resistencia")
n12a=int(input("Dime la primer medicion de la resistencia: "))
n12b=int(input("Dime la segunda medicion de la resistencia: "))
n12c=int(input("Dime la tercera medicion de la resistencia: "))
p12=((n12a+n12b+n12c)/3)
a12=np.array([n12a, n12b, n12c])
r12=((a12.max())-(a12.min()))
print()
print("Mediciones de las resistencias")
print(a1)
print(a2)
print(a3)
print(a4)
print(a5)
print(a6)
print(a7)
print(a8)
print(a9)
print(a10)
print(a11)
print(a12)
print()
print("PROM\t RANG")
o1="{}\t {}\t".format(p1, r1)
o2="{}\t {}\t".format(p2, r2)
o3="{}\t {}\t".format(p3, r3)
o4="{}\t {}\t".format(p4, r4)
o5="{}\t {}\t".format(p5, r5)
o6="{}\t {}\t".format(p6, r6)
o7="{}\t {}\t".format(p7, r7)
o8="{}\t {}\t".format(p8, r8)
o9="{}\t {}\t".format(p9, r9)
o10="{}\t {}\t".format(p10, r10)
o11="{}\t {}\t".format(p11, r11)
o12="{}\t {}\t".format(p12, r12)
print(o1)
print(o2)
print(o3)
print(o4)
print(o5)
print(o6)
print(o7)
print(o8)
print(o9)
print(o10)
print(o11)
print(o12)
a2=[0, 0, 1.88, 1.023, 0.729, 0.577, 0.483, 0.419, 0.373, 0.377, 0.308, 0.285, 0.266, 0.249, 0.235, 0.223]
d3=[0, 0, 0, 0, 0, 0, 0, 0.076, 0.136, 0.184, 0.223, 0.256, 0.284, 0.308, 0.329, 0.348]
d4=[0, 0, 3.268, 2.574, 2.282, 2.114, 2.004, 1.924, 1.864, 1.816, 1.777, 1.744, 1.717, 1.692, 1.671, 1.652]
pp=((p1+p2+p3+p4+p5+p6+p7+p8+p9+p10+p11+p12)/12)
pr=((r1+r2+r3+r4+r5+r6+r7+r8+r9+r10+r11+r12)/12)
lcip=(pp-(a2[3]*pr))
lcir=(d3[3]*pr)
lcsp=(pp+(a2[3]*pr))
lcsr=(d4[3]*pr)
print("\tPROM\t RANG")
o16="{}\t {}\t {}\t".format("PROM",pp,pr)
o17="{}\t {}\t {}\t".format("LCS",lcsp,lcsr)
o18="{}\t {}\t {}\t".format("LCI",lcip,lcir)
print(o16)
print(o17)
print(o18)
l=(p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12)
lpp=(pp,pp,pp,pp,pp,pp,pp,pp,pp,pp,pp,pp)
llcs=(lcsp,lcsp,lcsp,lcsp,lcsp,lcsp,lcsp,lcsp,lcsp,lcsp,lcsp,lcsp)
llci=(lcip,lcip,lcip,lcip,lcip,lcip,lcip,lcip,lcip,lcip,lcip,lcip)
plt.plot(l,label="Resistencias", marker="d")
plt.plot(lpp, label="PROM", marker="s")
plt.plot(llcs, label="LCS", marker="x")
plt.plot(llci, label="LCI", marker="h")
plt.legend(loc="upper center")
plt.show()
l1=(r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,r12)
l1pr=(pr,pr,pr,pr,pr,pr,pr,pr,pr,pr,pr,pr)
l1lcs=(lcsr,lcsr,lcsr,lcsr,lcsr,lcsr,lcsr,lcsr,lcsr,lcsr,lcsr,lcsr)
l1lci=(lcir,lcir,lcir,lcir,lcir,lcir,lcir,lcir,lcir,lcir,lcir,lcir)
plt.plot(l1, label="R", marker="d")
plt.plot(l1pr, label="PROM", marker="s")
plt.plot(l1lcs, label="LCS", marker="x")
plt.plot(l1lci, label="LCI", marker="h")
plt.legend(loc="upper center")
plt.show()

Segundo programa
En este programa se pide el numero de resistencias que se desea medir y el numero de veces que se desea promediar

from pylab import *
import numpy as np
r=int(input("Dime el numero de resistencias: "))
c=int(input("Dime cuantas veces quieres medirla: "))
r0=0
c0=0
a0=[]
for r0 in range (r):
for c0 in range (c):
a=int(input("Dime la muestra %d posicion %d: " %((r0+1),(c0+1))))
a0.append(a)
x=np.array(a0)
x=x.reshape(r,c)
print(x)
maximo=[]
minimo=[]
ra=[]
for y in range(r):
maxim=np.max(x[y])
maximo.append(maxim)
print("Maximos", maximo)
for z in range(r):
minim=np.min(x[z])
minimo.append(minim)
print("Minimos", minimo)
for g in range(r):
ran=((maximo[g])-(minimo[g]))
ra.append(ran)
rango=np.sum(ra)
rango1=(rango/r)
print("Rangos",ra)
print("Rango total", rango1)
promedio=[]
for d in range(r):
y1=np.sum(x[ran])
y11=y1/c
promedio.append(y11)
print("Promedios",promedio)
prom=np.sum(promedio)
prom1=(prom/r)
print("Promedio total: ", prom1)
a2=[0, 0, 1.88, 1.023, 0.729, 0.577, 0.483, 0.419, 0.373, 0.377, 0.308, 0.285, 0.266, 0.249, 0.235, 0.223]
d3=[0, 0, 0, 0, 0, 0, 0, 0.076, 0.136, 0.184, 0.223, 0.256, 0.284, 0.308, 0.329, 0.348]
d4=[0, 0, 3.268, 2.574, 2.282, 2.114, 2.004, 1.924, 1.864, 1.816, 1.777, 1.744, 1.717, 1.692, 1.671, 1.652]
lcip=(prom1-(a2[c]*rango1))
lcir=(d3[c]*rango1)
lcsp=(prom1+(a2[c]*rango1))
lcsr=(d4[c]*rango1)
print("\tPROM\t RANG")
o16="{}\t {}\t {}\t".format("PROM",prom1,rango1)
o17="{}\t {}\t {}\t".format("LCS",lcsp,lcsr)
o18="{}\t {}\t {}\t".format("LCI",lcip,lcir)
print(o16)
print(o17)
print(o18)
plt.plot(x,label="Resistencias", marker="d")
plt.plot(prom1, label="PROM", marker="s")
plt.plot(lcip, label="LCI", marker="x")
plt.plot(lcsp, label="LCS", marker="h")
plt.legend(loc="upper center")
plt.show()