用matlab编写关于cosa=（l4^2+l3^2+l1^2-l2^2)/(2l3*l4-2*l3)的程序，其中l4=57.2,，假设l1=常数，求l2和l3

根据公式：

cosa = (l4^2 + l3^2 + l1^2 - l2^2) / (2 * l3 * l4 - 2 * l3)

可以整理为：

cosa * (2 * l3 * l4 - 2 * l3) = l4^2 + l3^2 + l1^2 - l2^2

2 * l3 * l4 * cosa - 2 * l3 * cosa = l4^2 + l3^2 + l1^2 - l2^2

2 * l3 * (l4 * cosa - cosa) = l4^2 + l3^2 + l1^2 - l2^2

l3 * (2 * l4 * cosa - 2 * cosa) = l4^2 + l3^2 + l1^2 - l2^2

l3 * (l4 * cosa - cosa) = l4^2 / 2 + l3^2 / 2 + l1^2 / 2 - l2^2 / 2

l3 * cosa * (l4 - 1) = l4^2 / 2 + l3^2 / 2 + l1^2 / 2 - l2^2 / 2

l3 = (l4^2 / 2 + l3^2 / 2 + l1^2 / 2 - l2^2 / 2) / (cosa * (l4 - 1))

现在，我们已经得到了l3的表达式，可以使用该表达式和已知的l4的值（57.2）来解出l3的值。然后，我们可以使用已知的l1值（常数）和l3的值来解出l2的值。

下面是Matlab代码实现：

l4 = 57.2; % 已知l4的值
l1 = 10; % 已知l1的值
cosa = 0.8; % 给定的cosa值

l3 = (l4^2 / 2 + l1^2 / 2 - l2^2 / 2) / (cosa * (l4 - 1)); % 求解l3的值

l2 = sqrt(l4^2 + l3^2 + l1^2 - 2 * l3 * l4 * cosa); % 求解l2的值

fprintf('l2 = %f, l3 = %f\n', l2, l3); % 输出结果

请注意，这里假设了l1是常数，因此将其设置为了10进行演示。实际使用时，应该根据具体情况来确定l1的值。