解方程ax3+bx2+cx+d=0

honboo

解方程ax<SUP>3</SUP>+bx<SUP>2</SUP>+cx+d=0


如何用LISP编一个解1元3次方程的，只要求要它的正数根


精确到小数点后4位！

oyxx1023

就把X的通用求解公式做成LISP就行了啊.应该不麻烦的.

LONGXIN

三元的可能有些麻烦

leeyuan

相当于做函数y=ax3+bx2+cx+d与x轴的交点，用迭代法很容易做的

honboo

帮我做一个程序出来吧!

meflying

第一次写这种类型的程序,精度控制也不知道对不对,见笑了:

2. (defun c:test( / a b c d step dir m x y)
3.    (setq a (getreal "输入a:"))
4.    (setq b (getreal "输入b:"))
5.    (setq c (getreal "输入c:"))
6.    (setq d (getreal "输入d:"))
7.    (setq step 10)
8.    (setq dir 1)
9.    (setq m t)
10.    (setq i 0 x -100)
11.    (while m
12.        (setq x (+ x (* i dir step)))
13.        (setq y (+ (* a (expt x 3)) (* b (expt x 2)) (* c x) d))
14.        (if (> (* y dir) 0)
15.            (progn
16.   (if (< step 0.00001) (setq m nil))
17.   (setq dir (* dir -1))
18.   (setq step (/ step 2.0))
19.            )
20.        )      
21.        (setq i (1+ i))
22.    )
23.    (princ x)
24.    (princ)
25. )

(defun c:test( / a b c d step dir m x y)    (setq a (getreal "输入a:"))    (setq b (getreal "输入b:"))    (setq c (getreal "输入c:"))    (setq d (getreal "输入d:"))    (setq step 10)    (setq dir 1)    (setq m t)    (setq i 0 x -100)    (while m        (setq x (+ x (* i dir step)))        (setq y (+ (* a (expt x 3)) (* b (expt x 2)) (* c x) d))        (if (> (* y dir) 0)            (progn   (if (< step 0.00001) (setq m nil))   (setq dir (* dir -1))   (setq step (/ step 2.0))            )        )              (setq i (1+ i))    )    (princ x)    (princ) )

leeyuan

楼主应给出a、b、c、d的具体值，不然很有可能是该函数在x&gt;0的范围内无解

meflying

从0开始,往两边同时搜索就可以了吧

leeyuan

可是楼主只要正解，所以要楼主给出具体的参数值

honboo

a         b         c         d 没有确定值，但是它们之间的关系应该满足此方程拥有实根时判别式的要求！