高校数学の質問スレ Part438

[過去ﾛｸﾞ] 高校数学の質問スレ Part438 (1002ﾚｽ)
上下前次 1-新
通常表示 512ﾊﾞｲﾄ分割ﾚｽ栞

このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索歴削→次ｽﾚ栞削→次ｽﾚ過去ﾛｸﾞﾒﾆｭｰ

600: １３２人目の素数さん [sage] 2024/11/13(水) 07:40:38.84 ID:mLwXKDzl

>>597
Rに移植してみた。
options(digits=22)
re=matrix(nrow=20,ncol=3)
for(n in 1:20){
for(a in 1:3){
re[n,a]=sqrt(2^(4 %% a)*floor((sqrt(2)+1)^(2*n+1)/4)^2 + 2*(1 %% a) * floor((sqrt(2)+1)^(2*n+1)/4) + (1 %% a))
}
}
re

> re
[,1]             [,2]                [,3]
[1,]                3                4                5.00
[2,]               20               21               29.00
[3,]              119              120              169.00
[4,]              696              697              985.00
[5,]             4059             4060             5741.00
[6,]            23660            23661            33461.00
[7,]           137903           137904           195025.00
[8,]           803760           803761          1136689.00
[9,]          4684659          4684660          6625109.00
[10,]         27304196         27304197         38613965.00
[11,]        159140519        159140520        225058681.00
[12,]        927538920        927538921       1311738121.00
[13,]       5406093003       5406093004       7645370045.00
[14,]      31509019100      31509019101      44560482149.00
[15,]     183648021599     183648021600     259717522849.00
[16,]    1070379110496    1070379110497    1513744654945.00
[17,]    6238626641379    6238626641380    8822750406821.00
[18,]   36361380737780   36361380737781   51422757785981.00
[19,]  211929657785303  211929657785304  299713796309065.00
[20,] 1235216565974037 1235216565974038 1746860020068404.75

Rだとsqrtで丸め誤差がでてしまう。
＞wolfram優秀過ぎ
には同意。
整数解や分数解の計算はWolframが便利。

理工系ならWolfram言語、医系ならR言語が使えるのが普通。
マニアならHaskellかも。
昨今ではPython人気だけど、統計扱うならRだな。

http://rio2016.5ch.net/test/read.cgi/math/1723179542/600

601: １３２人目の素数さん [sage] 2024/11/13(水) 08:13:55.94 ID:mLwXKDzl

>>600
おまけ、Table[x=Floor[(Sqrt[2]+1)^(2n+1)/4];y=Ceiling[(Sqrt[2]+1)^(2n+1)/4];{x,y,Sqrt[x^2+y^2]},{n,1,20}]//Simplify
の方をRに移植。

f=\(n) c(x<-floor((sqrt(2)+1)^(2*n+1)/4),y<-ceiling((sqrt(2)+1)^(2*n+1)/4),sqrt(x^2+y^2))
t(sapply(1:20,f))

[,1]             [,2]                [,3]
[1,]                3                4                5.00
[2,]               20               21               29.00
[3,]              119              120              169.00
[4,]              696              697              985.00
[5,]             4059             4060             5741.00
[6,]            23660            23661            33461.00
[7,]           137903           137904           195025.00
[8,]           803760           803761          1136689.00
[9,]          4684659          4684660          6625109.00
[10,]         27304196         27304197         38613965.00
[11,]        159140519        159140520        225058681.00
[12,]        927538920        927538921       1311738121.00
[13,]       5406093003       5406093004       7645370045.00
[14,]      31509019100      31509019101      44560482149.00
[15,]     183648021599     183648021600     259717522849.00
[16,]    1070379110496    1070379110497    1513744654945.00
[17,]    6238626641379    6238626641380    8822750406821.00
[18,]   36361380737780   36361380737781   51422757785981.00
[19,]  211929657785303  211929657785304  299713796309065.00
[20,] 1235216565974037 1235216565974038 1746860020068404.75
>

http://rio2016.5ch.net/test/read.cgi/math/1723179542/601

上下前次 1-新書関写板覧索設栞歴

あと 400 ﾚｽあります
ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ AAｻﾑﾈｲﾙ

ぬこの手ぬこTOP 0.011s