1324.sagews

y=var('y')

plot(12.5/cos((2*y*pi)/360.0),(y,0,60))

v = [(i,12.5/cos((2*(i)*pi)/360.0)) for i in range(0,60,2)]
plot=plot_step_function(v)
plot.axes_labels(['Effective thickness [$\mu$m]','Degrator angle[degrees]'])
plot.show()

for i in range(0,60,2):
    equa=(12.5/cos((2*(i)*3.14159265)/360.0))
    print i,         equa

0 12.5000000000000 2 12.5076193037179 4 12.5305237259448 6 12.5688534943859 8 12.6228446561997 10 12.6928326481255 12 12.7792574351628 14 12.8826703659769 16 13.0037429470802 18 13.1432778014453 20 13.3022221540177 22 13.4816842810799 24 13.6829534784097 26 13.9075242524226 28 14.1571256294096 30 14.4337567247548 32 14.7397300361483 34 15.0777243494028 36 15.4508497106878 38 15.8627276790150 40 16.3175911057309 42 16.8204091073938 44 17.3770448729832 46 17.9944567282271 48 18.6809568534467 50 19.4465478126454 52 20.3033655415843 54 21.2662701772785 56 22.3536455876302 58 23.5884988913330

\newcommand{\Bold}[1]{\mathbf{#1}}12.5000000000000 \, \frac{1}{\cos\left(0.322222222222222 \, \pi\right)}

import numpy
deg=var('deg')
nurange = [(deg,acos(12.5/deg)*360.0/(2*pi)) for deg in numpy.arange(12.5,25,0.3)]
plot=plot_step_function(nurange)
plot.axes_labels(['Effective thickness [$\mu$m]','Degrator angle[degrees]'])
plot.show()

Traceback (most recent call last): File "<stdin>", line 1, in <module> File "_sage_input_11.py", line 12, in <module> exec compile(ur'plot(acos(_sage_const_12p5 /t)*_sage_const_360p0 /(_sage_const_2 *pi),(t,_sage_const_12p5 , _sage_const_25 ))' + '\n', '', 'single') File "", line 1, in <module> TypeError: 'Graphics' object is not callable

for deg in numpy.arange(12.5,25,0.3):
    print deg, acos(12.5/deg)*360.0/(2*3.14159265)

12.5 0.0 12.8 12.4292571467 13.1 17.4080176316 13.4 21.1187270131 13.7 24.1596241546 14.0 26.7655006073 14.3 29.0583087 14.6 31.1113901364 14.9 32.9731006201 15.2 34.6774106277 15.5 36.2493209243 15.8 37.7078969493 16.1 39.0680915284 16.4 40.3419010235 16.7 41.539131564 17.0 42.6679255429 17.3 43.7351345113 17.6 44.7465901803 17.9 45.7073057912 18.2 46.6216286663 18.5 47.4933577516 18.8 48.3258355555 19.1 49.1220210299 19.4 49.8845480391 19.7 50.615772778 20.0 51.3178126051 20.3 51.9925781295 20.6 52.64179994 20.9 53.2670510376 21.2 53.8697657916 21.5 54.4512560593 21.8 55.0127249746 22.1 55.5552788062 22.4 56.079937208 22.7 56.5876421208 23.0 57.0792655371 23.3 57.5556163027 23.6 58.0174460959 23.9 58.4654547055 24.2 58.9002947033 24.5 59.3225755952 24.8 59.7328675206

t=var('t')
nuplot=plot(acos(12.5/t)*360.0/(2*pi),(t,12.6,25))
nuplot.axes_labels(['Effective thickness [$\mu$m]','Degrator angle[degrees]'])
nuplot.show()

2+4

\newcommand{\Bold}[1]{\mathbf{#1}}6

x^3

x^3