〓Mathematica〓 - 暇つぶし2ch

〓Mathematica〓 ..

511:503
03/07/09 00:02
>>509
寝ようと思ってたが気づいたので。
多分 L がリストを返してるからだ。

512:１３２人目の素数さん
03/07/09 14:24
ndist2=NormalDistribution[0,4];
p1=Table[ Random[ndist2],{100}];

L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-({p1[[i]] - b}^2/(2 *a^2)))],{i,1,100}]

K1[a_,b_]=D[L[a,b],a]
K2[a_,b_]=D[L[a,b],b]

FindRoot[{R1[a,b]==0,R2[a,b]==0},{a,{1,2}},{b,{-2,4}}]

ソースは、これです。
>>511
あーそれは、どうすれば、解消できますかね。リストで返さずに関数として認識
してもらわないとつらいっすね・・・。

513:498
03/07/09 15:12
あっ最後の式は、
FindRoot[{K1[a,b]==0,K2[a,b]==0},{a,{1,2}},{b,{-2,4}}]
です。

514:１３２人目の素数さん
03/07/09 16:05
今課題をやっています。

誰か偏微分係数の出し方を教えてください。

ちなみにMathematika 4.1を使用しています。
ヘルプの使い方がわからなくて・・・

515:514
03/07/09 16:19
解決すますた。
失礼。

516:503
03/07/09 17:37
>>498
ソースはコピペしてそのまま張れ。
タイプミスもあるしな。

πってPiだよな？
<<Statistics`NormalDistribution`
してるんだよな？

原因はp1のまえにある'{'だ。'{'はリスト構築子。
　 L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-({p1[[i]] - b}^2/(2 *a^2)))],{i,1,100}]

↓におきかえろ。
　 L[a_,b_] = Product[(1/(Sqrt[2*π]*a))*Exp[(-((p1[[i]] - b)^2/(2 *a^2)))],{i,1,100}]

517:１３２人目の素数さん
03/07/09 21:13
このスレッドの人から見てＳ言語はどうですか？
社会学系だとmatlabやmathematicaより使われていますが…

518:498
03/07/10 16:44
>>516
ありがとうございました。それで正規分布の推定は、できました。ただ、本当に
求めたいのは、
L[a1_, a2_, b1_, b2_, b3_] = Sum[y1[[i]]*Log[F[a1 - (b1*
x1[[i]] + b2*x2[[i]] + b3* x3[[i]])]] +
y2[[i]]*
Log[F[a2 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])] - F[
a1 - (b1*x1[[i]] + b2*x2[[i]] +
b3*x3[[i]])]] +
y3[[i]]*Log[
1 - F[a2 - (b1*x1[[i]] + b2*x2[[i]] + b3*x3[[i]])]], {i, 1,
100}];
という結構複雑なもので、これも｛を（にしてみたんですが、やはりリストで返してる
みたいです。Sumの扱い方が悪いんでしょうか・・・・
てか、見にくくてすみません。

519:１３２人目の素数さん
03/07/10 19:28
>>518
a1,a2,b1,b2,b3,F,y1,... のどれかにリストが入ってんだろ…
同じことを何度も書くが、部分だけ張るなよ…

プログラミング自体初めてなのかも知れないが、レベル低すぎる。
一冊Mathematicaの本読め。

520:498
03/07/10 21:10
すみません、一遍に貼れないんで２回に分けますけど、これが全部です。。。

<<Statistics`ContinuousDistributions`

ndist=NormalDistribution[0,1];

x1=Table[ Random[ndist],{100}];

x2=Table[ Random[ndist],{100}];

x3=Table[ Random[ndist],{100}];

y=Table[2*x1[[i]]+3*x2[[i]]-x3[[i]],{i,100}];

y1=Table[If[y[[i]]<=-5,1,0],{i,100}];
y2=Table[If[-5<y[[i]]<=5,1,0],{i,100}];

y3=Table[If[5<y[[i]],1,0],{i,100}];

521:498
03/07/10 21:11
F[x_] = Integrate[Exp[z]/{1 + Exp[z]}^2, {z, -Infinity, x}];

L[a1_,a2_,b1_,b2_,b3_]=Sum[y1[[i]]*Log[F[a1-(
b1*x1[[i]]+b2*x2[[i]]+b3*
x3[[i]])]]+
y2[[i]]*Log[F[a2-(b1*x1[[i]]+b2*x2[[i]]+b3*x3[[i]])]-
F[a1-(b1*x1[[i]]+b2*x2[[i]]+
b3*x3[[i]])]]+y3[[i]]*Log[1-
F[a2-(b1*x1[[i]]+b2*x2[[i]]+b3*x3[[i]])]],{i,1,100}];

K1[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],a1];

K2[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],a2];

K3[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b1];

K4[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b2];

K5[a1_,a2_,b1_,b2_,b3_]=D[L[a1,a2,b1,b2,b3],b3];

FindRoot[
{K1[a1,a2,b1,b2,b3]\[Equal]0,
K2[a1,a2,b1,b2,b3]\[Equal]0,
K3[a1,a2,b1,b2,b3]\[Equal]0,
K4[a1,a2,b1,b2,b3]\[Equal]0,
K5[a1,a2,b1,b2,b3]\[Equal]0},
{a1,0},{a2,1},{b1,1},{b2,3},{b3,1}
]

522:503
03/07/10 21:20
>>521
F[x_]の定義を良く読め。

それからMathematicaでプログラミングするときは
出力を確認しながらやれ。
F[x_] の行をセミコロン無しで打つか、
F[x]とか試しに打ってみながらやってれば間違えないだろ？

あとは自分でやれ。

523:１３２人目の素数さん
03/07/11 14:51
Solve[]や、NSolve[]で、求めた複数の解のうち、
ひとつだけを使いたいのですが、
たとえばt = %[[-1]], t = Last[%]などとしても、
{x ->0.234hogehoge}と　t　に入ってしまいます。

この x->を取り除くにはどうすればいいのでしょうか？

524:523
03/07/11 15:00
x /. %
でできました。トンクス

525:１３２人目の素数さん
03/07/12 00:43
>521
えー、 Sum[<<>>,{i,1,100}] ってあるけど
数値やるなら NSum じゃ駄目なん？

どうせFindRoot なんだし

526:１３２人目の素数さん
03/07/14 23:27
あのー四捨五入するような関数ってないですか？切り捨てでもいいんですけど。
マニュアルで調べてみたんですが見つからなくて・・・・

527:１３２人目の素数さん
03/07/14 23:28
Roundじゃないの？

528:１３２人目の素数さん
03/07/14 23:36
0.5足して切り捨てろ。

529:山崎渉
03/07/15 12:38

　__∧＿∧_
　|（　　＾＾）|　＜寝るぽ（＾＾）
　|＼⌒⌒⌒＼
　＼ |⌒⌒⌒~|　　　　　　　　　山崎渉
　　 ~￣￣￣￣

530:１３２人目の素数さん
03/07/16 23:06
URLﾘﾝｸ(shop.ponycanyon.co.jp)

531:１３２人目の素数さん
03/07/17 07:26
age

532:１３２人目の素数さん
03/07/17 14:11
>>530
これなら3,800円か。安いな。

533:１３２人目の素数さん
03/07/17 16:26
質問なんですが、例えば
x=k/10
y=k+1
For[k=1,k<11,k++,
　　r=x+y]
みたいに繰り返し計算をさせる場合に、このrは、毎回更新されてしまうんで
その結果を格納してリストの形にしたいとするとどーすればいいですか?
勿論、これぐらい単純だとTable使えばいいんですが、やりたいのはもっと
ややこしいんで、結果をリストにする格納する方法を教えてください。
よろしくお願いします。

534:１３２人目の素数さん
03/07/17 19:13
>>533
単純であろうがなかろうがTableでいいはずだ。

どうしてもというなら
l ={}
For[k=1,k<11,k++,r=x+y;l = Prepend[r,l]]
等。

しかしmathematicaではFor文を使う必要はない(その上
遅いので使うべきでもない)ので再考することを勧める。

535:１３２人目の素数さん
03/07/18 14:39
とぇ

536:１３２人目の素数さん
03/07/22 12:25
FindMinimumを使って最小値とそのときの変数の値を求めると、
{234,{a -> 3,b -> 2}}
みたいな感じだと思うんですが、これをいっぱい出して、a,bの平均値を出すには
どうすればいいでしょうか？

537:１３２人目の素数さん
03/07/23 01:12
>>546
漠然としすぎ。

こういうことか？

<< Statistics`DescriptiveStatistics`

data = {x, y} /. (#[[2]] &) /@
　　　　 Flatten[Table[
　　　　　　 FindMinimum[Sin[x + y], {x, i}, {y, j}], {i, 0, 3}, {j, 0, 3}],
　　　　　 1];

Mean /@ Transpose @ data

538:１３２人目の素数さん
03/07/23 01:12
>>537
s/546/537

539:537=538
03/07/23 01:13
536駄・・・まちがえすぎだ。

540:536
03/07/24 13:07
すみません、537のプログラム見て解読してたんですけど、ちょっと良くわかんなくて。
たとえば、
FindMinimum[f[a,b],{a,1},{b,2}]
とした場合、{最小値,{a -> 0.89898},{b->0.1}]
みたいにした場合、a->の部分を消すのはどうすればいいですか？

541:536
03/07/24 14:36
なんか読みにくいですね・・えーと置換なんかでうまくa->の部分を消せない
かなあと・・・無理なんですかねえ。。。

542:１３２人目の素数さん
03/07/24 14:37
開店セール見に来てーーーーーーーーーー！DVDを見ようー！GO！
▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼
開店セール1枚　500円より→→→URLﾘﾝｸ(www.dvd-exp.com)<)　
アダルトDVDの価格はなーんと1枚900円→→→URLﾘﾝｸ(www.dvd-exp.com)<)　
セルビデオ店の人気商品が24時間お買い物が出来ます。
見やすい画面・注文簡単・時間をかけてゆっくり選んでください！
　　　　　　　　　　URLﾘﾝｸ(www.dvd-exp.com)←へGO！
　　　　　　　　　　　DVDエクスプレス

543:１３２人目の素数さん
03/07/24 15:18
>>541
簡単に消せるように"->"つきで出力されている。

s = {a -> 1, b -> 2}
{a,b} /. s

544:１３２人目の素数さん
03/07/24 22:08
s はリストみたいになっているから
s[[1]] や s[[2]] とかで取り出しても良いんじゃない？

545:１３２人目の素数さん
03/07/25 04:10
>>544
543で解決してます。”a->”みたいなのを取り除きたいってのが541の質問だろ。
ってか、房は黙ってろ（ｗ

546:536
03/07/25 12:02
>>543
ありがとうございました。出来ました。
>>544
それだと、a -> 0.89898などと出ます。

547:１３２人目の素数さん
03/07/25 13:33

s[[2,2]]
とやるんじゃ

548:１３２人目の素数さん
03/07/25 14:11
【クレジットカードのショッピング枠を現金化！】

クレジットカードで買い物をすると、商品代金の８０～９０％のキャッシュバック！！
１０万円の買い物をすれば８～９万円の現金がお手元に！！

■■朝日カード■■
URLﾘﾝｸ(www2.pekori.to)

549:１３２人目の素数さん
03/07/25 21:29
>>547
自分で一回やってみれば？違うのよくわかるから。
ほとんど使ったことないでしょ？分りもしないのにアドバイスしようとしない
方がいいじゃないのか？

550:１３２人目の素数さん
03/07/25 22:37
>549

残念ですが{a->1,b->2}[[2,2]] で 2 が出てくるよ。
[[]] は別に行列以外でも使える、数式構造の場所指定でした。

まあこのやり方は、-> の値を取り出すのにふさわしくないとは
僕も思いますがね。

551:１３２人目の素数さん
03/07/25 22:40
In[4]:=
{a -> 1, b -> 2} [[{1,2},2]]
Out[4]=
{1,2}

List[Rule[a, 1], Rule[b, 2]]

In[6]:=
FullForm[{{1,2},{a,b}}]
Out[6]//FullForm=
List[List[1,2],List[a,b]]

552:１３２人目の素数さん
03/07/25 22:52
>>550
質問の意図と違うから。その意味での違うなんだけど？なにがでしただよｗ

553:550
03/07/25 23:02
>552

確かに s[[2,2]] だと、本当に欲しい {1,2} は
出てきませんね。失礼スマソ。

ただ -> に対して、通常考える /. 以外の [[]] 指定方法に
目を付けた点は悪くはないのでは、とちょっと擁護
してみたかったのでした。

554:１３２人目の素数さん
03/07/26 10:30
>>553
> 　目を付けた点は悪くはないのでは、とちょっと擁護

何がやりたいのかわかりにくい、バグに弱い、
一般的な状況(a^2 + aを求めるとか)では使えない…
センス悪すぎるよ。

555:１３２人目の素数さん
03/07/26 20:47
リストから第一行目、第十一行目、第二十一行目・・・というふうに取り出したい
んですけど、どーすればいいですか？

556:１３２人目の素数さん
03/07/26 22:46
>554

Mathematicaのルールどおりだからバグに弱くない
おまえがセンス悪、
それとも正確ﾜﾙか

557:１３２人目の素数さん
03/07/26 23:07
　　　　　　　　　　　　　　　 _,.. ----　.._
　　　　　　　　　　　　　　,. '"　　　　　　　｀丶、
　　　　　　　　　　　／　　　　　　　　　　　　` ､
　　　　　　　 ,..-‐/　　　 ...:　 ,ｨ　,.ｉ .∧　,　　ヽ.
.　　　　　　　,:'　　.l　.::;',. :::;/..://:: /,':/　 ', l､ .i　　ヽ
. 　　　　　 ,'　　..::| .::;',' :;:','フ'７フ''7/　　 ',.ﾄ',_|,　, ',.',
　　　　　　　,' 　 .::::::!'''l/!:;'／　/'ﾞ　 /　　　'! ﾞ;:|:､.|､| 'l
. 　　　　 ,'. 　.:::::::{　ｌ'.l/　　､_　　_,.　　　　　　'l/',|.';|
　　　　　　ｌ　 :::::::::::';、ヾ　　　￣　　　 `‐-‐'/! ';. '
. 　　　　 !　:::::::::::/ `‐､　　　　　　　　ゝ　　　|'ﾞ　|
　　　　　　 | ::::::::／　　　＼　　　　､＿,　_.,.,_　ﾉ:::　!　
　　　　　　 |::::／. 　　　　_rl｀': ､_　　　///;ﾄ,ﾞ;:::::./ 　　　マターリと良スレ。
..　　　　　　`´　　　　　 /＼＼　 `i;┬:////ﾞlﾞl ヾ/　　　　　　　　　　数ヲタがうなぎを食べられますように・・・
　　　　　　　　 ,.:く::::::::`:､＼　〉lﾞ:l 　/　!.|
.　　　　　　　　／:.:.:.:＼:.:.:.:.`:､ｿ/:.:| 　　| |
　　　　　　　　 /.:.:.:.:.:.:.:.:.:＼:.:.:.:У:.:;l　　 /./
.　　　　　　　 /:.:.:.:.:.:.:.ｒ'´`‐,`､:/.,.:‐{　　 |　!`:､
　　　　　　 ,'.:.:.:.:.:.:.:.:.';_,ﾟ.,ﾉ.:./,:':.:.:.:',　　|　|`､:|
　　　　　　 !:.:.:.:.:.:.:.:.:.:.ﾞ､:.::／:.:.:.:.:.:.ヽ, /　,!:.:`､

558:１３２人目の素数さん
03/07/27 00:52
鰻食いてーー

559:１３２人目の素数さん
03/07/27 00:53
でも貧乏だからとても買えねーよ(;´Д`)

560:１３２人目の素数さん
03/07/27 01:07
>>556
> 　Mathematicaのルールどおりだからバグに弱くない
あきれますな。悪いこと言わんからもっと謙虚になれ。
プログラムくんだ経験自体無いだろ？

たとえばFindRootで解が見付からなくてそのまま式自体が
かえったらどうなる？

[[]]ってのはかなりError proneだぞ。

561:１３２人目の素数さん
03/07/27 09:37
>555
tensor_a[..., Table[1+10 i,{i,0,5}], ....]
てな感じ。ポイントは添字指定で数字ではなく
リストを使うこと。

562:550
03/07/27 09:53
>560
/. でも FindRoot そのままの式が返ってくると
あまりうれしくない気がするけど。。。

/. を使うと a^2+a /. みたいな純関数的使い方が出来て
見通しが良くなるのは確かに同意。
つまり /. と -> はペアで使うのが基本なのは当然だけど
Lisp 的な感覚として式の木構造を考えることが出来る [[]]も
無下に NG とするのもどうかも。。。。

あきらめ悪いっすかね

563:１３２人目の素数さん
03/07/27 10:59
Mathmatica 4を使っています。
Module[{i, j},
　ans = {};
　Do[
　　Do[
　　　ans = Append[ans, {i, j}], {j, i + 1, 3}
　　], {i, 2}
　]
　ans
]
を実行すると
{{Null, 2 Null}, {Null, 3 Null}, {2 Null, 3 Null}}
と変なNullという文字が引っ付いてくるのですが、
どうすれば、まともな答えが戻ってくるのでしょうか？

564:１３２人目の素数さん
03/07/27 13:15
>>562
エラーメッセージが出て、文字がそのままかえってくる
のと、
「何も言わずに、間違った答えがかえってくる」
のでは大違いだぞ。

全ては式である、ということの一例として[[]]がつかえることを
理解していることは必要だが、無闇に使うのは否定すべきだよ。

565:１３２人目の素数さん
03/07/27 13:19
>>563
　　]
　　ans
->
　　];
　　ans

566:１３２人目の素数さん
03/07/27 14:01
>>563
おれは>>534だから一応Tableによる解決を書いておく。
このほうがすくなくとも数倍速い。

UpperInds[n_Integer] := Flatten[Table[Table[{i,j},{j, i + 1, n}],{i,n - 1}],1]

567:563
03/07/27 14:45
>>565
ありがとうございます。解決ししました。

>>566
ありがとうございます。
じつは、「内側のループのなかで、i,jを使った条件判定をやった結果、条件
を満たした時のみ、i,jに依存する要素をansに加えていく」というのがやりたい
ことなので、ちょっとそのままでは教えていただいたやりかたは使えないよう
です。

568:１３２人目の素数さん
03/07/27 15:15
>>567
それなら
UpperInd[n_Integer] := Module[
　　 {ConvertFunc = Sin[#[[1]] + #[[2]]] &, CondFunc = EvenQ[#[[2]]] &},
　　 ConvertFunc /@ Flatten[Table[
　　　　　 Select[Table[{i, j}, {j, i + 1, n}], CondFunc], {i, n - 1}], 1]]
あるいは、
Table[If[..., ..., Gomi],...]としてあとでGomiを除くとか。

綺麗なのは、

ConvertFunc /@ Select[Flatten[Table[{i,j},{i,n},{j,n}],1], CondFunc]

(CondFuncのなかで#[[1]] < #[[2]]も判定)
みたいなやつかな。

569:１３２人目の素数さん
03/07/29 01:55
すみません、ものすごく初心者なんですけど教えてください。
非線形の最適化ってできるのですか？
お願いします。

570:１３２人目の素数さん
03/07/29 11:57
ウルフラムさん
セルオ－トマトンの原理で宇宙の起源を解明しようとしてるらしい
一体どうしちゃったんだろうか

571:１３２人目の素数さん
03/07/29 11:58
なんかMathematica5.0ての海外では売ってるらしいんですけど。

572:１３２人目の素数さん
03/07/29 22:13
a[0]=0
a[1]=1
a[n+2]=a[a+1]+a[n]
みたいな漸化式を使ってa[100]とかを求めるにはどうすればいいですか？

573:１３２人目の素数さん
03/07/29 22:39
age

574:１３２人目の素数さん
03/07/29 22:58
>>572
a[0] = 0;
a[1] = 1;
a[n_] := a[n - 1] + a[n - 2]

575:574
03/07/29 23:00
>>572
100は計算回数(2^100)が多すぎる求められんと思うぞ。
計算結果を覚えながらやれば別だが:
b[0] = 0;
b[1] = 1;
b[n_] := (b[n] = b[n - 1] + b[n - 2])
みたいに。

576:１３２人目の素数さん
03/07/29 23:04
レスありがとうございます！
早速やってみます

577:１３２人目の素数さん
03/07/29 23:06
出来ました！すごく感謝！

578:１３２人目の素数さん
03/07/30 02:25
x+y=1
xy=-1
というとても簡単な連立方程式を解かせたいのですが、
Solve[{x + y == 1, x - y == -1}, {x, y}]
とやっても上手くいきません。どうすればよいでしょうか？

579:１３２人目の素数さん
03/07/30 02:26
>>578
ミスったｗ
Solve[{x + y == 1, x y == -1}, {x, y}]
これでやって出来なかったです。お願いします。

580:１３２人目の素数さん
03/07/30 10:56
>>579
できるけど。

581:１３２人目の素数さん
03/07/30 11:47
買ったけど使い方がさっぱりわからない。友達に売ろ。

582:１３２人目の素数さん
03/07/30 11:53
ムカついたんでnyに流そっかな・・・

583:１３２人目の素数さん
03/07/30 14:31
>>575
b[1244]まではそれで出るんですけど、
b[1245]以降はHoldなんたらとか出て計算できないんですが
どうすればもっと大きい数まで表示出来るようになるでしょうか？

584:math.1st ◆M9pCfogc9g
03/07/30 14:51
Re:>583 再起呼び出しの関数は、ネストできる数に制限があることから、使わない方がよい。
使用メモリを増やして計算すると大きい数も表示できて、速い。

585:１３２人目の素数さん
03/07/30 15:01
>>584
ありがとうございます。
もうひとつ質問させてください。
n≦10の時
a[n_]:=a[n-1]+2
n＞10の時
a[n_]:=a[n-1]+1
のように、範囲を指定したい時はどうすれば良いでしょうか？

586:１３２人目の素数さん
03/07/30 21:34
>>585
そんな問題Cでやった方が速い。2^32進数の作り方もついでに勉強。
まったく典型的「ユーザさん」だね。
あなた甲南大学とやらの学生ですか？
そういう丸投げ的利用法は、あとあとためにならないですよ。

587:585
03/07/30 23:42
>>586
いえ、高校生です。

588:１３２人目の素数さん
03/07/30 23:55
高校生ならいいってもんでもないな。

589:１３２人目の素数さん
03/07/31 01:49
>>583
末尾再帰にすれば良い(演習問題)。
関数型言語では、関数は全て再帰で書く。

末尾再帰に書き換えないのであれば、
$RecursionLimit=3000
などと再帰の上限階数を増やす。

末尾再帰については自分で調べろ。

>>584
くだらんレスをするな。

>>585
例えば:
a[n_] /; n <= 10 := a[n - 1] + 2
a[n_] /; n > 10 := a[n - 1] + 1

>>586
丸なげ反対はわかるが、
> のように、範囲を指定したい時はどうすれば良いでしょうか？
これは言語仕様について聞いているんだろう。

590:１３２人目の素数さん
03/07/31 23:36
>>585
その後 FullForm で、出てくる式の構造を暇な時に見ると
将来、ひょっとしたら役に立つ可能性があるかもしれない。

591:585
03/08/01 02:39
>>589
ありがとうございます。
末尾再帰ってのは調べてみたんですけどよく分かりませんでした…。
n = 2000; a = 0; b = 1; While[n >= 1,
If[n ≠ 1, a = a + b; n--, Print[b]; Break[]];
If[n ≠ 1, b = a + b; n--, Print[a]; Break[]]]; If[n == 0, Print[0]]
結局上のを書いてみたんですが改良点が有ったらご指摘お願いします。

f[n_]:=a = 0; b = 1; While[n >= 1,
If[n ≠ 1, a = a + b; n--, Print[b]; Break[]];
If[n ≠ 1, b = a + b; n--, Print[a]; Break[]]]; If[n == 0, Print[0]]
あと、上のようにやってf[2000]とかやっても0とか出ちゃうんですけど、
何が間違ってるのでしょうか…。

592:１３２人目の素数さん
03/08/01 11:58
高校生でこんなことをやってるのか。
凄いな

593:
03/08/01 12:39
え～今までMathematicaなんて使ったことないのに、卒研で使えと言われまして・・・
さっぱりわからないんで教えていただけないでしょうか？

ワイブル分布において
平均Mean[dist] == 1と標準偏差StandardDeviation[dist] ==1/3の連立方程式を
αとβについてSolveで解こうとしたんですが、
『逆関数が使用されています．値は多価逆関数のため失われている可能性があります．』
と出て解けないんですよ。

これってどうすればいいのでしょうか？

594:１３２人目の素数さん
03/08/01 16:39
Play[Sin[700 t + 25 t Sin[350 t]], {t, 0, 4}]

595:１３２人目の素数さん
03/08/01 16:53
>>593
解けないんじゃなくて、そのメッセージの後に代表値は出てるよな？（ｗ

596:１３２人目の素数さん
03/08/01 17:07
Mathematicaとは関係ない話なんだけどwolfram reserchの日本語の
トップページ(URLﾘﾝｸ(www.wolfram.co.jp))ってさ、リロードするたびに違う
種類の画像とキャッチコピーが表示されるじゃん。キャッチコピーは
「世界の計算方法」とか「数学、科学、工学、教育における革新的言語」
とか色々あるんだけど。あれの中の一つに「技術計算の未来を定義」
ってのがあったのね。それを読んだときに「未来を予測する最良の方法
は、それを発明してしまうことだ」っていうAlan C. Kayの言葉を思い出し
て、「ああそうか、数学で何か発明するってはまさしく何か定義すること
だな」と思って少し感慨を覚えたんだよ。どうでもいい話だけど。

ってか英語の方のトップページ見たらMathematica5って書いてあるじゃん！
6/23に出たのか。全く知らなかった。

597:１３２人目の素数さん
03/08/01 17:31
例えば、
p1=Plot[Sin[x],{x,-10,10}];
というように最後にセミコロン「；」を付けてもグラフが表示されてしまうのですが
これを表示させないようにするにはどうすればよいのですか？

598:ビッグバン宇宙論は完全に間違いでした！！！！！
03/08/01 18:39
科学者よ、恥を知れ！
ビッグバン宇宙論は完全に間違いだった！
科学の原則を無視した、デタラメのインチキ理論だったのだ。
そして、そのビッグバン宇宙論の世界的な浸透は
アメリカ、ユダヤ・キリスト教勢力による世界支配のための思想的な戦略なのだ！
また、ビッグバン宇宙論の思想によって戦争が起こり、
貧富の差がひらき、終末的な絶望感が世界に蔓延しているのだ。
ビッグバン宇宙論は世界の平和を揺るがす、悪の元凶となっているのだ。
ビッグバン宇宙論とは、
「宇宙は『無』からビッグバン（大爆発）によって誕生した」という理論である。
この理論は、ユダヤ・キリスト教の創造神話（神が天地を創造した）そのものである。
ビッグバン宇宙論の実態は、科学理論ではなく宗教思想なのである。
『無』は科学的に証明できるものではなく、
そして、『無からの誕生』も科学では証明できるものではないのだ。
ビッグバン宇宙論が科学の正統であるという思想を、世界中の人々に
浸透させる戦略が成功したことにより、ユダヤ・キリスト教勢力の
世界における優位性が確立されていったのだ。（２０世紀に）
そして、その思想的支配の最大の例が、アメリカやイギリスによる
イラク戦争なのだ。
ビッグバン宇宙論の浸透により、世界中に終末思想（世界の終わり）が蔓延してしまっている。
そのことにより、自己中心的、せつな的、短絡的な考え方が社会に広がっている。
科学的に間違っているビッグバン宇宙論から脱却しなければならない。
そして、宇宙は無限だということを理解しなければならない。
人間は本当の宇宙観、世界観を構築し、
新しい時代に進んでいかなければならないのだ。
ビッグバン宇宙論が世界を支配している限り、平和な世界にはならないのだ。
そのことを科学者は重く受けとめるべきである。
平和な世界を！！！！！！！！！！！！！！！！！！！！！！！

599:１３２人目の素数さん
03/08/01 18:49
　　　　　　　　　　　,-rヽ、　　　　　　 __,.- 、
　　　　　　　　　　 /　）ｧﾍ　　　　　 / 　　 } 　
　　　　　　　／ ,ｨ　,.ﾆｨ　　　　 /´ヽo‐´
　　　　　　　／＿〈ﾚ'´∠-‐‐- ､__　/　　ｌ
　　　　ﾄ<,__,＞　 _ ,,,,, _　　　＜/　　　|
　　　　（＼_ ､__∠-､‐､‐ヽ　　､ .l　ｌ　 /
　　　　　〉／ _.ィ‐'´ヽﾄ､Alヽﾄ､ヽｌヽl　　　|
　　　イァﾉ.ｲ'"_　　　 ;==､ヽ　ｌ　　 L._
　　　　,ｲ |　　ｌ /lﾊ　　　}Viゝ } /|　　　ｌ:.:.ヽ
　　　／.:L.ヽﾊﾄﾞ上!　　　ﾞ‐'　ﾉﾄ､　　レ' {
　　〈-‐'´.-'´ >ゝ" ` _　....、´ ´7/､｀ー'´　.ｲ
　　　ヽ‐'´,-／ィヘ、　ヽ _ﾉ　　,彳　｀ｰ-‐'´ /
　　　　｀´l　　ﾉィﾍｨ＞､__,.. ‐'´ {　 l.　　　/
　　　　 l　 ! /〈__.L..､ﾉ __.-‐'´| 　! 　　〈
. 　　　　l　! ＞'´-､イ⌒'ｰ-‐'´　　、　　/
　　　　　　 | /＿,ィ´i　l. . : : : : . . ,' 　｀ｰ./
　　　　　〉､: : :.l _,イ: : : : : : : : !　　 /　
　　　　　　/. :! : 　´　　 : : : : : :/　　 / 　
　　　　　l: : : 　　　　　　　　　 /　
　　　　　 /. :　_. -‐- ._　　　. : : : :.〈
　　　　／_,.ィ´　　　　　｀ｰ-､;_: : : :|
　　　＼::ﾉ　　 j　__.　-‐ ‐- エ＞
　　　　　 /‐ '' "__,... -‐‐ｒｰ､-､_|　　
　　　　／-‐ｧ' ´/　 |　　l　ヽ　 ` 、
　　 .／　／　 /　　l 　　! 　`、　　ヽ、

600:１３２人目の素数さん
03/08/01 19:29
おい、そのAAはロリスレに貼れ

601:ぼるじょあ ◆yBEncckFOU
03/08/02 03:01
　　　　∧＿∧　 ∧＿∧
ﾋﾟｭ.ｰ　（　・３・） (　　＾＾）＜これからも僕たちを応援して下さいね（＾＾）。
　　＝〔~∪￣￣￣∪￣￣〕
　　＝ ◎―――◎ 　　　　　　　　　　　　　　　　　　　　　山崎渉&ぼるじょあ

602:１３２人目の素数さん
03/08/02 03:41
>>591
> 　n = 2000; a = 0; b = 1; While[n >= 1,
> 　...
論外なのでコメントできん。一冊以上Mathematicaの本を読め。

末尾再帰を用いた解答
f[n_] := f[n, 0, 1];
f[0, r1_ , r2_ ] := r1;
f[n_, r1_, r2_] := f[n - 1, r2, r1 + r2];

> 　f[n_]:=a = 0; b = 1; While[n >= 1,
a = 0;のセミコロンで文が切れている。
f[n_] : = a
b = 1;
と入力しているのと同じだ。

603:１３２人目の素数さん
03/08/04 00:48
>>602
論外ですかｗ
そうゆう本高そうで…。今度学校の図書館で詠んでみます。
ありがとうございました。

604:１３２人目の素数さん
03/08/04 08:23
597に応えてあげて

605:１３２人目の素数さん
03/08/04 20:18
>>604=>>597

606:１３２人目の素数さん
03/08/04 23:22
>595
日本国内向けに日本語文書作業中。
もっとも買う人がいないから強制バンドル
価格Up 中。

日本は非英語で高度な技術内容が学べ得るけど
その分仲介業者が入る余地が多くなってしまった。。。。

607:１３２人目の素数さん
03/08/05 23:48
Ver.5使ってます。
起動直後にガイドが出るくらいで
見た目はそんなに変わってない。

608:１３２人目の素数さん
03/08/05 23:49
★おまんこナビ登場！！★日本の美しきオマンコ★
URLﾘﾝｸ(endou.kir.jp)

609:１３２人目の素数さん
03/08/05 23:50
あなたの近所のあの娘が出てるかも
いつもすまして歩いてるあの娘が
乱れまくり！！
URLﾘﾝｸ(www.j-girlmovie.com)

610:１３２人目の素数さん
03/08/06 15:32
Sin[a°]=1/2なるaを求めるにはどうすればよいでしょうか？
Solve[Sin[a°]==1/2]では求められませんでした。

611:１３２人目の素数さん
03/08/07 00:11
Solve[Sin[a] == 1/2, a]

612:１３２人目の素数さん
03/08/07 14:32
>>611
違う。

613:１３２人目の素数さん
03/08/07 14:59
>>612
Solve[Sin[a Degree] == 1/2, a]
で
Solve[Sin[a] == 1/2, a]
と同じ程度に、求めることができているので、
質問者が何を求めているのかよくわからない。

Solve[Sin[a Degree] == 1/2, a] // N
ってのが欲しいのか。

614:１３２人目の素数さん
03/08/07 15:25
>>613
スペースが必要だったのか。ありがとう。

615:１３２人目の素数さん
03/08/09 12:52
>>614
そういう誤解か。
a\[Degree]はそれ自体がシンボルだから、
Solve[Sin[a\[Degree]] == 1/2, a\[Degree]]
はできるが、aは求められない。

スペース(かけ算)が見た目に嫌なら、
a\[InvisibleTimes]\[Degree]とすればよし。

616:１３２人目の素数さん
03/08/09 15:56
>>597
だれも書かないな。
$DisplayFunction = Identity

617:１３２人目の素数さん
03/08/12 22:43
>616
正確にいきましょう。

p1 = ListPlot[Table[{i, i*Random[]}, {i, 1, 8}],
PlotStyle -> RGBColor[1, 0, 0], DisplayFunction -> Identity];
p2 = Plot[Exp[x], {x, -2, 2},
PlotStyle -> RGBColor[0, 0, 1], DisplayFunction -> Identity];
Show[{p1, p2}, DisplayFunction -> $DisplayFunction];

618:１３２人目の素数さん
03/08/13 00:54
>>617
それのどこが正確なんだ？
PlotStyleなんて何の関係もないし。
# 意味がわかってないままコピペしているようにしかみえないのだが。

で、これのどこが正確でないんだ？
$DisplayFunction = Identity;
Plot[Sin[x],{x,0,1}];
# オプションの方がいい場合「も」あるがな。

619:１３２人目の素数さん
03/08/13 01:56
>>617-618
要はオプションを関数の中で指定しているかグローバルに指定しているかの違いだろ。

620:１３２人目の素数さん
03/08/13 03:57
>>619
違う。
$DisplayFunction は DisplayFunctionというオプションのデフォルト設定
を与えるものであって、オプションそのものではない。
オプションを設定することと、
(オプションが省略された時の)デフォルトを変更することは、厳密には別のことだ。

オプションを関数の外で設定するというのは
SetOptions[Plot, DisplayFunction -> Identity]
だ。

621:山崎渉
03/08/15 18:32
　　　 (⌒V⌒)
　　　│ ＾＾ │＜これからも僕を応援して下さいね（＾＾）。
　　⊂|　　　　|つ
　　　（＿）（＿）　　　　　　　　　　　　　　　　　　　　　　山崎パン

622:１３２人目の素数さん
03/08/18 18:59
おーい、Ver5持ってる奴
x^100-1==0
を解いて、根の絶対値を調べてくれ。
１にならないから。

623:tebe
03/08/18 19:06
半年間もかけてMathematicaの使い方を学ぶという
非常にエグイ授業を大学でやったよ。
Mathematicaって機能ありすぎて、僕みたいな素人
には宝の持ち腐れって感じだな

624:１３２人目の素数さん
03/08/18 19:11
｢今私のオッパイ見てたでしょう。｣ってそりゃ見ますって。
これだけの立派なオッパイが目の前にあったら。
淫乱女教師がフェラ、パイズリで同僚の先生を責めまくります。
挿入部のアップも豊富なのでオマンコもモロ見えですよ。
無料動画をどうぞ
URLﾘﾝｸ(www.exciteroom.com)

625:１３２人目の素数さん
03/08/18 19:19
Mathematicaの有名なバグ

Integrate[1/x^2, {x, -1, 1}]

古いバージョンだと２と評価してくれるんですが
新しいバージョンだと∞になってしまいます。
Ｖｅｒ５だとどうなんでしょうか？

626:１３２人目の素数さん
03/08/18 21:50
>>625
はぁ？
バグってるのはあなたの頭では？

627:１３２人目の素数さん
03/08/20 01:54
　1　　　　　　　　1
∫(1/x^2)dx＝-1/x｜＝-1/1-(-1/(-1))＝－２
　-1　　　　　　　-1

じゃないの？

628:１３２人目の素数さん
03/08/20 01:56
>>627
ずれた
　1　　　　　　　　　　1
∫(1/x^2)dx＝-1/x｜＝-1/1-(-1/(-1))＝－２
　-1　　　　　　　　　-1

じゃないの？

629:１３２人目の素数さん
03/08/20 01:59
>>628
1/x^2>0だから積分値は正だろ？

630:_
03/08/20 02:03
URLﾘﾝｸ(homepage.mac.com)

631:１３２人目の素数さん
03/08/20 02:06
>>625-629
ｘ＝０で発散しているから∞だろ

632:１３２人目の素数さん
03/08/20 08:37
625がまぬけだったということでよろしいでつか

633:１３２人目の素数さん
03/08/20 10:12
>>632
否、これは確かに古いバージョンでは２ではないが－２と評価しているぞ。
いまのバージョンだと∞に修正されてはいるが。

634:supermathmania ◆ViEu89Okng
03/08/20 10:50
発散積分の有限部分は確かに-2になる。
Integrate関数は超関数積分なども一部サポートしている。
これとは関係ないが、主値積分もオプション指定でできる。

ついでに、Integrate[Exp[1/(1-x^2)],{x,-1,1}]と
NIntegrate[Exp[1/(1-x^2)],{x,-1,1}]の結果を報告していただきたい。
期待しないで待っているぞ。

635:１３２人目の素数さん
03/08/20 18:05
>>625は古いバージョンと新しいバージョンの挙動の違いを言ってるだけだろ
積分値が正しい間違ってるの話をしているわけではないと思われ

636:132人目の素数さん
03/08/20 19:54
>>634
NIntegrate[Exp[1/(1-x^2)],{x,-0.998,0.998}]=3.88057*10^103でつた。
Integrate[Exp[1/(1-x^2)],{x,-1,1}]は、答えが帰ってきません（；＿；）
mathematica 4.2 for Linux-PPC

637:１３２人目の素数さん
03/08/20 20:50
ちなみにMapleの結果

int(exp(1/(1-x^2)),x);

/infinity \
| ----- 2 |
| \ hypergeom([1/2, _k1], [3/2], x )|
x | ) --------------------------------|
| / GAMMA(1 + _k1) |
| ----- |
\_k1 = 0 /

638:１３２人目の素数さん
03/08/21 01:41
オクで割れマセマティカを5000円で売ってるやつがいる。
タイーホされてほしい。

639:１３２人目の素数さん
03/08/21 01:55
>>638
通報すれ

640:１３２人目の素数さん
03/08/21 02:06
>>637
そんなの超幾何関数使って元の積分を表現しただけじゃん。

641:１３２人目の素数さん
03/08/21 03:33
0.11111111111111111111111111
って入力して実行すると
0.11111111111111111111111111
って出て、それをクリックしてエンター押すと
0.1111111111111111111111111099999999999964597
って出るんですがどうゆう意味があるのでしょうか？かなり気になるんですが。

642:１３２人目の素数さん
03/08/21 03:40
意味はない。Mathematicaが糞なだけ。（ｗ

643:１３２人目の素数さん
03/08/21 05:45
>>641
ネタで言ってるのかどうか知らんが・・・

機械精度の演算がどういう仕組みで行われているかを一度考えて見れ。

644:１３２人目の素数さん
03/08/21 07:45
それにしてもクソですな。そのくらい改善できそうなものだが。

645:supermathmania ◆ViEu89Okng
03/08/21 15:58
Integrate[Exp[1/(x^2-1),{x,-1,1}]
NIntegrate[Exp[1/(x^2-1),{x,-1,1}]
と云おうとしたら、間違えてしまった。
これなら計算できるのかな？

646:１３２人目の素数さん
03/08/21 19:01
>>644
あのさ、Mathematica では小数点を付けた数値は機械精度
で計算されるとマニュアルに書いてあるから、
0.11111111111111111111111111
と書いた時点でこれは
11111111111111111111111111/100000000000000000000000000
とは全然違う値になるわけだ。だから、もし
0.11111111111111111111111110000000000000000000
とかになったら、Mathematicaは機械精度の計算を唄っておきながら
実際にはサポートされていないことになる。

この前も同じような質問があったけど、今の大学って数値計算の理論
とか教えないんだろうか。もし習ったことが無いなら、そういう技術者や
科学者にとって必要最小限のことも教えない大学なんて行く意味がないよ。
>>641や>>642、>>644が通っている大学の学長は、こういう不祥事（と
言っても良いよね）に対する記者会見を開いても良いぐらいだと思う。鶴を
燃やした事件みたいにさ。高校生だったら、まあ仕方ないけど。

いや、今日は社会人になっても、２次方程式の解を計算するルーチンで
解の公式をそのまま使うような馬鹿に遭遇したんで、イライラしてるのかも
しれん。桁落ちが起きるだろ、アホが！お前の先輩の文型出身の奴でも
ちゃんと場合わけして回避してたぞ。

647:１３２人目の素数さん
03/08/21 20:13
NIntegrate[Exp[1/(x^2-1)],{x,-1,1}]
=0.443994

ですな。

648:１３２人目の素数さん
03/08/21 20:14
数値解析の講義なんか、あったって選択だろ。
そんなもんやってる時間があったら、
純粋数学やるわいな
数学専攻ならな

649:１３２人目の素数さん
03/08/21 20:35
>>646
俺は高校生です。仕方ないです。

650:１３２人目の素数さん
03/08/21 22:29
>>646
> 　あのさ、Mathematica では小数点を付けた数値は機械精度
> 　で計算されるとマニュアルに書いてあるから、

こんなことは書いてない。
おまえの書いた
　 0.11111111111111111111111111
はArbitrary-Precision Numberだ。
つまり、多倍長数の内部表現が二進であるか十進であるかという実装の問題だ。
偉そうなこと書くんだったらもっと勉強しろ。
一目見たらdoubleで書けないことぐらいわかれよ。

64bit以上のCPUつかってんだったら機械精度かも知れんがどっちにしろ
嘘だしな。

651:650
03/08/21 22:34
一応書いとくが、二進で結構だと考えてるよ。
# 上の馬鹿といっしょにされてはかなわん

652:１３２人目の素数さん
03/08/22 05:18
このスレ、人を罵倒するやつが多いな。

653:644
03/08/22 06:20
>>646
おいおいパターソン&ヘネシーくらい読んだことあるぜ。
0.1111111111111111111111111 != 111111111111111111111/1000000000000000000
(桁は適当ね)になっている実装がクソだと言ってるの。
機械精度なんて本来は一切要らないんだから。
そんなものこの時代のマシン特有の都合だろ？
もちろんMathematicaブックを読めば$MachinePrecisionや
N[]や12.3`なんかについてあれこれ載っている(はず)。
しかし、精度という基本的なことがらに付いてあれこれ頭を悩ませる実装および
言語体系が「クソですな。改善できそうなもんだが」というわけ。

掲示板の書き込みで憂さ晴らすのは不健康だね。

654:１３２人目の素数さん
03/08/22 06:49
よくわからんがマスマティカはくそってことでつか？

655:１３２人目の素数さん
03/08/22 06:50
>>653
６４６は書きすぎだと思うが、

> 機械精度なんて本来は一切要らないんだから。
> そんなものこの時代のマシン特有の都合だろ？

それは無いだろー

656:１３２人目の素数さん
03/08/22 11:57
>>646
知りもしないものを「クソ」あつかいか…
それに少し上のレスぐらい読め。
機械精度なんて関係ないんだよ。

計算機の上では有限の精度しか実現できないことをしっているよな？
0.1111111111111111111111111 == 111111111111111111111/1000000000000000000
にしたところで何か計算したとたんにほぼ確実に10進数表示での
数学的な等号は成立しなくなるんだぞ。
何か意味があるか？
精度の保証だけしとけば十分だろ。

Precision, Accuracyについてマニュアルを読め。

657:１３２人目の素数さん
03/08/22 12:09
>>646
この時代のマシンのことを無視してプログラム書いてどうすんだよ。

おまえが書くプログラムってどれも実用には耐えないんだろうな。

658:１３２人目の素数さん
03/08/22 20:32
私は高校生なので優しくしてください。
0.1111111111111111111111111を任意精度として扱ってもらう簡単な方法はないですか？
いちいち111111111111111111111/1000000000000000000のように入力しなおさいと正しい結果が得られないのは不便なんですけど。

659:644
03/08/22 21:53
さあ香ばしくなってまいりました。これも5発売フィーバーか？

>>658 に優しいアシストをすると、そういう小数でどんな計算したいのかを言えば、
より有益なレスがつくのではないかな？ていうか小数なの？有理数なの？
ちなみに手もとのMathematicaブックには「0.142が近似値でなく
厳密な値であるならば、142/1000と有理数として入力しなければならない」
という記述をハケーン。
あと0.11111111111111111111`200 で200桁の有効数字が、
0.1111111111111111111111''200で小数点以下200桁までの正確さが保証されるんだと。
それで十分かな？どうよ？

>>656
> 　にしたところで何か計算したとたんにほぼ確実に10進数表示での
> 　数学的な等号は成立しなくなるんだぞ。
おいおいそれこそ「計算」次第だろ。アタマ疲れてない？
>>657
> 　おまえが書くプログラムってどれも実用には耐えないんだろうな。
CやF90で書くのとMmaで書くのとには違いが出て当然だろ。言語3～4種しか知らん人？

660:１３２人目の素数さん
03/08/22 22:35
>>659
ありがとう。

661:１３２人目の素数さん
03/08/22 23:36
こんなところでしかストレス発散できないやつ。

662:１３２人目の素数さん
03/08/22 23:55
>>659
つっこみどころ満載だけどおまえにはどんな反論しても
無駄みたいだからやめとく

663:１３２人目の素数さん
03/08/23 00:51
＞　つっこみどころ満載だけどおまえにはどんな反論しても
というと
＞　さあ香ばしくなってまいりました。これも5発売フィーバーか？
の辺りか？(w >659のカキコそのものはべつに間違った内容はないんじゃないの？

ていうか5出るのか。洩れのマスマティカブック4.2版どうなるんだろう。
5対応版買った人に～とか＋'の項目だけ写させてもらおうかな。
みなさんどうしますか？ちなみに洩れはサイトライセンスで使ってる人ね。

664:１３２人目の素数さん
03/08/23 00:53
なんでこんなにｲﾀｲんだろう・・・
自覚ないのか？

665:１３２人目の素数さん
03/08/23 01:17
>>663
大分前に出たかと

666:１３２人目の素数さん
03/08/23 02:48
Mathematicaを持ってない貧乏人がひがんでるスレってここでつか？（ｗ

667:１３２人目の素数さん
03/08/23 14:08
>>659
>> 　にしたところで何か計算したとたんにほぼ確実に10進数表示での
>> 　数学的な等号は成立しなくなるんだぞ。
> 　おいおいそれこそ「計算」次第だろ。アタマ疲れてない？

だからほぼ確実に、ってかいてんだよ。日本語ぐらい読めるようになろうね。
割算もできないしSinやCosも計算できないし、かけ算したらどんどん
桁が増えてくし、どうすんの？

> 　CやF90で書くのとMmaで書くのとには違いが出て当然だろ。言語3～4種しか知らん人？

Mathematicaの設計の話をしてんだよ。

>>660
0.1111111111111111111111111は>>656で書いてあるように自動的に任意精度と
してあつかわれている。
例えば、
x = 0.1111111111111111111111111
y = 1000 x - 111
Precision[x]
Precision[Sin[x]]
Precision[y]
とか入力すると、Mathematicaが精度保証して計算してるのがわかるよ。

高校生だからわからないのだろうが、>>644はアホだぞ。

668:１３２人目の素数さん
03/08/23 15:31
なんかこのへんのフレームは、任意精度と厳密数値(Schemeでいうexactnum)との
混同に基づいているような気がする。
ということで、高校生君に
> 　0.11111111111111111111111111
> 　って出て、それをクリックしてエンター押すと
> 　0.1111111111111111111111111099999999999964597
> 　って出るんですがどうゆう意味があるのでしょうか？かなり気になるんですが。
の理由を分かりやすく教えてあげれ > 誰か
多倍長整数や高精度近似数は2^16または2^32進数(整数要素の配列)らしい。
あとMathematicaBookの「3.1数」には
　123.4 機械精度
　123.4567890123456789 システムによっては任意精度 <- なんじゃそりゃ
　123.4567890123456789` 機械精度
　123.456`200 200桁の任意精度 <- precision
　123.456``200 200桁の任意確度 <- accuracy
　1.234*^6 機械精度
　1.234`200*^6 200桁の精度 <- 任意精度ととっていいのか？
とあるね。しかしこの例が微妙に分かり辛い...
リファレンス「入力の文法」にはdigits.digitsは「近似数」としか書かれてない。
これは任意精度なの？機械精度？だれか教えてくれますか？
わたしは普段整数しか扱わないので知らないままだった。

669:１３２人目の素数さん
03/08/23 16:23
0.111... を 1/2^n の和で作ってみると一番近い数が
そのゴミがついた数字になる、というのが一番簡単な
説明だよね。
まあいくつ項を取るかというのは気になる点だろうけどね。

670:１３２人目の素数さん
03/08/23 16:54
どうも高校生です。
1-0.999999999999=9.99978*10^-13
ってなるのが嫌なんです。これは分数で書くしかないですよね。
999999999999/1000000000000って書いたりするには桁数数えなければなりません。
かなり苦痛なので分数に書き換えるプログラムを書いてみました。
For[i = 0, x ≠ IntegerPart[x], x *= 10; i++]; Rationalize[x]/10^i
xに厳密な数値に直したい小数を入れて、実行すればOKです。

671:
03/08/23 19:37
ある演算をさせて結果がたとえば
If[x>0,・・・,～～～]
のようになったとき、・・・の部分だけをある変数に代入するにはどうすればいいのでしょうか？

672:１３２人目の素数さん
03/08/23 20:33
>>668
> 　 123.4567890123456789 システムによっては任意精度 <- なんじゃそりゃ
わかってて書いてるのかも知れないが、システムによって機械精度の桁が違うから。
計算速度のために機械精度で扱える時はあつかおうとするという
Mathematicaの設計の結果。

> 　リファレンス「入力の文法」にはdigits.digitsは「近似数」としか書かれてない。
digits.digitsって？

>>670
そのプログラムバグってるぞ。
x = 0.9999999999999999999999999
で試してみ。
# というかRationalize[x,0]の不完全版を回りくどくやっているだけ。
欲しいのは
FromDigits[RealDigits[x], 10]
あたりか？

>>671
状況が良くわからないが、
If[x > 0, ..., ...][[2]]
とか、もしかしたら
Simplify[If[x > 0, ..., ...], x > 0]
とか。

673:１３２人目の素数さん
03/08/23 20:36
>>668
わかった > 　digits.digits
機械精度にしろ任意精度にしろ近似数、ということ。
# 内部表現が2進なら必然的にそうなるな。

次ページ