
代码如下:
Subscript[h, 0] = 5; theta0 = Pi/4; Subscript[c, 0] = 1;
Subscript[l, 0] = Subscript[h, 0]/Sin[theta0];
Subscript[l, 1] = Sqrt[((Subscript[l, 0]/2)^2) + (Subscript[c, 0]^2) ];
a = Subscript[h, 0] /Tan[theta0];
h = Subscript[h, 0] - x;
theta01 = ArcTan[h/a];
Subscript[l, 01] = h/Sin[theta01];
thetamid = ArcCos[Subscript[l, 01]/(2*Subscript[l, 1] )]


