A Sample Session using SymPy#
In this first part, we do all of the examples in the SymPy tutorial (https://github.com/sympy/sympy/wiki/Tutorial), but using Sage instead of SymPy.
sage: a = Rational((1,2))
sage: a
1/2
sage: a*2
1
sage: Rational(2)^50 / Rational(10)^50
1/88817841970012523233890533447265625
sage: 1.0/2
0.500000000000000
sage: 1/2
1/2
sage: pi^2
pi^2
sage: float(pi)
3.141592653589793
sage: RealField(200)(pi)
3.1415926535897932384626433832795028841971693993751058209749
sage: float(pi + exp(1))
5.85987448204883...
sage: oo != 2
True
>>> from sage.all import *
>>> a = Rational((Integer(1),Integer(2)))
>>> a
1/2
>>> a*Integer(2)
1
>>> Rational(Integer(2))**Integer(50) / Rational(Integer(10))**Integer(50)
1/88817841970012523233890533447265625
>>> RealNumber('1.0')/Integer(2)
0.500000000000000
>>> Integer(1)/Integer(2)
1/2
>>> pi**Integer(2)
pi^2
>>> float(pi)
3.141592653589793
>>> RealField(Integer(200))(pi)
3.1415926535897932384626433832795028841971693993751058209749
>>> float(pi + exp(Integer(1)))
5.85987448204883...
>>> oo != Integer(2)
True
sage: var('x y')
(x, y)
sage: x + y + x - y
2*x
sage: (x+y)^2
(x + y)^2
sage: ((x+y)^2).expand()
x^2 + 2*x*y + y^2
sage: ((x+y)^2).subs(x=1)
(y + 1)^2
sage: ((x+y)^2).subs(x=y)
4*y^2
>>> from sage.all import *
>>> var('x y')
(x, y)
>>> x + y + x - y
2*x
>>> (x+y)**Integer(2)
(x + y)^2
>>> ((x+y)**Integer(2)).expand()
x^2 + 2*x*y + y^2
>>> ((x+y)**Integer(2)).subs(x=Integer(1))
(y + 1)^2
>>> ((x+y)**Integer(2)).subs(x=y)
4*y^2
sage: limit(sin(x)/x, x=0)
1
sage: limit(x, x=oo)
+Infinity
sage: limit((5^x + 3^x)^(1/x), x=oo)
5
>>> from sage.all import *
>>> limit(sin(x)/x, x=Integer(0))
1
>>> limit(x, x=oo)
+Infinity
>>> limit((Integer(5)**x + Integer(3)**x)**(Integer(1)/x), x=oo)
5
sage: diff(sin(x), x)
cos(x)
sage: diff(sin(2*x), x)
2*cos(2*x)
sage: diff(tan(x), x)
tan(x)^2 + 1
sage: limit((tan(x+y) - tan(x))/y, y=0)
cos(x)^(-2)
sage: diff(sin(2*x), x, 1)
2*cos(2*x)
sage: diff(sin(2*x), x, 2)
-4*sin(2*x)
sage: diff(sin(2*x), x, 3)
-8*cos(2*x)
>>> from sage.all import *
>>> diff(sin(x), x)
cos(x)
>>> diff(sin(Integer(2)*x), x)
2*cos(2*x)
>>> diff(tan(x), x)
tan(x)^2 + 1
>>> limit((tan(x+y) - tan(x))/y, y=Integer(0))
cos(x)^(-2)
>>> diff(sin(Integer(2)*x), x, Integer(1))
2*cos(2*x)
>>> diff(sin(Integer(2)*x), x, Integer(2))
-4*sin(2*x)
>>> diff(sin(Integer(2)*x), x, Integer(3))
-8*cos(2*x)
sage: cos(x).taylor(x,0,10)
-1/3628800*x^10 + 1/40320*x^8 - 1/720*x^6 + 1/24*x^4 - 1/2*x^2 + 1
sage: (1/cos(x)).taylor(x,0,10)
50521/3628800*x^10 + 277/8064*x^8 + 61/720*x^6 + 5/24*x^4 + 1/2*x^2 + 1
>>> from sage.all import *
>>> cos(x).taylor(x,Integer(0),Integer(10))
-1/3628800*x^10 + 1/40320*x^8 - 1/720*x^6 + 1/24*x^4 - 1/2*x^2 + 1
>>> (Integer(1)/cos(x)).taylor(x,Integer(0),Integer(10))
50521/3628800*x^10 + 277/8064*x^8 + 61/720*x^6 + 5/24*x^4 + 1/2*x^2 + 1
sage: matrix([[1,0], [0,1]])
[1 0]
[0 1]
sage: var('x y')
(x, y)
sage: A = matrix([[1,x], [y,1]])
sage: A
[1 x]
[y 1]
sage: A^2
[x*y + 1 2*x]
[ 2*y x*y + 1]
sage: R.<x,y> = QQ[]
sage: A = matrix([[1,x], [y,1]])
sage: A^10
[x^5*y^5 + 45*x^4*y^4 + 210*x^3*y^3 + 210*x^2*y^2 + 45*x*y + 1 10*x^5*y^4 + 120*x^4*y^3 + 252*x^3*y^2 + 120*x^2*y + 10*x]
[ 10*x^4*y^5 + 120*x^3*y^4 + 252*x^2*y^3 + 120*x*y^2 + 10*y x^5*y^5 + 45*x^4*y^4 + 210*x^3*y^3 + 210*x^2*y^2 + 45*x*y + 1]
sage: var('x y')
(x, y)
>>> from sage.all import *
>>> matrix([[Integer(1),Integer(0)], [Integer(0),Integer(1)]])
[1 0]
[0 1]
>>> var('x y')
(x, y)
>>> A = matrix([[Integer(1),x], [y,Integer(1)]])
>>> A
[1 x]
[y 1]
>>> A**Integer(2)
[x*y + 1 2*x]
[ 2*y x*y + 1]
>>> R = QQ['x, y']; (x, y,) = R._first_ngens(2)
>>> A = matrix([[Integer(1),x], [y,Integer(1)]])
>>> A**Integer(10)
[x^5*y^5 + 45*x^4*y^4 + 210*x^3*y^3 + 210*x^2*y^2 + 45*x*y + 1 10*x^5*y^4 + 120*x^4*y^3 + 252*x^3*y^2 + 120*x^2*y + 10*x]
[ 10*x^4*y^5 + 120*x^3*y^4 + 252*x^2*y^3 + 120*x*y^2 + 10*y x^5*y^5 + 45*x^4*y^4 + 210*x^3*y^3 + 210*x^2*y^2 + 45*x*y + 1]
>>> var('x y')
(x, y)
And here are some actual tests of sympy:
sage: from sympy import Symbol, cos, sympify, pprint # needs sympy
sage: from sympy.abc import x # needs sympy
>>> from sage.all import *
>>> from sympy import Symbol, cos, sympify, pprint # needs sympy
>>> from sympy.abc import x # needs sympy
sage: e = (1/cos(x)^3)._sympy_(); e # needs sympy
cos(x)**(-3)
sage: f = e.series(x, 0, int(10)); f # needs sympy
1 + 3*x**2/2 + 11*x**4/8 + 241*x**6/240 + 8651*x**8/13440 + O(x**10)
>>> from sage.all import *
>>> e = (Integer(1)/cos(x)**Integer(3))._sympy_(); e # needs sympy
cos(x)**(-3)
>>> f = e.series(x, Integer(0), int(Integer(10))); f # needs sympy
1 + 3*x**2/2 + 11*x**4/8 + 241*x**6/240 + 8651*x**8/13440 + O(x**10)
And the pretty-printer. Since unicode characters are not working on some architectures, we disable it:
sage: from sympy.printing import pprint_use_unicode # needs sympy
sage: prev_use = pprint_use_unicode(False) # needs sympy
sage: pprint(e) # needs sympy
1
-------
3
cos (x)
sage: pprint(f) # needs sympy
2 4 6 8
3*x 11*x 241*x 8651*x / 10\
1 + ---- + ----- + ------ + ------- + O\x /
2 8 240 13440
sage: pprint_use_unicode(prev_use) # needs sympy
False
>>> from sage.all import *
>>> from sympy.printing import pprint_use_unicode # needs sympy
>>> prev_use = pprint_use_unicode(False) # needs sympy
>>> pprint(e) # needs sympy
1
-------
3
cos (x)
>>> pprint(f) # needs sympy
2 4 6 8
3*x 11*x 241*x 8651*x / 10\
1 + ---- + ----- + ------ + ------- + O\x /
2 8 240 13440
>>> pprint_use_unicode(prev_use) # needs sympy
False
And the functionality to convert from sympy format to Sage format:
sage: e._sage_() # needs sympy
cos(x)^(-3)
sage: e._sage_().taylor(x._sage_(), 0, 8) # needs sympy
8651/13440*x^8 + 241/240*x^6 + 11/8*x^4 + 3/2*x^2 + 1
sage: f._sage_() # needs sympy
8651/13440*x^8 + 241/240*x^6 + 11/8*x^4 + 3/2*x^2 + Order(x^10) + 1
>>> from sage.all import *
>>> e._sage_() # needs sympy
cos(x)^(-3)
>>> e._sage_().taylor(x._sage_(), Integer(0), Integer(8)) # needs sympy
8651/13440*x^8 + 241/240*x^6 + 11/8*x^4 + 3/2*x^2 + 1
>>> f._sage_() # needs sympy
8651/13440*x^8 + 241/240*x^6 + 11/8*x^4 + 3/2*x^2 + Order(x^10) + 1
Mixing SymPy with Sage:
sage: # needs sympy
sage: import sympy
sage: var("x")._sympy_() + var("y")._sympy_()
x + y
sage: o = var("omega")
sage: s = sympy.Symbol("x")
sage: t1 = s + o
sage: t2 = o + s
sage: type(t1)
<class 'sympy.core.add.Add'>
sage: type(t2)
<class 'sage.symbolic.expression.Expression'>
sage: t1, t2
(omega + x, omega + x)
sage: e = sympy.sin(var("y")) + sage.functions.trig.cos(sympy.Symbol("x"))
sage: type(e)
<class 'sympy.core.add.Add'>
sage: e
sin(y) + cos(x)
sage: e=e._sage_()
sage: type(e)
<class 'sage.symbolic.expression.Expression'>
sage: e
cos(x) + sin(y)
sage: e = sage.functions.trig.cos(var("y")**3)**4+var("x")**2
sage: e = e._sympy_()
sage: e
x**2 + cos(y**3)**4
>>> from sage.all import *
>>> # needs sympy
>>> import sympy
>>> var("x")._sympy_() + var("y")._sympy_()
x + y
>>> o = var("omega")
>>> s = sympy.Symbol("x")
>>> t1 = s + o
>>> t2 = o + s
>>> type(t1)
<class 'sympy.core.add.Add'>
>>> type(t2)
<class 'sage.symbolic.expression.Expression'>
>>> t1, t2
(omega + x, omega + x)
>>> e = sympy.sin(var("y")) + sage.functions.trig.cos(sympy.Symbol("x"))
>>> type(e)
<class 'sympy.core.add.Add'>
>>> e
sin(y) + cos(x)
>>> e=e._sage_()
>>> type(e)
<class 'sage.symbolic.expression.Expression'>
>>> e
cos(x) + sin(y)
>>> e = sage.functions.trig.cos(var("y")**Integer(3))**Integer(4)+var("x")**Integer(2)
>>> e = e._sympy_()
>>> e
x**2 + cos(y**3)**4
sage: a = sympy.Matrix([1, 2, 3]) # needs sympy
sage: a[1] # needs sympy
2
>>> from sage.all import *
>>> a = sympy.Matrix([Integer(1), Integer(2), Integer(3)]) # needs sympy
>>> a[Integer(1)] # needs sympy
2
sage: sympify(1.5) # needs sympy
1.50000000000000
sage: sympify(2) # needs sympy
2
sage: sympify(-2) # needs sympy
-2
>>> from sage.all import *
>>> sympify(RealNumber('1.5')) # needs sympy
1.50000000000000
>>> sympify(Integer(2)) # needs sympy
2
>>> sympify(-Integer(2)) # needs sympy
-2