想求解一个式子很长的方程的数值解，有没有办法让输入简单点

很长的式子写起来很麻烦，也容易出错，很想用一些符合代表式子里面的一部分，比方说用ad代替sin(x^2)，这样式子就看起来简洁一些了，但求解的时候总是报错，不知道该怎么办
u=4;
p=7;
R=5;
u1=pi/3.5;
u2=pi/16;
T=136;
th=acos((u^2+p^2-R^2)/(2*u*p));
K1=tan(th-pi/2+u1);
K2=tan(th-pi/2+u2);
syms th1
om=2*pi*R*cos(th1)/T;
x=R*sin(th1)*cos(om)+p*sin(om);
y=-R*sin(th1)*sin(om)+p*cos(om);
fun1=@(th1)y-r*sin(th)-K1*(x-R*cos(th));
[x,fval,exitflag]=fsolve(fun1,x0,options)
求解未知数是th1的一个方程