| 1 | initial version |
LA(x) = x - (x^3.0-1.0)/(3.0*x^2.0)
S = solve(LA==0,x)
n=0
for s in S:
a=solve(LA == s.rhs(),x)
n=n+1
point([a[i].rhs().n().real(),a[i].rhs().n().imag()] for i in range(n))
Does this help?
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.