Ticket #11934: simplify_error.sage

File simplify_error.sage, 1.8 KB (added by mjo, 11 years ago)

A sage file that displays and attempts to plot the simplification

Line 
1x = SR.symbol('x', domain='real')
2
3f = QQ(0.25)*(sqrt(2) - 2)*(x + 1)*x**3 - QQ(3)/QQ(8)*(sqrt(2) - 2)*(x + 1)*x**2 - QQ(0.25)*(sqrt(2) - 2)*(2*(3*sqrt(2) - 2)*x**2 - 2*(sqrt(2) - 1)*x + sqrt(-8*(4*sqrt(2) - 7)*x**4 + 16*(3*sqrt(2) - 5)*x**3 - 8*(2*sqrt(2) - 3)*x - 4*x**2 + 4) - 4*sqrt(2) + 2)*x**3/(sqrt(2)*x**2 + sqrt(2)*x - 2*sqrt(2)) - 1/24*(x + 1)**3*(x**3 - 3*x + 2) + QQ(3)/QQ(8)*(sqrt(2) - 2)*(2*(3*sqrt(2) - 2)*x**2 - 2*(sqrt(2) - 1)*x + sqrt(-8*(4*sqrt(2) - 7)*x**4 + 16*(3*sqrt(2) - 5)*x**3 - 8*(2*sqrt(2) - 3)*x - 4*x**2 + 4) - 4*sqrt(2) + 2)*x**2/(sqrt(2)*x**2 + sqrt(2)*x - 2*sqrt(2)) - QQ(1)/QQ(16)*(x + 1)**2*(2*(sqrt(2) - 3)*x**3 - (3*sqrt(2) - 8)*x**2 + 2*x + sqrt(2) - 4) + QQ(1)/QQ(8)*(x + 1)*sqrt(2) + QQ(1)/QQ(96)*(x**3 - 3*x + 2)*(2*(3*sqrt(2) - 2)*x**2 - 2*(sqrt(2) - 1)*x + sqrt(-8*(4*sqrt(2) - 7)*x**4 + 16*(3*sqrt(2) - 5)*x**3 - 8*(2*sqrt(2) - 3)*x - 4*x**2 + 4) - 4*sqrt(2) + 2)**3/(sqrt(2)*x**2 + sqrt(2)*x - 2*sqrt(2))**3 + 1/32*(2*(3*sqrt(2) - 2)*x**2 - 2*(sqrt(2) - 1)*x + sqrt(-8*(4*sqrt(2) - 7)*x**4 + 16*(3*sqrt(2) - 5)*x**3 - 8*(2*sqrt(2) - 3)*x - 4*x**2 + 4) - 4*sqrt(2) + 2)**2*(2*(sqrt(2) - 3)*x**3 - (3*sqrt(2) - 8)*x**2 + 2*x + sqrt(2) - 4)/(sqrt(2)*x**2 + sqrt(2)*x - 2*sqrt(2))**2 - QQ(0.25)*x - QQ(1)/QQ(8)*(2*(3*sqrt(2) - 2)*x**2 - 2*(sqrt(2) - 1)*x + sqrt(-8*(4*sqrt(2) - 7)*x**4 + 16*(3*sqrt(2) - 5)*x**3 - 8*(2*sqrt(2) - 3)*x - 4*x**2 + 4) - 4*sqrt(2) + 2)*sqrt(2)/(sqrt(2)*x**2 + sqrt(2)*x - 2*sqrt(2)) + QQ(0.25)*(2*(3*sqrt(2) - 2)*x**2 - 2*(sqrt(2) - 1)*x + sqrt(-8*(4*sqrt(2) - 7)*x**4 + 16*(3*sqrt(2) - 5)*x**3 - 8*(2*sqrt(2) - 3)*x - 4*x**2 + 4) - 4*sqrt(2) + 2)/(sqrt(2)*x**2 + sqrt(2)*x - 2*sqrt(2)) - QQ(0.25)
4
5print "f:"
6print f
7print ""
8
9g = f.full_simplify()
10print "f.full_simplify():"
11print g
12print ""
13
14p1  = plot(f, x, -1, 1)
15p1.show()
16
17p2 = plot(g, x, -1, 1, color='red')
18p2.show()