NIST: Roszman1 — Defectos cuánticos yodo [Media]

[ ]

Xv = [-4868.68,-4868.09,-4867.41,-3375.19,-3373.14,-3372.03,-2473.74,-2472.35,-2469.45,-1894.65,-1893.40,-1497.24,-1495.85,-1493.41,-1208.68,-1206.18,-1206.04,-997.92,-996.61,-996.31,-834.94,-834.66,-710.03,-530.16,-464.17]
Yv = [0.252429,0.252141,0.251809,0.297989,0.296257,0.295319,0.339603,0.337731,0.333820,0.389510,0.386998,0.438864,0.434887,0.427893,0.471568,0.461699,0.461144,0.513532,0.506641,0.505062,0.535648,0.533726,0.568064,0.612886,0.624169]
b0 = [0.2,-0.000005,1200,-150]

res = []
ajuste_no_lineal(Xv, Yv, 17, b0, res)
b1 = res[1]
b2 = res[2]
b3 = res[3]
b4 = res[4]

b2 = res[2]
b3 = res[3]
b4 = res[4]
mostrar "b1 (cert 0.2020):"
mostrar res[1]
mostrar "b3 (cert 1204.5):"
mostrar res[3]
mostrar "RSS (cert 4.948e-4):"
mostrar res[5]
n_graf = 25.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]
i = 1
para i de 1 hasta 25
    Y_fit[i] = b1 - b2*Xv[i] - arcotangente(b3/(Xv[i]-res[4])) / 3.14159265358979
finpara
graficar_dispersion(Xv, Yv, "Datos observados")
graficar_xy(Xv, Y_fit, "Ajuste NIST")

Salida Jota: b1 (cert 0.2020): Salida Jota: 0.201 Salida Jota: b3 (cert 1204.5): Salida Jota: 1209.380 Salida Jota: RSS (cert 4.948e-4): Salida Jota: 0.000495

Shift + Enter ejecuta la celda actual

Ayuda de Jota