错误使用 \ 矩阵维度必须一致。 出错 untitled1 (第 25 行) p = A\(-K*diff(sin(x))/mu); % 注意这里的正负号

这个错误是因为矩阵 A 的大小为 N × N，而右侧的向量的大小为 (N-2) × 1，两者的维度不一致。解决这个问题的方法是将右侧的向量大小改为 N × 1，其中第一个和最后一个元素分别为 0 和 -50.6366，即：

b = zeros(N,1);
b(2:N-1) = -K*diff(sin(x))/mu;
b(1) = 0;
b(N) = -50.6366;
p = A\b;

这样就可以使用 \ 运算符求解线性方程组了。