解三角方程时，怎么出错，程序如何修改？？？

笨笨

northwolves 发表于 2023-8-20 13:17

先生你好，能否把秒（6.26429……）的小数点再向后显示20位？？？

Jack315

笨笨 发表于 2023-8-20 13:40
先生你好，能否把秒（6.26429……）的小数点再向后显示20位？？？

先求解方程：

1. Solve[(Sin[x/2])^2 == 1/\[Pi] - 1/4, x]

复制代码

然后四个解中的每一个值可以按下列步骤来得到所需的精度：
1. 计算出角度

1. 角度 = 2 (ArcSin[1/2 Sqrt[(4 - \[Pi])/\[Pi]]])/Degree

复制代码

注：分母上的一个小圆圈代表的就是 Degree 。


2. 计算度——取出整数部分

1. 度 = IntegerPart[N[角度]]

复制代码

注：数值计算命令 N 中现在不需要给精度要求。


3. 计算分——也是取出整数部分

1. 分 = IntegerPart[N[(角度 - 度) 60]]

复制代码

注：“角度-度”后剩下的就是分，
乘 60 取整数部分，同样不需要给精度要求。


4. 计算秒
秒 = N[((角度 - 度) 60 - 分) 60, 20]
注：“((角度 - 度) 60 - 分) 60” 就是秒，
想要多少精度就在 N 命令中指定，
这里给出的是小数后 20 位精度。


最后的结果为：
30 度 18 分 6.264 291 701 702 226 949 60 秒

笨笨

Jack315 发表于 2023-8-20 16:07
先求解方程：

你好，该楼解答不是我想要的，楼上版主northwolves这个代码才是我想要的：

1. DMSList[180/Pi*
   
2.   N[x /. Solve[(Sin[x/2])^2 == 1/\[Pi] - 1/4 && 0 < x < Pi/2,
   
3.       x][[1]]]]

复制代码

请问这个代码怎么修改才能把秒（6.26429……）的小数点再向后显示20位？？？

northwolves

改的话一起都改了：

1. N[DMSList[180/Pi*(x/.Solve[(Sin[x/2])^2==1/\[Pi]-1/4&&0<x<Pi/2,x][[1]])],50]]

复制代码

{30.000000000000000000000000000000000000000000000000,18.000000000000000000000000000000000000000000000000,6.2642917017022269495924530117549292501495143448715}

笨笨

northwolves 发表于 2023-8-20 17:36
改的话一起都改了：

终于搞定了，谢谢版主一路指导：

1. N[DMSList[(x /.
   
2.      Solve[(Sin[x/2])^2 == 1/\[Pi] - 1/4 && 0 < x < \[Pi]/2, x][[1]])/
   
3.    Degree], 50]

复制代码