clear,clc syms a b L syms x positive a=120 L=125 b = x/a eqn = a*(1+(8*b^2)/3-(32*b^4)/5+(256*b^6)/7)-L == 0 solve(eqn,x,'MaxDegree',6)