code2.sagews

%load code.sage

show(symbolic_cos_theta(100))

$$\cos\left(98.831194218 \, t\right) + \cos\left(95.870634228 \, t\right) + \cos\left(94.651344041 \, t\right) + \cos\left(92.491899271 \, t\right) + \cos\left(88.809111208 \, t\right) + \cos\left(87.425274613 \, t\right) + \cos\left(84.735492981 \, t\right) + \cos\left(82.910380854 \, t\right) + \cos\left(79.33737502 \, t\right) + \cos\left(77.144840069 \, t\right) + \cos\left(75.704690699 \, t\right) + \cos\left(72.067157674 \, t\right) + \cos\left(69.546401711 \, t\right) + \cos\left(67.079810529 \, t\right) + \cos\left(65.112544048 \, t\right) + \cos\left(60.831778525 \, t\right) + \cos\left(59.347044003 \, t\right) + \cos\left(56.446247697 \, t\right) + \cos\left(52.970321478 \, t\right) + \cos\left(49.773832478 \, t\right) + \cos\left(48.005150881 \, t\right) + \cos\left(43.327073281 \, t\right) + \cos\left(40.918719012 \, t\right) + \cos\left(37.586178159 \, t\right) + \cos\left(32.935061588 \, t\right) + \cos\left(30.424876126 \, t\right) + \cos\left(25.01085758 \, t\right) + \cos\left(21.022039639 \, t\right) + \cos\left(14.134725142 \, t\right) + 1$$

g = plot(symbolic_cos_theta(100), 0, 40)

%time g = plot_symbolic_cos_theta(100, 0, 40)

CPU time: 3.42 s, Wall time: 3.43 s

g

g = symbolic_cos_theta(1000)
t = var('t')
h = g(t=log(t))
%time plot(h, 5, 40)

CPU time: 223.68 s, Wall time: 224.22 s