NIST: Thurber — Movilidad electrones [Alta]

[ ]

Xv = [-3.067,-2.981,-2.921,-2.912,-2.84,-2.797,-2.702,-2.699,-2.633,-2.481,-2.363,-2.322,-1.501,-1.46,-1.274,-1.212,-1.1,-1.046,-0.915,-0.714,-0.566,-0.545,-0.4,-0.309,-0.109,-0.103,0.01,0.119,0.377,0.79,0.963,1.006,1.115,1.572,1.841,2.047,2.2]
Yv = [80.574,84.248,87.264,87.195,89.076,89.608,89.868,90.101,92.405,95.854,100.696,101.06,401.672,390.724,567.534,635.316,733.054,759.087,894.206,990.785,1090.109,1080.914,1122.643,1178.351,1260.531,1273.514,1288.339,1327.543,1353.863,1414.509,1425.208,1421.384,1442.962,1464.35,1468.705,1447.894,1457.628]
b0 = [1000,1000,400,40,0.7,0.3,0.03]

res = []
ajuste_no_lineal(Xv, Yv, 11, b0, res)
b1 = res[1]
b2 = res[2]
b3 = res[3]
b4 = res[4]
b5 = res[5]
b6 = res[6]
b7 = res[7]

b2 = res[2]
b3 = res[3]
b4 = res[4]
b5 = res[5]
b6 = res[6]
b7 = res[7]
mostrar "b1 (cert 1288.14):"
mostrar res[1]
mostrar "RSS (cert 5642.7):"
mostrar res[8]
n_graf = 37.0
Y_fit = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
i = 1
para i de 1 hasta 37
    Y_fit[i] = (b1 + b2*Xv[i] + b3*Xv[i]*Xv[i] + b4*Xv[i]*Xv[i]*Xv[i]) / (1.0 + b5*Xv[i] + b6*Xv[i]*Xv[i] + b7*Xv[i]*Xv[i]*Xv[i])
finpara
graficar_dispersion(Xv, Yv, "Datos observados")
graficar_xy(Xv, Y_fit, "Ajuste NIST")

Salida Jota: b1 (cert 1288.14): Salida Jota: 1288.140 Salida Jota: RSS (cert 5642.7): Salida Jota: 5642.708

Shift + Enter ejecuta la celda actual

Ayuda de Jota