clc;clear all;close all;
syms ra R h delta0 mu k rbd0=(rb-ra)*tan(delta0);OC3=sqrt(R^2+(d0+h)^2);OC2=OC3-k*d0;R1=sqrt(OC2^2-(d0+h)^2);equ1=(mu*R1-rb)*((mu*R1)^2+d0^2)^(3/2)+ra*rb*mu*R1*d0;solve('equ1=0','rb')
各位大哥大姐帮帮喽!1