Revision history [back]

You could expand around t=0 instead, for t = w + 0.

var('t')
bmac = b.substitute(w=t-0.1).taylor(t,0,3)
cmac = c.substitute(w=t-0.1).taylor(t,0,3)
print(bmac)
print(cmac)

gives

(77.59833119183804 - 60.0058944656944I)t^3 - (8.236463209769589 - 3.574042330177453I)t^2 + (0.8150178993589599 - 0.1780426145916193I)t + 0.0810824806541207

(77.59833119183804 - 60.0058944656944I)t^3 - (8.236463209769589 - 3.574042330177453I)t^2 + (0.8150178993589599 - 0.1780426145916193I)t + 0.0810824806541207

​

You could expand around t=0 instead, for t = w + 0.0.1

var('t')
bmac = b.substitute(w=t-0.1).taylor(t,0,3)
cmac = c.substitute(w=t-0.1).taylor(t,0,3)
print(bmac)
print(cmac)

gives

(77.59833119183804 - 60.0058944656944I)t^3 - (8.236463209769589 - 3.574042330177453I)t^2 + (0.8150178993589599 - 0.1780426145916193I)t + 0.0810824806541207

(77.59833119183804 - 60.0058944656944I)t^3 - (8.236463209769589 - 3.574042330177453I)t^2 + (0.8150178993589599 - 0.1780426145916193I)t + 0.0810824806541207

​