4350.sagews

num = 10
pesos = [1.96, 3.92, 5.88, 7.84, 9.80, 11.76, 13.72, 15.68, 17.64, 19.6]
flechas = [0.00055, 0.00112, 0.00170, 0.00228, 0.00292, 0.00349, 0.00409, 0.00453, 0.00506, 0.00560]

a = (num*(sum(pesos[i]*flechas[i] for i in range(num))) - (sum(pesos[i] for i in range(num)))*(sum(flechas[i] for i in range(num))))/ (num*(sum(pesos[i]^2 for i in range(num))) - (sum(pesos[i] for i in range(num)))^2)

b = ((sum(pesos[i]^2 for i in range(num)))*(sum(flechas[i] for i in range(num))) - (sum(pesos[i] for i in range(num)))*(sum(pesos[i]*flechas[i] for i in range(num))))/ ((num*(sum(pesos[i]^2 for i in range(num))))-(sum(pesos[i] for i in range(num)))^2)

y=a*x+b

n(a)

0.000288126159554731

n(b)

0.0000279999999999986

espesor = 0.0046
anchura = 0.01264
longitud = 0.475
inercia = (espesor^3)*(anchura)/12
E = (longitud^3)/((48)*(10^9)*(a)*(inercia))
n(E)

75.5818604722929

graph1 = list_plot(zip(pesos, flechas), axes_labels=['$F\;(N)$','$\delta\;(m)$'], size=20, fontsize=20, color='black')
graph1.save('graph1.pdf')
graph1

graph2 = list_plot(zip(pesos, flechas), axes_labels=['$F\;(N)$','$\delta\;(m)$'], size=20, fontsize=20, color='black')+plot(y,(x,0,20),color='black')
graph2.save('graph2.pdf')
graph2