Matlab编写程序用二分法求解非线性方程的根
用二分法求方程x*x-x-1=0的正根，要求精确到小数点后四位。（matlab）
l1
计算公式
f(ak)*f(bk)<0;
bk-ak=1/2k-1*(b-a);
a1<=a2<=…
<=ak<=…,b1<=b2<=…<=bk<=…。
l2
算法分析
设f(x)∈C[a,b],且f(a)f(b)0,则a2=x1,b2=b1,否则a2=a1,b2=x1。得到[a2,b2]满足：f(a2)f(b2)=a1,b2<=b1。以[a2,b2]取代[a1,b1]，继续以上过程，直到精度达到要求为止。
l3
源程序
function
f1=fun(x);
f1=x-cos(x);
function
[x,k]=erfen(a,b,s)
%a,b为根区间，s为精度
a=0;b=1;s=1e-4;k=0;
while
abs(a-b)>s
x=(a+b)/2;
if
fun(a)*fun(x)<0
b=x;
else
a=x;
end
k=k+1;
end
x=(a+b)/2
%x为方程的解
k
%
k为计算次数
实验结果讨论和分析
本题使用二分法得到的x=0.7391，满足基本要求，题目要求精确到小数点后四位，告诉了本题二分法得应达到得精确度；计算次数为14，二分法收敛性很好，收敛速度不快。
Matlab编写程序用二分法求解非线性方程的根
用二分法求方程x*x-x-1=0的正根，要求精确到小数点后四位。（matlab）
l1 计算公式
f(ak)*f(bk)<0;
bk-ak=1/2k-1*(b-a);
a1<=a2<=… <=ak<=…,b1<=b2<=…<=bk<=…。
l2 算法分析
设f(x)∈C[a,b],且f(a)f(b)0,则a2=x1,b2=b1,否则a2=a1,b2=x1。得到[a2,b2]满足：f(a2)f(b2)=a1,b2<=b1。以[a2,b2]取代[a1,b1]，继续以上过程，直到精度达到要求为止。
l3 源程序
function f1=fun(x);
f1=x-cos(x);
function [x,k]=erfen(a,b,s) %a,b为根区间，s为精度
a=0;b=1;s=1e-4;k=0;
while abs(a-b)>s
x=(a+b)/2;
if fun(a)*fun(x)<0
b=x;
else
a=x;
end
k=k+1;
end
x=(a+b)/2 %x为方程的解
k % k为计算次数
实验结果讨论和分析
本题使用二分法得到的x=0.7391，满足基本要求，题目要求精确到小数点后四位，告诉了本题二分法得应达到得精确度；计算次数为14，二分法收敛性很好，收敛速度不快。
matlab二分法求方程的根
matlab源程序如下：function erfenfa(a,b)%a，b为区间，s=(a+b)/2;，while b-a>1e-5 if fun(a)*fun(s)>0。 a=s; elseif fun(a)*fun(s)<0function y=fun(x)二分法 即一分为二的方法。设[a，b]为R的紧区间， 逐次二分法就是造出如下的区间序列：a0=a，b0=b，且对任一自然数n，[an+1，bn+1]或者等于[an，cn]，或者等于[cn，bn]，其中cn表示[an，bn]的中点。一般地，对于函数f(x),如果存在实数c,当x=c时，若f(c)=0,那么把x=c叫做函数f(x)的零点。解方程即要求f(x)的所有零点。先找到a、b属于区间（x，y），使f(a)，f(b)异号，说明在区间(a,b)内一定有零点，然后求f[(a+b)/2],现在假设f(a)0,a<b如果f[(a+b)/2]=0，该点就是零点，如果f[(a+b)/2]<0,则在区间（(a+b)/2，b)内有零点，(a+b)/2赋给a，从①开始继续使用中点函数值判断。如果f[(a+b)/2]>0，则在区间(a,(a+b)/2)内有零点，(a+b)/2赋给b，从①开始继续使用中点函数值判断。通过每次把f(x)的零点所在小区间收缩一半的方法，使区间的两个端点逐步迫近函数的零点，以求得零点的近似值，这种方法叫做二分法。
如何用MATLAB 求解黎卡提代数方程？
Unbalanced or misused parentheses or brackets.是提示错误括号不对称，缺了或者多了．
lqr函数的格式有这几种
[K,S,e] = lqr(SYS,Q,R)
[K,S,e] = lqr(SYS,Q,R,N)
[K,S,e] = lqr(A,B,Q,R,N)
楼主[K,P,E]=lqr[A,B,Q,R,N]，应该把后面改成小括号
[K,P,E]=lqr（A,B,Q,R,N）
我有疑问哦，楼主你求最优反馈增益矩阵K
自接用
K=lqr(A,B,Q,R) 就可以算出来K啦。
K =
-0.0008 0.0187 0.0065
怎样使用matlab解下面的代数方程？急急。。。。。
这个答案不就是y=f(r)的形式吗？只是比较长而已。。。
可以这样解决：
syms a b c d e;
solve('2*b^2=a^2+c^2+d^2-d*((4*r^2-y^2)^(1/2)*cos(e)+y*sin(e))-2*c*d*((r^2-d^2*(cos(e))^2)^(1/2)*cos(e)+d*cos(e)*sin(e))/r+c*((4*r^2-y^2)^(1/2)* (r^2-d^2*(cos(e))^2)^(1/2)+y*d*cos(e))')
Warning: Warning, solutions may have been lost
ans =
-(-2*b^2*r+a^2*r+c^2*r+d^2*r+d*(-2*c^2*r^2*d*cos(e)+c^3*r*(r^2-d^2*cos(e)^2)^(1/2)+2*b^2*r*d*cos(e)-d^3*r*cos(e)+a^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)+d^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)-a^2*r*d*cos(e)-2*b^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)-c^2*r*d*cos(e)-2*d^3*c*sin(e)^3+2*d^3*c*sin(e)+2*d^2*c*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-2*c^2*d^2*cos(e)*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+2*cos(e)^3*c^2*d^3-(-4*d^6*c^2*cos(e)^6-d^2*c^4*r^2-d^2*a^4*r^2+4*d^4*b^2*r^2-4*d^2*b^4*r^2-2*d^2*a^2*r^2*c^2-2*d^4*c^2*r^2-2*d^4*a^2*r^2+4*d^2*b^2*r^2*c^2+2*d^4*r^2*cos(e)^2*a^2+4*d^2*b^2*r^2*a^2-4*c^4*r*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+4*c^4*r*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)+8*c^4*r^2*d^3*cos(e)^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)-8*c^4*r^2*d^4*cos(e)^4+4*d^6*c^2*sin(e)^2+24*b^2*r*d^3*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+d^6*r^2*cos(e)^2+4*d^6*r*cos(e)*c*sin(e)^3-4*b^2*r^2*d^4*cos(e)^2-4*b^2*r^2*d^2*cos(e)^2*a^2-4*d^5*r*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)-4*b^2*r^2*d^2*cos(e)^2*c^2-8*b^2*r*d^4*cos(e)*c*sin(e)^3+4*b^4*r^2*d^2*cos(e)^2+8*b^2*r*d^3*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)-2*d^4*r^2*cos(e)^2*c^2+4*c^5*r*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3-4*a^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+4*a^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)+8*r^2*c^4*d^4*sin(e)^4-12*d^5*r*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+2*a^2*r^2*d^2*cos(e)^2*c^2+4*a^2*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3+a^4*r^2*d^2*cos(e)^2+4*d^5*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)+4*d^5*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3-4*d^5*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+4*a^2*r*d^4*cos(e)*c*sin(e)^3+c^4*r^2*d^2*cos(e)^2-4*a^2*r*d^3*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)+8*d^5*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^5-12*a^2*r*d^3*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*cos(e)^6*c^4*d^6-8*b^2*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3-4*c^3*r*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)+4*c^3*r*d^4*cos(e)*sin(e)^3+4*r^6*c^4*d^2*cos(e)^2+4*d^6*c^2*sin(e)^6-12*c^4*r*d^3*cos(e)^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*b^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+8*d^6*c^3*sin(e)^3*cos(e)^3-8*b^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)-32*d^6*c^3*sin(e)*cos(e)^3-8*c*d^2*cos(e)*sin(e)*b^2*r^2*a^2+12*d^6*c^2*cos(e)^2-8*r^2*c^4*d^3*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)-8*c^4*d^5*cos(e)^4*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*r^2*c^4*d^3*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+8*d^5*c^3*sin(e)^2*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-24*d^5*c^3*sin(e)^4*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*d^5*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-8*d^5*c^2*sin(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-d^6*c^2*r^2*cos(e)^2-4*r^6*c^2*d^2*cos(e)^2-c^6*r^2*d^2*cos(e)^2+4*d^4*c^2*cos(e)^4*r^2+8*d^6*c^3*cos(e)^5*sin(e)+10*r^2*c^4*d^4*cos(e)^2-4*c^4*d^6*cos(e)^4*sin(e)^2-8*c^3*d^2*cos(e)*sin(e)*b^2*r^2-8*c^3*d^4*cos(e)^3*sin(e)*r^4-8*r^6*c^3*d^2*sin(e)*cos(e)-8*c*d^4*cos(e)*sin(e)*b^2*r^2+8*d^5*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+8*c*d^2*cos(e)*sin(e)*b^4*r^2-24*c^3*d^5*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)+2*c*d^2*cos(e)*sin(e)*a^4*r^2+16*c^3*d^5*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*c*d^4*cos(e)*sin(e)*a^2*r^2-8*c^2*d^4*cos(e)^2*sin(e)^2*a^2*r-16*c^3*d^6*cos(e)*sin(e)^5+16*c^2*d^4*cos(e)^2*sin(e)^2*b^2*r+16*c^3*d^6*cos(e)*sin(e)+8*c*d^4*cos(e)^3*sin(e)*r^4+2*c*d^6*cos(e)*sin(e)*r^2+2*c^5*d^2*cos(e)*sin(e)*r^2+16*c^2*d^4*cos(e)^2*sin(e)^2*r^4+4*c^3*d^2*cos(e)*sin(e)*a^2*r^2+4*d^4*r^4*c^2*cos(e)^2-8*c^3*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*r^4-8*c^2*d^6*cos(e)^2*sin(e)^2*r-16*c*d^4*cos(e)*sin(e)*r^4-8*c^4*d^4*cos(e)^2*sin(e)^2*r-4*b^4*r^2*c^2*d^2*cos(e)^2-8*d^3*b^2*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*c*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*r^4+4*d^5*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*d^3*c^3*r*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-8*c*d^3*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)*r^4+16*c^2*d^3*cos(e)^2*(r^2-d^2*cos(e)^2)^(1/2)*r^4*sin(e)-8*d^5*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)-8*c^3*r^2*d^4*cos(e)*sin(e)^3+12*c^3*r^2*d^4*cos(e)*sin(e)-4*d^4*r^4*cos(e)^2+4*d^3*a^2*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*a^2*r^2*c^2*b^2*d^2*cos(e)^2+4*a^2*r*c^3*d^4*cos(e)^3*sin(e)-a^4*r^2*c^2*d^2*cos(e)^2+4*c^5*r*d^4*cos(e)^3*sin(e)-8*r^2*c^4*d^4+4*r^6*c^2*d^2+4*c^4*r^2*b^2*d^2*cos(e)^2-2*d^4*a^2*r^2*c^2*cos(e)^2+4*d^4*b^2*r^2*c^2*cos(e)^2-2*c^4*r^2*a^2*d^2*cos(e)^2+4*d^4*r^4+4*c^4*d^4*cos(e)^2*sin(e)^2*r^2+4*c^3*r*d^6*cos(e)^3*sin(e)-8*d^6*c^2-8*b^2*r*c^3*d^4*cos(e)^3*sin(e)-d^6*r^2)^(1/2))/(r^2*c^2-2*c*d^2*cos(e)*sin(e)-2*c*d*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+d^2)*cos(e)-2*c*d*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-2*c*d^2*cos(e)*sin(e)-c*(-2*c^2*r^2*d*cos(e)+c^3*r*(r^2-d^2*cos(e)^2)^(1/2)+2*b^2*r*d*cos(e)-d^3*r*cos(e)+a^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)+d^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)-a^2*r*d*cos(e)-2*b^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)-c^2*r*d*cos(e)-2*d^3*c*sin(e)^3+2*d^3*c*sin(e)+2*d^2*c*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-2*c^2*d^2*cos(e)*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+2*cos(e)^3*c^2*d^3-(-4*d^6*c^2*cos(e)^6-d^2*c^4*r^2-d^2*a^4*r^2+4*d^4*b^2*r^2-4*d^2*b^4*r^2-2*d^2*a^2*r^2*c^2-2*d^4*c^2*r^2-2*d^4*a^2*r^2+4*d^2*b^2*r^2*c^2+2*d^4*r^2*cos(e)^2*a^2+4*d^2*b^2*r^2*a^2-4*c^4*r*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+4*c^4*r*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)+8*c^4*r^2*d^3*cos(e)^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)-8*c^4*r^2*d^4*cos(e)^4+4*d^6*c^2*sin(e)^2+24*b^2*r*d^3*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+d^6*r^2*cos(e)^2+4*d^6*r*cos(e)*c*sin(e)^3-4*b^2*r^2*d^4*cos(e)^2-4*b^2*r^2*d^2*cos(e)^2*a^2-4*d^5*r*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)-4*b^2*r^2*d^2*cos(e)^2*c^2-8*b^2*r*d^4*cos(e)*c*sin(e)^3+4*b^4*r^2*d^2*cos(e)^2+8*b^2*r*d^3*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)-2*d^4*r^2*cos(e)^2*c^2+4*c^5*r*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3-4*a^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+4*a^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)+8*r^2*c^4*d^4*sin(e)^4-12*d^5*r*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+2*a^2*r^2*d^2*cos(e)^2*c^2+4*a^2*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3+a^4*r^2*d^2*cos(e)^2+4*d^5*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)+4*d^5*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3-4*d^5*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+4*a^2*r*d^4*cos(e)*c*sin(e)^3+c^4*r^2*d^2*cos(e)^2-4*a^2*r*d^3*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)+8*d^5*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^5-12*a^2*r*d^3*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*cos(e)^6*c^4*d^6-8*b^2*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3-4*c^3*r*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)+4*c^3*r*d^4*cos(e)*sin(e)^3+4*r^6*c^4*d^2*cos(e)^2+4*d^6*c^2*sin(e)^6-12*c^4*r*d^3*cos(e)^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*b^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+8*d^6*c^3*sin(e)^3*cos(e)^3-8*b^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)-32*d^6*c^3*sin(e)*cos(e)^3-8*c*d^2*cos(e)*sin(e)*b^2*r^2*a^2+12*d^6*c^2*cos(e)^2-8*r^2*c^4*d^3*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)-8*c^4*d^5*cos(e)^4*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*r^2*c^4*d^3*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+8*d^5*c^3*sin(e)^2*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-24*d^5*c^3*sin(e)^4*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*d^5*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-8*d^5*c^2*sin(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-d^6*c^2*r^2*cos(e)^2-4*r^6*c^2*d^2*cos(e)^2-c^6*r^2*d^2*cos(e)^2+4*d^4*c^2*cos(e)^4*r^2+8*d^6*c^3*cos(e)^5*sin(e)+10*r^2*c^4*d^4*cos(e)^2-4*c^4*d^6*cos(e)^4*sin(e)^2-8*c^3*d^2*cos(e)*sin(e)*b^2*r^2-8*c^3*d^4*cos(e)^3*sin(e)*r^4-8*r^6*c^3*d^2*sin(e)*cos(e)-8*c*d^4*cos(e)*sin(e)*b^2*r^2+8*d^5*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+8*c*d^2*cos(e)*sin(e)*b^4*r^2-24*c^3*d^5*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)+2*c*d^2*cos(e)*sin(e)*a^4*r^2+16*c^3*d^5*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*c*d^4*cos(e)*sin(e)*a^2*r^2-8*c^2*d^4*cos(e)^2*sin(e)^2*a^2*r-16*c^3*d^6*cos(e)*sin(e)^5+16*c^2*d^4*cos(e)^2*sin(e)^2*b^2*r+16*c^3*d^6*cos(e)*sin(e)+8*c*d^4*cos(e)^3*sin(e)*r^4+2*c*d^6*cos(e)*sin(e)*r^2+2*c^5*d^2*cos(e)*sin(e)*r^2+16*c^2*d^4*cos(e)^2*sin(e)^2*r^4+4*c^3*d^2*cos(e)*sin(e)*a^2*r^2+4*d^4*r^4*c^2*cos(e)^2-8*c^3*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*r^4-8*c^2*d^6*cos(e)^2*sin(e)^2*r-16*c*d^4*cos(e)*sin(e)*r^4-8*c^4*d^4*cos(e)^2*sin(e)^2*r-4*b^4*r^2*c^2*d^2*cos(e)^2-8*d^3*b^2*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*c*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*r^4+4*d^5*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*d^3*c^3*r*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-8*c*d^3*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)*r^4+16*c^2*d^3*cos(e)^2*(r^2-d^2*cos(e)^2)^(1/2)*r^4*sin(e)-8*d^5*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)-8*c^3*r^2*d^4*cos(e)*sin(e)^3+12*c^3*r^2*d^4*cos(e)*sin(e)-4*d^4*r^4*cos(e)^2+4*d^3*a^2*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*a^2*r^2*c^2*b^2*d^2*cos(e)^2+4*a^2*r*c^3*d^4*cos(e)^3*sin(e)-a^4*r^2*c^2*d^2*cos(e)^2+4*c^5*r*d^4*cos(e)^3*sin(e)-8*r^2*c^4*d^4+4*r^6*c^2*d^2+4*c^4*r^2*b^2*d^2*cos(e)^2-2*d^4*a^2*r^2*c^2*cos(e)^2+4*d^4*b^2*r^2*c^2*cos(e)^2-2*c^4*r^2*a^2*d^2*cos(e)^2+4*d^4*r^4+4*c^4*d^4*cos(e)^2*sin(e)^2*r^2+4*c^3*r*d^6*cos(e)^3*sin(e)-8*d^6*c^2-8*b^2*r*c^3*d^4*cos(e)^3*sin(e)-d^6*r^2)^(1/2))/(r^2*c^2-2*c*d^2*cos(e)*sin(e)-2*c*d*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+d^2)*(r^2-d^2*cos(e)^2)^(1/2))/d/r/(-sin(e)+c*cos(e))
-(-2*b^2*r+a^2*r+c^2*r+d^2*r+d*(-2*c^2*r^2*d*cos(e)+c^3*r*(r^2-d^2*cos(e)^2)^(1/2)+2*b^2*r*d*cos(e)-d^3*r*cos(e)+a^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)+d^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)-a^2*r*d*cos(e)-2*b^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)-c^2*r*d*cos(e)-2*d^3*c*sin(e)^3+2*d^3*c*sin(e)+2*d^2*c*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-2*c^2*d^2*cos(e)*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+2*cos(e)^3*c^2*d^3+(-4*d^6*c^2*cos(e)^6-d^2*c^4*r^2-d^2*a^4*r^2+4*d^4*b^2*r^2-4*d^2*b^4*r^2-2*d^2*a^2*r^2*c^2-2*d^4*c^2*r^2-2*d^4*a^2*r^2+4*d^2*b^2*r^2*c^2+2*d^4*r^2*cos(e)^2*a^2+4*d^2*b^2*r^2*a^2-4*c^4*r*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+4*c^4*r*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)+8*c^4*r^2*d^3*cos(e)^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)-8*c^4*r^2*d^4*cos(e)^4+4*d^6*c^2*sin(e)^2+24*b^2*r*d^3*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+d^6*r^2*cos(e)^2+4*d^6*r*cos(e)*c*sin(e)^3-4*b^2*r^2*d^4*cos(e)^2-4*b^2*r^2*d^2*cos(e)^2*a^2-4*d^5*r*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)-4*b^2*r^2*d^2*cos(e)^2*c^2-8*b^2*r*d^4*cos(e)*c*sin(e)^3+4*b^4*r^2*d^2*cos(e)^2+8*b^2*r*d^3*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)-2*d^4*r^2*cos(e)^2*c^2+4*c^5*r*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3-4*a^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+4*a^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)+8*r^2*c^4*d^4*sin(e)^4-12*d^5*r*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+2*a^2*r^2*d^2*cos(e)^2*c^2+4*a^2*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3+a^4*r^2*d^2*cos(e)^2+4*d^5*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)+4*d^5*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3-4*d^5*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+4*a^2*r*d^4*cos(e)*c*sin(e)^3+c^4*r^2*d^2*cos(e)^2-4*a^2*r*d^3*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)+8*d^5*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^5-12*a^2*r*d^3*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*cos(e)^6*c^4*d^6-8*b^2*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3-4*c^3*r*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)+4*c^3*r*d^4*cos(e)*sin(e)^3+4*r^6*c^4*d^2*cos(e)^2+4*d^6*c^2*sin(e)^6-12*c^4*r*d^3*cos(e)^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*b^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+8*d^6*c^3*sin(e)^3*cos(e)^3-8*b^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)-32*d^6*c^3*sin(e)*cos(e)^3-8*c*d^2*cos(e)*sin(e)*b^2*r^2*a^2+12*d^6*c^2*cos(e)^2-8*r^2*c^4*d^3*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)-8*c^4*d^5*cos(e)^4*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*r^2*c^4*d^3*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+8*d^5*c^3*sin(e)^2*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-24*d^5*c^3*sin(e)^4*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*d^5*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-8*d^5*c^2*sin(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-d^6*c^2*r^2*cos(e)^2-4*r^6*c^2*d^2*cos(e)^2-c^6*r^2*d^2*cos(e)^2+4*d^4*c^2*cos(e)^4*r^2+8*d^6*c^3*cos(e)^5*sin(e)+10*r^2*c^4*d^4*cos(e)^2-4*c^4*d^6*cos(e)^4*sin(e)^2-8*c^3*d^2*cos(e)*sin(e)*b^2*r^2-8*c^3*d^4*cos(e)^3*sin(e)*r^4-8*r^6*c^3*d^2*sin(e)*cos(e)-8*c*d^4*cos(e)*sin(e)*b^2*r^2+8*d^5*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+8*c*d^2*cos(e)*sin(e)*b^4*r^2-24*c^3*d^5*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)+2*c*d^2*cos(e)*sin(e)*a^4*r^2+16*c^3*d^5*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*c*d^4*cos(e)*sin(e)*a^2*r^2-8*c^2*d^4*cos(e)^2*sin(e)^2*a^2*r-16*c^3*d^6*cos(e)*sin(e)^5+16*c^2*d^4*cos(e)^2*sin(e)^2*b^2*r+16*c^3*d^6*cos(e)*sin(e)+8*c*d^4*cos(e)^3*sin(e)*r^4+2*c*d^6*cos(e)*sin(e)*r^2+2*c^5*d^2*cos(e)*sin(e)*r^2+16*c^2*d^4*cos(e)^2*sin(e)^2*r^4+4*c^3*d^2*cos(e)*sin(e)*a^2*r^2+4*d^4*r^4*c^2*cos(e)^2-8*c^3*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*r^4-8*c^2*d^6*cos(e)^2*sin(e)^2*r-16*c*d^4*cos(e)*sin(e)*r^4-8*c^4*d^4*cos(e)^2*sin(e)^2*r-4*b^4*r^2*c^2*d^2*cos(e)^2-8*d^3*b^2*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*c*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*r^4+4*d^5*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*d^3*c^3*r*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-8*c*d^3*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)*r^4+16*c^2*d^3*cos(e)^2*(r^2-d^2*cos(e)^2)^(1/2)*r^4*sin(e)-8*d^5*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)-8*c^3*r^2*d^4*cos(e)*sin(e)^3+12*c^3*r^2*d^4*cos(e)*sin(e)-4*d^4*r^4*cos(e)^2+4*d^3*a^2*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*a^2*r^2*c^2*b^2*d^2*cos(e)^2+4*a^2*r*c^3*d^4*cos(e)^3*sin(e)-a^4*r^2*c^2*d^2*cos(e)^2+4*c^5*r*d^4*cos(e)^3*sin(e)-8*r^2*c^4*d^4+4*r^6*c^2*d^2+4*c^4*r^2*b^2*d^2*cos(e)^2-2*d^4*a^2*r^2*c^2*cos(e)^2+4*d^4*b^2*r^2*c^2*cos(e)^2-2*c^4*r^2*a^2*d^2*cos(e)^2+4*d^4*r^4+4*c^4*d^4*cos(e)^2*sin(e)^2*r^2+4*c^3*r*d^6*cos(e)^3*sin(e)-8*d^6*c^2-8*b^2*r*c^3*d^4*cos(e)^3*sin(e)-d^6*r^2)^(1/2))/(r^2*c^2-2*c*d^2*cos(e)*sin(e)-2*c*d*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+d^2)*cos(e)-2*c*d*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-2*c*d^2*cos(e)*sin(e)-c*(-2*c^2*r^2*d*cos(e)+c^3*r*(r^2-d^2*cos(e)^2)^(1/2)+2*b^2*r*d*cos(e)-d^3*r*cos(e)+a^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)+d^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)-a^2*r*d*cos(e)-2*b^2*r*c*(r^2-d^2*cos(e)^2)^(1/2)-c^2*r*d*cos(e)-2*d^3*c*sin(e)^3+2*d^3*c*sin(e)+2*d^2*c*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-2*c^2*d^2*cos(e)*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+2*cos(e)^3*c^2*d^3+(-4*d^6*c^2*cos(e)^6-d^2*c^4*r^2-d^2*a^4*r^2+4*d^4*b^2*r^2-4*d^2*b^4*r^2-2*d^2*a^2*r^2*c^2-2*d^4*c^2*r^2-2*d^4*a^2*r^2+4*d^2*b^2*r^2*c^2+2*d^4*r^2*cos(e)^2*a^2+4*d^2*b^2*r^2*a^2-4*c^4*r*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+4*c^4*r*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)+8*c^4*r^2*d^3*cos(e)^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)-8*c^4*r^2*d^4*cos(e)^4+4*d^6*c^2*sin(e)^2+24*b^2*r*d^3*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+d^6*r^2*cos(e)^2+4*d^6*r*cos(e)*c*sin(e)^3-4*b^2*r^2*d^4*cos(e)^2-4*b^2*r^2*d^2*cos(e)^2*a^2-4*d^5*r*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)-4*b^2*r^2*d^2*cos(e)^2*c^2-8*b^2*r*d^4*cos(e)*c*sin(e)^3+4*b^4*r^2*d^2*cos(e)^2+8*b^2*r*d^3*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)-2*d^4*r^2*cos(e)^2*c^2+4*c^5*r*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3-4*a^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+4*a^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)+8*r^2*c^4*d^4*sin(e)^4-12*d^5*r*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+2*a^2*r^2*d^2*cos(e)^2*c^2+4*a^2*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3+a^4*r^2*d^2*cos(e)^2+4*d^5*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)+4*d^5*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3-4*d^5*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+4*a^2*r*d^4*cos(e)*c*sin(e)^3+c^4*r^2*d^2*cos(e)^2-4*a^2*r*d^3*cos(e)^3*c*(r^2-d^2*cos(e)^2)^(1/2)+8*d^5*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^5-12*a^2*r*d^3*cos(e)^2*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*cos(e)^6*c^4*d^6-8*b^2*r*c^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^3*d^3-4*c^3*r*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)+4*c^3*r*d^4*cos(e)*sin(e)^3+4*r^6*c^4*d^2*cos(e)^2+4*d^6*c^2*sin(e)^6-12*c^4*r*d^3*cos(e)^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*b^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)^3+8*d^6*c^3*sin(e)^3*cos(e)^3-8*b^2*r*c^2*(r^2-d^2*cos(e)^2)^(1/2)*d^3*sin(e)-32*d^6*c^3*sin(e)*cos(e)^3-8*c*d^2*cos(e)*sin(e)*b^2*r^2*a^2+12*d^6*c^2*cos(e)^2-8*r^2*c^4*d^3*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)-8*c^4*d^5*cos(e)^4*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*r^2*c^4*d^3*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+8*d^5*c^3*sin(e)^2*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-24*d^5*c^3*sin(e)^4*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*d^5*c^2*sin(e)*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-8*d^5*c^2*sin(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*cos(e)^2-d^6*c^2*r^2*cos(e)^2-4*r^6*c^2*d^2*cos(e)^2-c^6*r^2*d^2*cos(e)^2+4*d^4*c^2*cos(e)^4*r^2+8*d^6*c^3*cos(e)^5*sin(e)+10*r^2*c^4*d^4*cos(e)^2-4*c^4*d^6*cos(e)^4*sin(e)^2-8*c^3*d^2*cos(e)*sin(e)*b^2*r^2-8*c^3*d^4*cos(e)^3*sin(e)*r^4-8*r^6*c^3*d^2*sin(e)*cos(e)-8*c*d^4*cos(e)*sin(e)*b^2*r^2+8*d^5*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)^3+8*c*d^2*cos(e)*sin(e)*b^4*r^2-24*c^3*d^5*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)+2*c*d^2*cos(e)*sin(e)*a^4*r^2+16*c^3*d^5*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*c*d^4*cos(e)*sin(e)*a^2*r^2-8*c^2*d^4*cos(e)^2*sin(e)^2*a^2*r-16*c^3*d^6*cos(e)*sin(e)^5+16*c^2*d^4*cos(e)^2*sin(e)^2*b^2*r+16*c^3*d^6*cos(e)*sin(e)+8*c*d^4*cos(e)^3*sin(e)*r^4+2*c*d^6*cos(e)*sin(e)*r^2+2*c^5*d^2*cos(e)*sin(e)*r^2+16*c^2*d^4*cos(e)^2*sin(e)^2*r^4+4*c^3*d^2*cos(e)*sin(e)*a^2*r^2+4*d^4*r^4*c^2*cos(e)^2-8*c^3*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*r^4-8*c^2*d^6*cos(e)^2*sin(e)^2*r-16*c*d^4*cos(e)*sin(e)*r^4-8*c^4*d^4*cos(e)^2*sin(e)^2*r-4*b^4*r^2*c^2*d^2*cos(e)^2-8*d^3*b^2*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+8*c*d^3*cos(e)^3*(r^2-d^2*cos(e)^2)^(1/2)*r^4+4*d^5*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*d^3*c^3*r*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)-8*c*d^3*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)*r^4+16*c^2*d^3*cos(e)^2*(r^2-d^2*cos(e)^2)^(1/2)*r^4*sin(e)-8*d^5*c^2*(r^2-d^2*cos(e)^2)^(1/2)*sin(e)-8*c^3*r^2*d^4*cos(e)*sin(e)^3+12*c^3*r^2*d^4*cos(e)*sin(e)-4*d^4*r^4*cos(e)^2+4*d^3*a^2*r*c*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+4*a^2*r^2*c^2*b^2*d^2*cos(e)^2+4*a^2*r*c^3*d^4*cos(e)^3*sin(e)-a^4*r^2*c^2*d^2*cos(e)^2+4*c^5*r*d^4*cos(e)^3*sin(e)-8*r^2*c^4*d^4+4*r^6*c^2*d^2+4*c^4*r^2*b^2*d^2*cos(e)^2-2*d^4*a^2*r^2*c^2*cos(e)^2+4*d^4*b^2*r^2*c^2*cos(e)^2-2*c^4*r^2*a^2*d^2*cos(e)^2+4*d^4*r^4+4*c^4*d^4*cos(e)^2*sin(e)^2*r^2+4*c^3*r*d^6*cos(e)^3*sin(e)-8*d^6*c^2-8*b^2*r*c^3*d^4*cos(e)^3*sin(e)-d^6*r^2)^(1/2))/(r^2*c^2-2*c*d^2*cos(e)*sin(e)-2*c*d*cos(e)*(r^2-d^2*cos(e)^2)^(1/2)+d^2)*(r^2-d^2*cos(e)^2)^(1/2))/d/r/(-sin(e)+c*cos(e))
我觉得你的目的不是解这个方程，不管怎样，告你怎么搞吧
