Mercurialsympy / file revision
summary | shortlog | changelog | tags | files | changeset | file | revisions | annotate | diff | raw
sympy/series/tests/test_demidovich.py
author Vinzent Steinberg <Vinzent.Steinberg@gmail.com>
Thu Nov 20 02:14:29 2008 +0100 (2 weeks ago)
changeset 1916 42fa25e25e6e
parent 16960d5613d7b037
permissions -rw-r--r--
add a comment why a "working" test is xfailed 
     1 from sympy import limit, Symbol, oo, sqrt, Rational, log, exp, cos, sin, tan, \
     2     pi, asin, together
     3 
     4 """
     5 (*) in problem number means that the number is relative to the book "Anti-demidovich,
     6 problemas resueltos, Ed. URSS"
     7 
     8 """
     9 
    10 x = Symbol("x", real=True, unbounded=True)
    11 
    12 def test_leadterm():
    13     assert (3+2*x**(log(3)/log(2)-1)).leadterm(x)==(3,0)
    14 
    15 def sqrt3(x):
    16     return x**Rational(1,3)
    17 
    18 def sqrt4(x):
    19     return x**Rational(1,4)
    20 
    21 def test_Limits_simple_0():
    22     assert limit((2**(x+1)+3**(x+1))/(2**x+3**x),x,oo)==3  #175
    23 
    24 def test_Limits_simple_1():
    25     assert limit((x+1)*(x+2)*(x+3)/x**3,x, oo)==1  #172
    26     assert limit(sqrt(x+1)-sqrt(x),x,oo)==0  #179
    27     assert limit((2*x-3)*(3*x+5)*(4*x-6)/(3*x**3+x-1),x,oo)==8  #Primjer 1
    28     assert limit(x/sqrt3(x**3+10),x,oo)==1  #Primjer 2
    29     assert limit((x+1)**2/(x**2+1),x,oo)==1  #181
    30 
    31 def test_Limits_simple_2():
    32     assert limit(1000*x/(x**2-1),x,oo)==0  #182
    33     assert limit((x**2-5*x+1)/(3*x+7),x,oo)==oo  #183
    34     assert limit((2*x**2-x+3)/(x**3-8*x+5),x,oo)==0  #184
    35     assert limit((2*x**2-3*x-4)/sqrt(x**4+1),x,oo)==2  #186
    36     assert limit((2*x+3)/(x+sqrt3(x)),x,oo)==2  #187
    37     assert limit(x**2/(10+x*sqrt(x)),x,oo)==oo  #188
    38     assert limit(sqrt3(x**2+1)/(x+1),x,oo)==0  #189
    39     assert limit(sqrt(x)/sqrt(x+sqrt(x+sqrt(x))),x,oo)==1  #190
    40 
    41 def test_Limits_simple_3a():
    42     a = Symbol('a', real=True)
    43     #issue 414
    44     assert together(limit((x**2-(a+1)*x+a)/(x**3-a**3),x,a)) == \
    45             (a-1)/(3*a**2)  #196
    46 
    47 def test_Limits_simple_3b():
    48     h = Symbol("h")
    49     assert limit(((x+h)**3-x**3)/h,h,0)==3*x**2  #197
    50     assert limit((1/(1-x)-3/(1-x**3)),x,1)==-1  #198
    51     assert limit((sqrt(1+x)-1)/(sqrt3(1+x)-1),x,0)==Rational(3)/2  #Primer 4
    52     assert limit((sqrt(x)-1)/(x-1),x,1)==Rational(1)/2  #199
    53     assert limit((sqrt(x)-8)/(sqrt3(x)-4),x,64)==3  #200
    54     assert limit((sqrt3(x)-1)/(sqrt4(x)-1),x,1)==Rational(4)/3  #201
    55     assert limit((sqrt3(x**2)-2*sqrt3(x)+1)/(x-1)**2,x,1)==Rational(1)/9  #202
    56 
    57 def test_Limits_simple_4a():
    58     a = Symbol('a', real=True)
    59     assert limit((sqrt(x)-sqrt(a))/(x-a),x,a)==1/(2*sqrt(a))  #Primer 5
    60     assert limit((sqrt(x)-1)/(sqrt3(x)-1),x,1)==Rational(3)/2  #205
    61     assert limit((sqrt(1+x)-sqrt(1-x))/x,x,0)==1  #207
    62     assert limit(sqrt(x**2-5*x+6)-x,x,oo)==-Rational(5)/2  #213
    63 
    64 def test_limits_simple_4aa():
    65     assert limit(x*(sqrt(x**2+1)-x),x,oo)==Rational(1)/2  #214
    66 
    67 def test_Limits_simple_4b():
    68     #issue 412
    69     assert limit(x-sqrt3(x**3-1),x,oo)==0  #215
    70 
    71 def test_Limits_simple_4c():
    72     assert limit(log(1+exp(x))/x,x,-oo)==0  #267a
    73     assert limit(log(1+exp(x))/x,x,oo)==1  #267b
    74 
    75 def test_bounded():
    76     assert limit(sin(x)/x, x, oo) == 0 #216b
    77     assert limit(x*sin(1/x), x, 0) == 0 #227a
    78 
    79 def test_f1a():
    80     h = Symbol("h")
    81     #issue 409:
    82     assert limit((sin(2*x)/x)**(1+x),x,0) == 2 #Primer 7
    83 
    84 def test_f1a2():
    85     #issue 410:
    86     assert limit(((x-1)/(x+1))**x,x,oo) == exp(-2) #Primer 9
    87 
    88 def test_f1b():
    89     m = Symbol("m")
    90     n = Symbol("n")
    91     h = Symbol("h")
    92     a = Symbol("a")
    93     assert limit(sin(x)/x,x,2) == sin(2)/2 #216a
    94     assert limit(sin(3*x)/x,x,0) == 3 #217
    95     assert limit(sin(5*x)/sin(2*x),x,0) == Rational(5)/2 #218
    96     assert limit(sin(pi*x)/sin(3*pi*x),x,0) == Rational(1)/3 #219
    97     assert limit(x*sin(pi/x),x,oo) == pi #220
    98     assert limit((1-cos(x))/x**2,x,0) == Rational(1,2) #221
    99     assert limit(x*sin(1/x),x,oo) == 1 #227b
   100     assert limit((cos(m*x)-cos(n*x))/x**2,x,0) == ((n**2-m**2)/2) #232
   101     assert limit((tan(x)-sin(x))/x**3,x,0) == Rational(1,2) #233
   102     assert limit((x-sin(2*x))/(x+sin(3*x)),x,0) == -Rational(1,4) #237
   103     assert limit((1-sqrt(cos(x)))/x**2,x,0) == Rational(1,4) #239
   104     assert limit((sqrt(1+sin(x))-sqrt(1-sin(x)))/x,x,0) == 1 #240
   105 
   106     assert limit((1+h/x)**x,x,oo) == exp(h) #Primer 9
   107     assert limit((sin(x)-sin(a))/(x-a),x,a) == cos(a) #222, *176
   108     assert limit((cos(x)-cos(a))/(x-a),x,a) == -sin(a) #223
   109     assert limit((sin(x+h)-sin(x))/h,h,0) == cos(x) #225
   110 
   111 def test_f2a():
   112     assert limit(((x+1)/(2*x+1))**(x**2),x,oo) == 0 #Primer 8
   113 
   114 def test_f2():
   115     assert limit((sqrt(cos(x))-sqrt3(cos(x)))/(sin(x)**2),x,0) == -Rational(1,12) #*184
   116 
   117 def test_f3():
   118     a = Symbol('a', real=True)
   119     #issue 405
   120     assert limit(asin(a*x)/x, x, 0) == a
sympy