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")
