ammornil написа:Може да има и по-кратко (по-хитро решение), но бакалският план е следният:...
Аз викам с груба сила да пробваме:
In [57]: var("x,m,a,k")
...: U = x**4 +(3*m-5)*x**2 +(m+1)**2
...: x0 = a; x1 = a+k; x2=a+2*k; x3=a+3*k
In [58]: [[s[m],solve(U.subs(m,s[m]))] for s in solve([U.subs(x,q) for q in (x0,x1,x2,x3)])]
...:
Out[58]:
[[-3*a**2/2 - a*sqrt(5*a**2 + 32)/2 - 1,
[{a: -x},
{a: x},
{a: -3*x/2 - sqrt(5*x**2 + 32)/2},
{a: -3*x/2 + sqrt(5*x**2 + 32)/2},
{a: 3*x/2 - sqrt(5*x**2 + 32)/2},
{a: 3*x/2 + sqrt(5*x**2 + 32)/2}]],
[-3*a**2/2 + a*sqrt(5*a**2 + 32)/2 - 1,
[{a: -x},
{a: x},
{a: -3*x/2 - sqrt(5*x**2 + 32)/2},
{a: -3*x/2 + sqrt(5*x**2 + 32)/2},
{a: 3*x/2 - sqrt(5*x**2 + 32)/2},
{a: 3*x/2 + sqrt(5*x**2 + 32)/2}]],
[-1, [0, -2*sqrt(2), 2*sqrt(2)]],
[3/5, [-2*sqrt(10)/5, 2*sqrt(10)/5]],
[3/5, [-2*sqrt(10)/5, 2*sqrt(10)/5]],
[-25, [-6*sqrt(2), -2*sqrt(2), 2*sqrt(2), 6*sqrt(2)]],
[-25, [-6*sqrt(2), -2*sqrt(2), 2*sqrt(2), 6*sqrt(2)]],
[5/19, [-6*sqrt(38)/19, -2*sqrt(38)/19, 2*sqrt(38)/19, 6*sqrt(38)/19]],
[5/19, [-6*sqrt(38)/19, -2*sqrt(38)/19, 2*sqrt(38)/19, 6*sqrt(38)/19]],
[7, [-2*sqrt(2)*I, 2*sqrt(2)*I]],
[7, [-2*sqrt(2)*I, 2*sqrt(2)*I]]]
-1 даде аритметична редица от 3 елемента, но 0 е двоен и е спорен.