def write_table(ex,f_re): global inner inner+=ocd.table_groesse(g,var,ex+offsets,"g") inner+=ocd.table_groesse(b,var,ex+offsets,"b") inner+=ocd.table_groesse(f,var,ex+offsets,"e") inner+=ocd.table_groesse(f,var,ex+offsets,"f") inner+=ocd.table_groesse(beta,var,ex+offsets,ur"beta") f_r= ocd.eval_expr(f,var,ex+offsets,"f") f_mean=np.mean(f_r.x.magnitude) f_std=np.std(f_r.x.magnitude) inner+=ur"Mittelwert: %.3f\\"%f_mean inner+=ur"Standartabweichung: %.3f\\" % np.std(f_r.x.magnitude) inner+=ur"gemessene Brechkraft: $\phi=\frac{1}{f}=%.3f +- %.3f\\ $"% ((1/f_mean),f_std/f_mean**2) inner+=ur"theoretische Brechkraft: %.3f\\ "% (1/f_re)
beta = b / g f = g*b/ (g+b) inner+=ur"Experiment II 1. \\ " write_table(ex2,0.0574) inner+=ur"Experiment II 2. \\ " write_table(ex3,0.0574) inner+=ur"Experiment III 1.\\ " B=sy.Symbol("B") B_orginal=sy.Symbol("B_orginal") b__=sy.Symbol("b__") b_=e_-b__+xs-xk var+=[B,b__] beta=B/B_orginal inner+=ocd.table_groesse(b_,var,ex4+offsets,ur"b'") inner+=ocd.table_groesse(b_,var,ex5+offsets,ur"b'") inner+=ocd.table_groesse(beta,var,ex4+offsets,ur"beta") inner+=ocd.table_groesse(beta,var,ex5+offsets,ur"beta") plt.figure() (a,b,Sa,Sb,Sy)=ocd.plot_var(b_,beta,var,ex4+offsets,True) (a,b,Sa,Sb,Sy)=ocd.plot_var(b_,beta,var,ex5+offsets,True) plt.xlim(0,0.9) plt.ylim(0,7) plt.xlabel(ur"b' in m ") plt.ylabel(ur"$\beta$") plt.title(ur"gemessene Brechkraft: $\phi=(%.3f \pm %.3f) \frac{1}{m}$ theoretische Brechkraft: $\phi=%.3f \frac{1}{m}$"%(b,(Sb),(1/0.11429))) plt.show() f1=sy.Symbol("f1") fges=sy.Symbol("fges")