Как мне избавиться от вертикальной асимптоты в Python sympy? - PullRequest
0 голосов
/ 14 ноября 2018

Я пытаюсь вычислить вторую производную от f (x) = x ** 2 / (x-1), используя пакет Sympy, но вертикальная асимптота выглядела так, как будто x приближался к 1. Как мне от этого избавиться? Вот мой код:

import sympy as sp
import numpy as np
import pylab as pyl 

t = sp.symbols('t')
def fp(t):
    return t**2 / (t -1)
f = sp.lambdify(t, fp(t), modules=['numpy'])

dp = sp.diff(fp(t))
d2p = sp.diff(dp)
d2 = sp.lambdify(t, d2p, modules=['numpy'])

x = np.linspace(-7.5, 10, 10**3)

pyl.plot(x, d2(x), 'b-')

pyl.ylim(-6, 10)
pyl.grid()

pyl.show()

Вертикальная асимптота при x = 1

1 Ответ

0 голосов
/ 14 ноября 2018

Я бы попробовал вторую производную в том же вызове:

d2p = sp.diff(t**2 / (t-1), t, 2)

и посмотрю, имеет ли это значение.Вы также можете оценить d2p и посмотреть, каково значение d2p(1).

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...