Hiya Burcin, Thanks for reply! For my cubic polynomials p0 and p1, I did get the following to work: p0=p.series(x==xD[0],4) p1=q.series(x==xD[1],4) print p0; p1
x |--> 1.0600000000000001 + 0.59499999999999997*x + (-0.08500000000000002)*x^3 x |--> 1.5699999999999998 + 0.34000000000000008*(x - 1) + (-0.25500000000000006)*(x - 1)^2 + 0.08500000000000002*(x - 1)^3 (I am guessing that here (a) you must add one more degree than you want and (b) if the derivative is not zero, you get the Order thing indicating a remainder??) I did like the taylor command. It was easy to understand and the syntax was intuitive (good for people like me). I don't mind the Order thing so much, but the x==1, I don't like. BTW: The wolframalpha command is: series[0.085*x^3 - 0.510*x^2 + 1.105*x + 0.890,(x,1,3)] I know next to nothing about programming so I hope someone else offers to help :) Thanks again. Linda -- You received this message because you are subscribed to the Google Groups "sage-support" group. To post to this group, send email to sage-support@googlegroups.com. To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com. Visit this group at http://groups.google.com/group/sage-support?hl=en.