为什么这个方程解不出来？

clc;
clear all;
close all;

syms ra R h delta0 mu k rb
d0=(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')