/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/ /* [ Created with wxMaxima version 11.04.0 ] */ /* [wxMaxima: comment start ] L. 1.9, Anikushina [wxMaxima: comment end ] */ /* [wxMaxima: comment start ] Analytics [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ expression_g: 4*(%pi)^2 * ((1/12)*'l^2 + 'a^2) / ('a * 'Te^2); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ expression_err_g: ratsimp((abs(diff(expression_g, 'a)/expression_g))*'Dl) + ratsimp((abs(diff(expression_g, 'Te)/expression_g/'Nosc))*'Dt) + ratsimp((abs(diff(expression_g, 'l)/expression_g))*'Dl); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Calculations [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ fpprintprec: 3$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Units [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ meter:1$ s:1$ kg:1$ g:9.81*meter/s^2$ cm:meter/100$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Constant data [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ l: 149*cm$ Dl: 0.5*cm$ Dt: 0.2*s$ Nosc: 40$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Functions [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ f_Te(t):=t/''Nosc$ f_Tt(g,a):=2 * %pi *((1/12*''l^2+a^2)/(g*a))^0.5$ f_g(a,Te):= ''expression_g$ f_err_g(a,Te):= ''expression_err_g$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Experimental data [wxMaxima: comment end ] */ /* [wxMaxima: comment start ] Number of measurements [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Nmeas: 16$ /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ s1: zeromatrix(Nmeas,6)$ /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ s1[1]:[8*cm, 133*s, 0, 0, 0, 0]$ s1[2]:[12*cm, 110*s, 0, 0, 0, 0]$ s1[3]:[16*cm, 98.5*s, 0, 0, 0, 0]$ s1[4]:[20*cm, 90.5*s, 0, 0, 0, 0]$ s1[5]:[24*cm, 82.5*s, 0, 0, 0, 0]$ s1[6]:[28*cm, 82*s, 0, 0, 0, 0]$ s1[7]:[32*cm, 80*s, 0, 0, 0, 0]$ s1[8]:[36*cm, 78.5*s, 0, 0, 0, 0]$ s1[9]:[40*cm, 77.5*s, 0, 0, 0, 0]$ s1[10]:[44*cm, 77.5*s, 0, 0, 0, 0]$ s1[11]:[48*cm, 77.5*s, 0, 0, 0, 0]$ s1[12]:[52*cm, 77.5*s, 0, 0, 0, 0]$ s1[13]:[56*cm, 78*s, 0, 0, 0, 0]$ s1[14]:[60*cm, 79*s, 0, 0, 0, 0]$ s1[15]:[64*cm, 79*s, 0, 0, 0, 0]$ s1[16]:[68*cm, 80*s, 0, 0, 0, 0]$ for i thru Nmeas do ( s1[i][3]: f_Te(s1[i][2]), s1[i][4]: f_Tt(''g, s1[i][1]), s1[i][5]: f_g(s1[i][1], s1[i][3]), s1[i][6]: f_err_g(s1[i][1], s1[i][3]) ); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] a, t, Te, Tt, g, err_g [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ float(s1); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ submatrix(s1,2,4,5)$ data1:maplist(copy,%)$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Here we change realative errors to absolute errors in data1 [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ for i thru Nmeas do (data1[i][3]:data1[i][2] * data1[i][3]); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Experimental value of g [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ T_vs_a:submatrix(s1,2,4,5)$ lsquares_estimates(T_vs_a,['x,y], y=f_Tt(ggg, 'x),[ggg])$ g_exp:ggg,%$ float(g_exp); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ wxdraw2d( xlabel="a, m", ylabel="T, s", xaxis=true, xaxis_type=dots/*solid*/, grid=true, color=red, key="Experimental", explicit(f_Tt(''g_exp,'a),'a,0.1,0.7), color=green, key="Theoretical", explicit(f_Tt(''g,'a),'a,0.1,0.7), color=blue, key="Experimental", errors(data1))$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Function 'set_draw_defaults' may be used to set default values for wxdraw2d. [wxMaxima: comment end ] */ /* Maxima can't load/batch files which end with a comment! */ "Created with wxMaxima"$