passagemath-repl 10.5.28__py3-none-any.whl → 10.5.30__py3-none-any.whl

sage/tests/books/computational-mathematics-with-sagemath/sol/combinat_doctest.py

@@ -0,0 +1,218 @@
+# sage_setup: distribution = sagemath-repl
+# sage.doctest: needs sage.combinat sage.graphs sage.modules
+"""
+This file (./sol/combinat_doctest.sage) was *autogenerated* from ./sol/combinat.tex,
+with sagetex.sty version 2011/05/27 v2.3.1.
+It contains the contents of all the sageexample environments from this file.
+You should be able to doctest this file with:
+sage -t ./sol/combinat_doctest.sage
+It is always safe to delete this file; it is not used in typesetting your
+document.
+
+Sage example in ./sol/combinat.tex, line 15::
+
+  sage: Suits = FiniteEnumeratedSet(
+  ....:            ["Hearts", "Diamonds", "Spades", "Clubs"])
+  sage: Values  = FiniteEnumeratedSet([2, 3, 4, 5, 6, 7, 8, 9, 10,
+  ....:                                "Jack", "Queen", "King", "Ace"])
+  sage: FourOfaKind = cartesian_product([Arrangements(Values,2), Suits])
+
+Sage example in ./sol/combinat.tex, line 34::
+
+  sage: FourOfaKind.list()
+  [([2, 3], 'Hearts'),
+   ([2, 3], 'Diamonds'),
+  ...
+   (['Ace', 'King'], 'Clubs')]
+
+Sage example in ./sol/combinat.tex, line 51::
+
+  sage: FourOfaKind.cardinality()
+  624
+
+Sage example in ./sol/combinat.tex, line 66::
+
+  sage: Cards = cartesian_product([Values, Suits])
+  sage: Hands = Subsets(Cards, 5)
+  sage: FourOfaKind.cardinality() / Hands.cardinality()
+  1/4165
+
+Sage example in ./sol/combinat.tex, line 86::
+
+  sage: StraightFlush = cartesian_product([range(1, 11), Suits])
+  sage: StraightFlush.cardinality()
+  40
+
+Sage example in ./sol/combinat.tex, line 98::
+
+  sage: AllFlush = cartesian_product([Subsets(Values,5),Suits])
+  sage: AllFlush.cardinality() - StraightFlush.cardinality()
+  5108
+
+Sage example in ./sol/combinat.tex, line 113::
+
+  sage: _ / Hands.cardinality()
+  1277/649740
+  sage: float(_)
+  0.001965401545233478
+
+Sage example in ./sol/combinat.tex, line 153::
+
+  sage: Word(['a','b','b','a','a','b','a']).evaluation_dict()
+  {'a': 4, 'b': 3}
+
+Sage example in ./sol/combinat.tex, line 170::
+
+  sage: def is_full_hand(hand):
+  ....:     suits = Word([value for (value, suit) in hand])
+  ....:     repetitions = sorted(suits.evaluation_dict().values())
+  ....:     return repetitions == [2,3]
+  sage: is_full_hand({(5, 'Diamonds'), (6, 'Diamonds'), (6, 'Hearts'),
+  ....:               (5, 'Spades'), (1, 'Spades')})
+  False
+  sage: is_full_hand({(3, 'Clubs'), (3, 'Spades'), (3, 'Hearts'),
+  ....:               (2, 'Clubs'), (2, 'Spades')})
+  True
+
+Sage example in ./sol/combinat.tex, line 200::
+
+  sage: def estimate_proportion(S, predicate, n):
+  ....:     count = 0
+  ....:     for i in range(n):
+  ....:         if predicate(S.random_element()):
+  ....:             count += 1
+  ....:     return count/n
+
+Sage example in ./sol/combinat.tex, line 212::
+
+  sage: float(estimate_proportion(Hands, is_full_hand, 10000)) # random
+  0.0014
+
+Sage example in ./sol/combinat.tex, line 229::
+
+  sage: FullHands = cartesian_product([Arrangements(Values, 2),
+  ....:     Subsets(Suits, 3), Subsets(Suits, 2)])
+
+Sage example in ./sol/combinat.tex, line 242::
+
+  sage: [sorted(v) for v in FullHands.first()]
+  [[2, 3], ['Diamonds', 'Hearts', 'Spades'], ['Diamonds', 'Hearts']]
+
+Sage example in ./sol/combinat.tex, line 252::
+
+  sage: float(FullHands.cardinality() / Hands.cardinality())
+  0.0014405762304921968
+
+Sage example in ./sol/combinat.tex, line 312::
+
+  sage: C = Compositions(5,length=3)
+  sage: C.cardinality
+  <bound method ..._cardinality_from_iterator ...>
+
+Sage example in ./sol/combinat.tex, line 322::
+
+  sage: IntegerVectors(5,3).list()
+  [[5, 0, 0], [4, 1, 0], [4, 0, 1], [3, 2, 0], [3, 1, 1], [3, 0, 2],
+   ...
+   [0, 4, 1], [0, 3, 2], [0, 2, 3], [0, 1, 4], [0, 0, 5]]
+
+Sage example in ./sol/combinat.tex, line 328::
+
+  sage: OrderedSetPartitions(3).cardinality()
+  13
+  sage: OrderedSetPartitions(3).list()
+  [[{1}, {2}, {3}], [{1}, {3}, {2}], [{2}, {1}, {3}], [{3}, {1}, {2}],
+   ...
+   [{1, 2}, {3}], [{1, 3}, {2}], [{2, 3}, {1}], [{1, 2, 3}]]
+  sage: OrderedSetPartitions(3,2).random_element()    # random
+  [{1, 3}, {2}]
+
+Sage example in ./sol/combinat.tex, line 338::
+
+  sage: StandardTableaux([3,2]).cardinality()
+  5
+  sage: StandardTableaux([3,2]).an_element()
+  [[1, 3, 5], [2, 4]]
+
+Sage example in ./sol/combinat.tex, line 350::
+
+  sage: list(AlternatingSignMatrices(1))
+  [[1]]
+  sage: list(AlternatingSignMatrices(2))
+  [
+  [1 0]  [0 1]
+  [0 1], [1 0]
+  ]
+
+Sage example in ./sol/combinat.tex, line 361::
+
+  sage: list(AlternatingSignMatrices(3))
+  [
+  [1 0 0]  [0 1 0]  [1 0 0]  [ 0  1  0]  [0 0 1]  [0 1 0]  [0 0 1]
+  [0 1 0]  [1 0 0]  [0 0 1]  [ 1 -1  1]  [1 0 0]  [0 0 1]  [0 1 0]
+  [0 0 1], [0 0 1], [0 1 0], [ 0  1  0], [0 1 0], [1 0 0], [1 0 0]
+  ]
+
+Sage example in ./sol/combinat.tex, line 383::
+
+  sage: GF(2)^5
+  Vector space of dimension 5 over Finite Field of size 2
+  sage: _.cardinality()
+  32
+
+Sage example in ./sol/combinat.tex, line 400::
+
+  sage: (2^3-2^0)*(2^3-2^1)*(2^3-2^2)
+  168
+
+Sage example in ./sol/combinat.tex, line 406::
+
+  sage: GL(3,2)
+  General Linear Group of degree 3 over Finite Field of size 2
+  sage: _.cardinality()
+  168
+
+Sage example in ./sol/combinat.tex, line 423::
+
+  sage: from sage.combinat.q_analogues import q_factorial
+  sage: q = 2; n = 3
+  sage: q^(n*(n-1)/2) * (q-1)^n * q_factorial(n,q)
+  168
+  sage: q = 3; n = 5
+  sage: GL(n, q).cardinality()
+  475566474240
+  sage: q^(n*(n-1)/2) * (q-1)^n * q_factorial(n,q)
+  475566474240
+
+Sage example in ./sol/combinat.tex, line 537::
+
+  sage: def C(n):
+  ....:     if n == 1:
+  ....:         yield BinaryTree()
+  ....:     elif n > 1:
+  ....:         for k in range(1,n):
+  ....:             for t1 in C(k):
+  ....:                 for t2 in C(n-k):
+  ....:                     yield BinaryTree([t1,t2])
+
+Sage example in ./sol/combinat.tex, line 549::
+
+  sage: list(C(1))
+  [.]
+  sage: list(C(2))
+  [[., .]]
+  sage: list(C(3))
+  [[., [., .]],
+   [[., .], .]]
+  sage: list(C(4))
+  [[., [., [., .]]],
+   [., [[., .], .]],
+   [[., .], [., .]],
+   [[., [., .]], .],
+   [[[., .], .], .]]
+
+Sage example in ./sol/combinat.tex, line 566::
+
+  sage: [len(list(C(n))) for n in range(9)]
+  [0, 1, 1, 2, 5, 14, 42, 132, 429]
+"""

sage/tests/books/computational-mathematics-with-sagemath/sol/domaines_doctest.py

@@ -0,0 +1,60 @@
+# sage_setup: distribution = sagemath-repl
+# sage.doctest: needs sage.symbolic
+"""
+This file (./sol/domaines_doctest.sage) was *autogenerated* from ./sol/domaines.tex,
+with sagetex.sty version 2011/05/27 v2.3.1.
+It contains the contents of all the sageexample environments from this file.
+You should be able to doctest this file with:
+sage -t ./sol/domaines_doctest.sage
+It is always safe to delete this file; it is not used in typesetting your
+document.
+
+Sage example in ./sol/domaines.tex, line 4::
+
+  sage: def ndigits(x): return x.ndigits()
+  sage: o = 720; ndigits(o)
+  3
+
+Sage example in ./sol/domaines.tex, line 28::
+
+  sage: a = Reals(17)(pi); b = Reals(42)(pi)
+  sage: type(a) == type(b)
+  True
+  sage: parent(a), parent(b)
+  (Real Field with 17 bits of precision, Real Field with 42 bits of precision)
+
+Sage example in ./sol/domaines.tex, line 46::
+
+  sage: a = 0.1; b = 0.1*1
+  sage: type(a), type(b)
+  (<... 'sage.rings.real_mpfr.RealLiteral'>, <... 'sage.rings.real_mpfr.RealNumber'>)
+  sage: parent(a) == parent(b)
+  True
+
+Sage example in ./sol/domaines.tex, line 70::
+
+  sage: Reals(100)(a)-1/10
+  0.00000000000000000000000000000
+  sage: Reals(100)(b)-1/10
+  5.5511151231257629805955278152e-18
+
+Sage example in ./sol/domaines.tex, line 80::
+
+  sage: E = CombinatorialFreeModule(QQ, [1,2,3])
+  sage: H = Hom(E,E); H.rename('H')
+  sage: C = E.category(); C
+  Category of finite dimensional vector spaces with basis over Rational Field
+  sage: phi1 = E.module_morphism(on_basis=lambda i: E.term(i), codomain=E)
+  sage: phi2 = E.module_morphism(on_basis=lambda i: E.term(i),
+  ....:                          triangular='lower', codomain=E)
+  sage: phi3 = E.module_morphism(diagonal=lambda i: 1, codomain=E,
+  ....:                          category=C)
+  sage: phi1.parent() == phi2.parent() == phi3.parent() == H
+  True
+  sage: type(phi1)
+  <class 'sage.modules.with_basis.morphism.ModuleMorphismByLinearity_with_category'>
+  sage: type(phi2)
+  <class 'sage.modules.with_basis.morphism.TriangularModuleMorphismByLinearity_with_category'>
+  sage: type(phi3)
+  <class 'sage.modules.with_basis.morphism.DiagonalModuleMorphism_with_category'>
+"""

sage/tests/books/computational-mathematics-with-sagemath/sol/float_doctest.py

@@ -0,0 +1,142 @@
+# sage_setup: distribution = sagemath-repl
+# sage.doctest: needs sage.symbolic
+"""
+This file (./sol/float_doctest.sage) was *autogenerated* from ./sol/float.tex,
+with sagetex.sty version 2011/05/27 v2.3.1.
+It contains the contents of all the sageexample environments from this file.
+You should be able to doctest this file with:
+sage -t ./sol/float_doctest.sage
+It is always safe to delete this file; it is not used in typesetting your
+document.
+
+Sage example in ./sol/float.tex, line 17::
+
+  sage: R100=RealField(100)
+  sage: x=R100(10^30)
+  sage: x>2^99
+  True
+  sage: x<2^100
+  True
+
+Sage example in ./sol/float.tex, line 48::
+
+  sage: e=2^100
+  sage: s1=10^30
+  sage: significand=[]
+  sage: nbdigits=0 # number of significant digits
+  sage: while s1>0:
+  ....:    e/=2
+  ....:    if e<=s1:
+  ....:        significand.append(1)
+  ....:        s1-=e
+  ....:    else:
+  ....:        significand.append(0)
+  ....:    nbdigits+=1
+  sage: print(significand)
+  [1, 1, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0,
+   1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0,
+   0, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1, 0, 1, 1, 0, 1,
+   1, 1, 1, 0, 1, 0, 1, 0, 0, 1]
+  sage: print("number of significant digits: " + str(nbdigits))
+  number of significant digits: 70
+
+Sage example in ./sol/float.tex, line 86::
+
+  sage: R100=RealField(100)
+  sage: x=R100(10^30)
+  sage: s,m,e = x.sign_mantissa_exponent()
+  sage: s,m,e
+  (1, 1000000000000000000000000000000, 0)
+
+Sage example in ./sol/float.tex, line 126::
+
+  sage: var("u0 u1 u2 alpha beta gamma n")
+  (u0, u1, u2, alpha, beta, gamma, n)
+  sage: recurrence = lambda a,b:  111-1130/a+3000/(a*b)
+  sage: gener1 = lambda n: (alpha*100^n+beta*6^n+gamma*5^n)
+  sage: solGen =  lambda n: gener1(n+1)/gener1(n)
+
+Sage example in ./sol/float.tex, line 137::
+
+  sage: u2 = recurrence(u1,u0)
+  sage: s = [u2==solGen(2),u1==solGen(1),u0==solGen(0)]
+  sage: t = [s[i].substitute(u0=2,u1=-4) for i in range(0,3)]
+
+Sage example in ./sol/float.tex, line 143::
+
+  sage: solve(t,alpha,beta,gamma)
+  [[alpha == 0, beta == -3/4*r1, gamma == r1]]
+
+Sage example in ./sol/float.tex, line 157::
+
+  sage: alpha=0
+  sage: beta = -3/4*gamma
+  sage: final=solGen(n)-recurrence(solGen(n-1),solGen(n-2))
+  sage: final.simplify_full()
+  0
+
+Sage example in ./sol/float.tex, line 179::
+
+  sage: def recur(x1, x0):
+  ....:     return 111 - 1130/x1 + 3000/(x0*x1)
+
+Sage example in ./sol/float.tex, line 190::
+
+  sage: u0 = 2.
+  sage: u1 = -4.
+  sage: for i in range(1,25):
+  ....:     x = recur(u1,u0)
+  ....:     print((i, x))
+  ....:     u0 = u1
+  ....:     u1 = x
+  (1, 18.5000000000000)
+  (2, 9.37837837837838)
+  (3, 7.80115273775217)
+  (4, 7.15441448097533)
+  (5, 6.80678473692481)
+  (6, 6.59263276872179)
+  ..................
+  (23, 99.9999986592167)
+  (24, 99.9999999193218)
+
+Sage example in ./sol/float.tex, line 229::
+
+  sage: var("x")
+  x
+  sage: solve(x==recurrence(x,x),x)
+  [x == 100, x == 5, x == 6]
+
+Sage example in ./sol/float.tex, line 256::
+
+  sage: RL = RealField(5000)
+  sage: u0 = RL(2)
+  sage: u1 = RL(-4)
+  sage: for i in range(1,2500):
+  ....:     x = recur(u1,u0)
+  ....:     u0 = u1
+  ....:     u1= x
+  sage: x
+  100.00000000000000000000000000000000000000000000000000000...
+
+Sage example in ./sol/float.tex, line 281::
+
+  sage: u0 = 2
+  sage: u1 = -4
+  sage: for i in range(1,2500):
+  ....:     x = recur(u1,u0)
+  ....:     u0 = u1
+  ....:     u1 = x
+  sage: float(x)
+  6.0
+
+Sage example in ./sol/float.tex, line 325::
+
+  sage: f = lambda x: x^2
+  sage: g = lambda x: x*x
+  sage: sage.rings.real_mpfi.printing_style = 'brackets'
+  sage: I = RIF(-1,1)
+  sage: f(I)
+  [0.0000000000000000 .. 1.0000000000000000]
+  sage: g(I)
+  [-1.0000000000000000 .. 1.0000000000000000]
+"""

sage/tests/books/computational-mathematics-with-sagemath/sol/graphique_doctest.py

@@ -0,0 +1,106 @@
+# sage_setup: distribution = sagemath-repl
+# sage.doctest: needs sage.symbolic
+"""
+This file (./sol/graphique_doctest.sage) was *autogenerated* from ./sol/graphique.tex,
+with sagetex.sty version 2011/05/27 v2.3.1.
+It contains the contents of all the sageexample environments from this file.
+You should be able to doctest this file with:
+sage -t ./sol/graphique_doctest.sage
+It is always safe to delete this file; it is not used in typesetting your
+document.
+
+Sage example in ./sol/graphique.tex, line 3::
+
+  sage: reset()
+
+Sage example in ./sol/graphique.tex, line 10::
+
+  sage: t = var('t'); liste = [a + cos(t) for a in srange(0, 2, 0.1)]
+  sage: g = polar_plot(liste, (t, 0, 2 * pi)); g.show(aspect_ratio = 1)
+
+Sage example in ./sol/graphique.tex, line 38::
+
+  sage: f = lambda x: abs(x**2 - 1/4)
+  sage: def liste_pts(u0, n):
+  ....:     u = u0; liste = [[u0,0]]
+  ....:     for k in range(n):
+  ....:         v, u = u, f(u)
+  ....:         liste.extend([[v,u], [u,u]])
+  ....:     return(liste)
+  sage: g = line(liste_pts(1.1, 8), rgbcolor = (.9,0,0))
+  sage: g += line(liste_pts(-.4, 8), rgbcolor = (.01,0,0))
+  sage: g += line(liste_pts(1.3, 3), rgbcolor = (.5,0,0))
+  sage: g += plot(f, -1, 3, rgbcolor = 'blue')
+  sage: g += plot(x, -1, 3, rgbcolor = 'green')
+  sage: g.show(aspect_ratio = 1, ymin = -.2, ymax = 3)
+
+Sage example in ./sol/graphique.tex, line 78::
+
+  sage: x = var('x'); y = function('y')
+  sage: DE = x^2 * diff(y(x), x) - y(x) == 0
+  sage: desolve(DE, y(x))
+  _C*e^(-1/x)
+  sage: g = plot([c*e^(-1/x) for c in srange(-8, 8, 0.4)], (x, -3, 3))
+  sage: y = var('y')
+  sage: g += plot_vector_field((x^2, y), (x,-3,3), (y,-5,5))
+  sage: g.show()  # not tested, known bug, see :issue:`32657`
+
+Sage example in ./sol/graphique.tex, line 124::
+
+  sage: from sage.calculus.desolvers import desolve_system_rk4
+  sage: f = lambda x, y: [a*x-b*x*y,-c*y+d*b*x*y]
+  sage: x, y, t = var('x, y, t')
+  sage: a, b, c, d = 1., 0.1, 1.5, 0.75
+  sage: P = desolve_system_rk4(f(x,y), [x,y],\
+  ....:        ics=[0,10,5], ivar=t, end_points=15)
+  sage: Ql = [[i,j] for i,j,k in P]; p = line(Ql, color='red')
+  sage: p += text("Rabbits", (12,37), fontsize=10, color='red')
+  sage: Qr = [[i,k] for i,j,k in P]; p += line(Qr, color='blue')
+  sage: p += text("Foxes", (12,7), fontsize=10, color='blue')
+  sage: p.axes_labels(["time", "population"])
+  sage: p.show(gridlines = True)
+
+Sage example in ./sol/graphique.tex, line 154::
+
+  sage: n = 10;  L = srange(6, 18, 12 / n); R = srange(3, 9, 6 / n)
+  sage: def g(x, y): v = vector(f(x, y)); return v / v.norm()
+  sage: q = plot_vector_field(g(x, y), (x, 0, 60), (y, 0, 36))
+  sage: for j in range(n):
+  ....:     P = desolve_system_rk4(f(x,y), [x,y],
+  ....:            ics=[0,L[j],R[j]], ivar=t, end_points=15)
+  ....:     Q = [[j,k] for i,j,k in P]
+  ....:     q += line(Q, color=hue(.8-j/(2*n)))
+  sage: q.axes_labels(["rabbits", "foxes"]); q.show()
+
+Sage example in ./sol/graphique.tex, line 185::
+
+  sage: from scipy import integrate
+  sage: def dX_dt(X, t=0): return [X[1], 0.5*X[1] - X[0] - X[1]^3]
+  sage: t = srange(0, 40, 0.01);  x0 = srange(-2, 2, 0.1); y0 = 2.5
+  sage: CI = [[i, y0] for i in x0] + [[i, -y0] for i in x0]
+  sage: def g(x, y): v = vector(dX_dt([x, y])); return v / v.norm()
+  sage: x, y = var('x, y'); n = len(CI)
+  sage: q = plot_vector_field(g(x, y), (x, -3, 3), (y, -y0, y0))
+  sage: for j in range(n):  # long time
+  ....:     X = integrate.odeint(dX_dt, CI[j], t)
+  ....:     q += line(X, color=(1.7*j/(4*n),1.5*j/(4*n),1-3*j/(8*n)))
+  sage: X = integrate.odeint(dX_dt, [0.01,0], t)
+  sage: q += line(X, color = 'red'); q.show()
+
+Sage example in ./sol/graphique.tex, line 234::
+
+  sage: from scipy import integrate
+  sage: t = srange(0, 40, 0.2)
+  sage: n = 35; CI_cart = [[4, .2 * i] for i in range(n)]
+  sage: CI = list(map(lambda x: [sqrt(x[0]^2+x[1]^2),
+  ....:      pi - arctan(x[1]/x[0])], CI_cart))
+  sage: for alpha in [0.1, 0.5, 1, 1.25]:                      # long time
+  ....:     dX_dt = lambda X, t=0: [cos(X[1])*(1-1/X[0]^2),
+  ....:             -sin(X[1]) * (1/X[0]+1/X[0]^3) + 2*alpha/X[0]^2]
+  ....:     q = circle((0, 0), 1, fill=True, rgbcolor='purple')
+  ....:     for j in range(n):
+  ....:         X = integrate.odeint(dX_dt, CI[j], t)
+  ....:         Y = [[u[0]*cos(u[1]), u[0]*sin(u[1])] for u in X]
+  ....:         q += line(Y, xmin = -4, xmax = 4, color='blue')
+  ....:     q.show(aspect_ratio = 1, axes = False)
+"""

sage/tests/books/computational-mathematics-with-sagemath/sol/graphtheory_doctest.py

@@ -0,0 +1,54 @@
+# sage_setup: distribution = sagemath-repl
+# sage.doctest: needs sage.graphs
+"""
+This file (./sol/graphtheory_doctest.sage) was *autogenerated* from ./sol/graphtheory.tex,
+with sagetex.sty version 2011/05/27 v2.3.1.
+It contains the contents of all the sageexample environments from this file.
+You should be able to doctest this file with:
+sage -t ./sol/graphtheory_doctest.sage
+It is always safe to delete this file; it is not used in typesetting your
+document.
+
+Sage example in ./sol/graphtheory.tex, line 5::
+
+  sage: def circulant(n, d):
+  ....:    g = Graph(n)
+  ....:    for u in range(n):
+  ....:        for c in range(d):
+  ....:            g.add_edge(u,(u+c)%n)
+  ....:    return g
+
+Sage example in ./sol/graphtheory.tex, line 19::
+
+  sage: def kneser(n, k):
+  ....:    g = Graph()
+  ....:    g.add_vertices(Subsets(n,k))
+  ....:    for u in g:
+  ....:        for v in g:
+  ....:            if not u & v:
+  ....:                g.add_edge(u,v)
+  ....:    return g
+
+Sage example in ./sol/graphtheory.tex, line 33::
+
+  sage: def kneser(n, k):
+  ....:    g = Graph()
+  ....:    sommets = Set(range(n))
+  ....:    g.add_vertices(Subsets(sommets,k))
+  ....:    for u in g:
+  ....:        for v in Subsets(sommets - u,k):
+  ....:            g.add_edge(u,v)
+  ....:    return g
+
+Sage example in ./sol/graphtheory.tex, line 59::
+
+  sage: g = graphs.PetersenGraph()
+  sage: def optimal_order(g):
+  ....:    order = []
+  ....:    for color_class in sorted(g.coloring()):
+  ....:        for v in color_class:
+  ....:            order.append(v)
+  ....:    return order
+  sage: optimal_order(g)
+  [0, 2, 6, 1, 3, 5, 9, 4, 7, 8]
+"""

sage/tests/books/computational-mathematics-with-sagemath/sol/integration_doctest.py

@@ -0,0 +1,57 @@
+# sage_setup: distribution = sagemath-repl
+# sage.doctest: needs sage.symbolic
+"""
+This file (./sol/integration_doctest.sage) was *autogenerated* from ./sol/integration.tex,
+with sagetex.sty version 2011/05/27 v2.3.1.
+It contains the contents of all the sageexample environments from this file.
+You should be able to doctest this file with:
+sage -t ./sol/integration_doctest.sage
+It is always safe to delete this file; it is not used in typesetting your
+document.
+
+Sage example in ./sol/integration.tex, line 22::
+
+  sage: x = var('x')
+  sage: def NCRule(n):
+  ....:     P = prod([x - j for j in range(n)])
+  ....:     return [integrate(P / (x-i), x, 0, n-1) \
+  ....:             / (P/(x-i)).subs(x=i) for i in range(n)]
+
+Sage example in ./sol/integration.tex, line 35::
+
+  sage: def QuadNC(f, a, b, n):
+  ....:     W = NCRule(n)
+  ....:     ret = 0
+  ....:     for i in range(n):
+  ....:         ret += f(a + (b-a)/(n-1)*i) * W[i]
+  ....:     return (b-a)/(n-1)*ret
+
+Sage example in ./sol/integration.tex, line 49::
+
+  sage: QuadNC(lambda u: 1, 0, 1, 12)
+  1
+  sage: N(QuadNC(sin, 0, pi, 10))
+  1.99999989482634
+
+Sage example in ./sol/integration.tex, line 59::
+
+  sage: numerical_integral(x * log(1+x), 0, 1)
+  (0.25, 2.7755575615628914e-15)
+  sage: N(QuadNC(lambda x: x * log(1+x), 0, 1, 19)) # abs tol 2e-14
+  0.250000000000001
+  sage: numerical_integral(sqrt(1-x^2), 0, 1) # abs tol 2e-14
+  (0.7853981677264822, 9.042725224536535e-07)
+  sage: N(pi/4)
+  0.785398163397448
+  sage: N(QuadNC(lambda x: sqrt(1-x^2), 0, 1, 20))
+  0.784586419900198
+
+Sage example in ./sol/integration.tex, line 74::
+
+  sage: [N(QuadNC(lambda x: x * log(1+x), 0, 1, n) - 1/4) # abs tol 3e-16
+  ....:  for n in [2, 8, 16]]
+  [0.0965735902799726, 1.17408932933522e-7, 2.13449050101566e-13]
+  sage: [N(QuadNC(lambda x: sqrt(1-x^2), 0, 1, n) - pi/4)    # abs tol 1e-13
+  ....:  for n in [2, 8, 16]]
+  [-0.285398163397448, -0.00524656673640445, -0.00125482109302663]
+"""

sage/tests/books/computational-mathematics-with-sagemath/sol/linalg_doctest.py

@@ -0,0 +1,57 @@
+# sage_setup: distribution = sagemath-repl
+# sage.doctest: needs numpy sage.libs.pari sage.modules
+"""
+This file (./sol/linalg_doctest.sage) was *autogenerated* from ./sol/linalg.tex,
+with sagetex.sty version 2011/05/27 v2.3.1.
+It contains the contents of all the sageexample environments from this file.
+You should be able to doctest this file with:
+sage -t ./sol/linalg_doctest.sage
+It is always safe to delete this file; it is not used in typesetting your
+document.
+
+Sage example in ./sol/linalg.tex, line 89::
+
+  sage: A = matrix(GF(7),[[0,0,3,0,0],[1,0,6,0,0],[0,1,5,0,0],
+  ....:                   [0,0,0,0,5],[0,0,0,1,5]])
+  sage: P = A.minpoly(); P
+  x^5 + 4*x^4 + 3*x^2 + 3*x + 1
+  sage: P.factor()
+  (x^2 + 2*x + 2) * (x^3 + 2*x^2 + x + 4)
+
+Sage example in ./sol/linalg.tex, line 100::
+
+  sage: e1 = identity_matrix(GF(7),5)[0]
+  sage: e4 = identity_matrix(GF(7),5)[3]
+  sage: A.transpose().maxspin(e1)
+  [(1, 0, 0, 0, 0), (0, 1, 0, 0, 0), (0, 0, 1, 0, 0)]
+  sage: A.transpose().maxspin(e4)
+  [(0, 0, 0, 1, 0), (0, 0, 0, 0, 1)]
+  sage: A.transpose().maxspin(e1 + e4)
+  [(1, 0, 0, 1, 0), (0, 1, 0, 0, 1), (0, 0, 1, 5, 5),
+  (3, 6, 5, 4, 2), (1, 5, 3, 3, 0)]
+
+Sage example in ./sol/linalg.tex, line 168::
+
+  sage: def Similar(A, B):
+  ....:     F1, U1 = A.frobenius_form(2)
+  ....:     F2, U2 = B.frobenius_form(2)
+  ....:     if F1 == F2:
+  ....:         return True, ~U2*U1
+  ....:     else:
+  ....:         return False, F1 - F2
+  sage: B = matrix(ZZ, [[0,1,4,0,4],[4,-2,0,-4,-2],[0,0,0,2,1],
+  ....:                 [-4,2,2,0,-1],[-4,-2,1,2,0]])
+  sage: U = matrix(ZZ, [[3,3,-9,-14,40],[-1,-2,4,2,1],[2,4,-7,-1,-13],
+  ....:                 [-1,0,1,4,-15],[-4,-13,26,8,30]])
+  sage: A = (U^-1 * B * U).change_ring(ZZ)
+  sage: ok, V = Similar(A, B); ok
+  True
+  sage: V
+  [                 1    2824643/1601680   -6818729/1601680 -43439399/11211760  73108601/11211760]
+  [                 0      342591/320336     -695773/320336  -2360063/11211760  -10291875/2242352]
+  [                 0     -367393/640672      673091/640672    -888723/4484704   15889341/4484704]
+  [                 0     661457/3203360    -565971/3203360  13485411/22423520 -69159661/22423520]
+  [                 0   -4846439/3203360    7915157/3203360 -32420037/22423520 285914347/22423520]
+  sage: ok, V = Similar(2*A, B); ok
+  False
+"""