Python 怎么解方程啊？

V2EX = way to explore

V2EX 是一个关于分享和探索的地方

现在注册

已注册用户请登录

推荐学习书目

› Learn Python the Hard Way

Python Sites

› PyPI - Python Package Index

› http://diveintopython.org/toc/index.html

› Pocoo

值得关注的项目

› PyPy

› Celery

› Jinja2

› Read the Docs

› gevent

› pyenv

› virtualenv

› Stackless Python

› Beautiful Soup

› 结巴中文分词

› Green Unicorn

› Sentry

› Shovel

› Pyflakes

› pytest

Python 编程

› pep8 Checker

Styles

› PEP 8

› Google Python Style Guide

› Code Style from The Hitchhiker's Guide

这是一个创建于 2820 天前的主题，其中的信息可能已经有所发展或是发生改变。

用的 sympy 中的 solve，这个库是不是只能解有解的方程？

p1 = 0.8
p2 = 1
r1 = 500
r2 = 3500
x = symbols('x')
print('in')
s=solve(p1*exp(-x*r2)-p2*exp(-x*r1)+p2-p1, x)
print(s)

一直输出不了 s。。。像这种方程用 python 怎么解？

方程

Python

exp

17 条回复 • 2018-04-24 06:00:27 +08:00

est

2018 年 4 月 23 日

又来一个廖学峰课后作业系列？

flyzy2005

2018 年 4 月 23 日 via iPhone

@est。。不是 = = 要算个方程

l1093178

2018 年 4 月 23 日

WolframAlpha 了解一下？

flyzy2005

2018 年 4 月 23 日 via iPhone

@l1093178 谢谢算出来了

princelai

2018 年 4 月 23 日 via Android

前两天不是刚问完，sympy

flyzy2005

2018 年 4 月 23 日 via iPhone

@princelai 我用了那个没解算不出来好像换一个 x-2 很快算出来

asdf123101

2018 年 4 月 23 日 via Android

自己写个二分法？知道导数可以用牛顿。这种方程自己写一个很快的。

sciooga

2018 年 4 月 23 日

这样行不行？

import random

while True:
x = random.random()
然后将 x 带入方程两边如果相等就 break

reself

2018 年 4 月 23 日 via Android

@sciooga 哈哈猴子算法

flyzy2005

2018 年 4 月 23 日 via iPhone

@sciooga 😂😂😂 用了那个 wolframalpha

xxm459259

2018 年 4 月 23 日

sympy 啊。。。

flyzy2005

2018 年 4 月 23 日 via iPhone

@pandachow 我就用的这个库啊……你看我代码卡在 solve 算不出来

idlerlestat

2018 年 4 月 23 日 via iPhone

sympy 是给解析解，你这个方程没有解析解啊。

flyzy2005

2018 年 4 月 23 日 via iPhone

@idlerlestat 嗯。。所以我来问应该用什么

idlerlestat

2018 年 4 月 23 日

@flyzy2005 numpy 和 scipy.optimize.fsolve 就行

inflationaaron

2018 年 4 月 23 日 via iPhone

这种需求用 Mathematica 啊

20150517

2018 年 4 月 24 日 via Android

不是我说，这种最好用 tensorflow，做个 variable 然后用模型能帮你算出系数