Как найти неизвестную в интеграле?

Question

[artshelom]

Использую библиотеку sympy и scipy, для нахождения интеграла и неизвестной через метод nsolve
Неизвестная C (если С подставить и убрать nsolve, все считает):

C = Symbol('C')
      print(nsolve(integrate.quad(lambda p:(((p**2-1)*(p-f))*(2*C**2)).subs('f', self.f).subs('k', self.k), math.sqrt((((self.h * self.f * self.gamma) / C) + 1)), self.f+math.sqrt(1+self.f**2))), C, 1)

А ошибка:
can't convert expression to float
Как сделать поиск неизвестной в интеграле через подбор?

Я интегрирую по p, нужно чтобы (integralP[1]-integralP[0]-545.7245820872807)=0
Вот такой код:

self.k = 1.4281480067421142
        self.f = 0.7002075382097097
        self.s = 1.63044141640584
        self.h = 27.5620864142761
        C = Symbol('C')
        lambdaP=lambda p:(((p**2-1)*(p-f))*(2*C**2)).subs('f', self.f).subs('k', self.k)
        integralP=integrate.quad(lambdaP, math.sqrt((((self.h * self.f * self.gamma) / C) + 1)), self.f+math.sqrt(1+self.f**2))
        print(integralP[1]-integralP[0]-545.7245820872807)

Comments

[Andy_U]

Огрызок вашего кода совершенно нерабочий, даже если "убрать" nsolve, по куче разных причин сразу.

1) отсутствуют import'ы;
2) Не определен self;
3) Что вы подразумевали, применяя subs к числовому выражению, возвращаемому лямбдой, в результате чего и вылезает ошибка, я понять не в силах.

В общем, вы бы написали человеческим языком, чего хотели посчитать.

[Andy_U]

Так, теперь стало немного понятнее, хотя код, который вы привели, все еще полностью нерабочий (вы, кстати еще gamma не определили). Объясните, пожалуйста, вот это свое желание:

Я интегрирую по p, нужно чтобы (integralP[1]-integralP[0]-545.7245820872807)=0

Вы читали документацию к функции scipy.integhrate.quad? Какой смысл искать величину интеграла минус ошибка интегрирования? 0 и 1 - это не пределы интегрирования - они явно заданы вторым и третьим аргументом quad.

[Andy_U]

Вдогонку - нафига quad использовать? Интеграл берется аналитически. А 2*с**2 вообще выносится за скобки.

[artshelom]

Andy_U, мне надо понять как найти неизвестную величину в интеграле. Когда я пытаюсь найти неизвестную в интеграле, то получаю ошибку выше.
Так как найти неизвестную в интеграле??
quat, для нахождения определенного интеграла, а как правильно делать?

[Andy_U]

artshelom, Слона надо есть по частям.

1) Что значит "найти неизвестную величину в интеграле"? Вы можете написать уравнение, которое хотите решить? Хоть картинкой, хот в формате latex. Ваши слова никому непонятны.

2) Если компьютер неизвестную величину "не знает" - интеграл любым численным методом не посчитать. scipy.integrate.quad - именно такой метод. Т.е. подинтегральная функция должна принимать числовой аргумент (или даже несколько - см. примеры к quad) и возвращать число, и пределы интегрирования должны быть числовыми. Такие фокусы с sympy, какие вы пытаетесь использовать - недопустимы. Это понятно? Попробуйте на простых примерах.

3) Вот это (integralP[1]-integralP[0]) - полная фигня. См.выше.

3) Чему равна self. gamma? И вообще, если используется self - речь идет о фрагменте метода какого-то класса. А телепаты в отпуске.

Answer 1

[MechanicZelenyy]

Мне лень было разбирать ваш крокодил, так что попробую угадать --- может надо сделать лямбда выражение с C в качестве аргумента?

Comments

[artshelom]

Изменил, не помогло. Если лямду делаю С