Mathematicaプログラミング　質問箱　その1

MATLABの質問スレがあるのにMathematicaの質問スレがなかったので勃てました

本家
http://www.wolfram.com/index.ja.html
http://www.jip.co.jp/si/soft/mathematica/
Mathematicaユーザー会（有料）
http://www.jp-mathusers.org/

終了

もってねーYO!

高杉

海外逝き

単に数式組み合わせるだけだからプログラミングなんていう
ほどのものでは無いよな。まともに数値計算したい香具師なら
Cでゴリゴリやるし、記号計算ならそれこそ何を計算したいのか
どういう式を使えばいいかが分かっているはずだからな。
つまり、>>1はそんなもんに頼る前に数学でももっと勉強しろ。

Mathmaticaのスレなら数学板にあるけど

>>7
ここはMathematicaではなくMathematicaプログラミング

またmathematicaスレですか。。。。
どうせコンピュータで数式処理できることで何かちょっと勘違いしちゃった
卒研生かあほな院生が立てたｽﾚだろうな〜。研究室に一冊くらい凸版の関連本が
あるでしょ。その程度のもの見て理解できなければプログラミング以前に学力の
問題ですよ。マニュアルもカナーリ詳しい。mathematicaのプログラミングって
言っても、紙に数式書くのとほとんど同じなのに一体何を質問したいんですか？
アルゴリズムとか数値計算ですか？残念ながらそういう用途には向いていませんよ。
確かに使いやすいですがmathematicaじゃトロくて時間がかかりすぎるので
本気でやるならC/C++あたりを使った方がいいですよ。

実は英語の数学用語が分からなくてスレを立てた、に一票

それで>>1の質問は何なの？
何を質問したくてスレ作った？
板ちがいだが、分かる範囲で答えてやるぞ。
ただし、最新版は使ってない。それでよければいいよ。

Mathematica使いこなしてるようなのは、この板にいないだろ？
方向性が狭い

つーか、mathematicaは単なる数学ソフトなんだが。

>>13
プログラミングの為の構文が用意されているじゃないか。
グラフを表示したりするのもわりと楽だし、シミュレーションしたりするのにも
いちいちCで書いていたら手間が掛かりすぎる。

mathematica programing nb

で検索し、援助交際　叶う筈　　( *.nb Getし）

いちいちCで書いていたら人生無駄。

在るもの　*.nb は　美味しく　イタダク　べし

　高　杉

mathematica programing nb

で検索し、援助交際　叶う筈　　( *.nb Getし）

いちいちCで書いていたら人生無駄。

在るもの　*.nb は　美味しく　イタダク　べし


追記；　絞りきれないなら　目的を　明確にして下さい

参考　.nb を　紹介します

mathematicaで出来るシュミレーションなんて小規模なもの。
何をシュミレーションしたいんだよ？はっきり言ってみろ。
グラフを表示したりするのはそこいらの表計算ソフトでも出来るぞ。
mathematicaのプログラミングの為の構文なんてたかが知れてる。
mathematicaが単なる数学ソフトであることに気が付いていない
時点でアウト。研究者には向かないな。もうageるな。

>>19
お前面白いな。

mathematicaなんていうのは、一般的な数学計算を対象にしているから
シュミレーションのように極度に特化したものを本気でやりたいなら
専用のプログラムをCかなんかで作って行うべし。でないと、
スピードも遅いし、シュミレーションの検証もまともに出来ない。
まあ、結果を整理する時には便利かもしれないが。

>>19
見ていてちょっと恥ずかしい。

>>21
君が考えているシミュレーションってなんか偏ってるよね。

何をシュミレーションしたい？ 明言して！
@@@@@@@＠＠＠＠＠＠＠＠＠＠＠＠＠＠＠

　　　　　　　　　　　研究者　様

　　　　　　　　スグ悪意にトル　癖　ヤメテ

Mathematicaが買えないからうらやましいんだよ、きっと。

繰り返し計算してみるとかなり遅いのはわかるけど、
そういうことはただ向いていないだけなんだよー　にゃんにゃん

繰り返し計算してみるとかなり遅いのはわかるけど
　　　　とありますが、　私マツーワ　　（数秒間です）で

　　　　　　以下から　Getした　*.nb は　叶う。

　　　　　　　　　超繰り返し計算の例ですが　....
　　　　　　　　　@@@＠＠＠＠＠＠＠＠＠＠＠＠＠＠＠

http://www.diku.dk/forskning/image/teaching/Courses/e01.426/lectures.html

Getされて　試され　Cかなんかで作った　のと

比較　備忘録　を　此処に　記述　して　下さい。

　　　　　　　　　研究者　様

たとえばねー　う〜ん…
FIRフィルタを音楽全体にかけたりとかすると違うと思うですー　にゃんにゃん

>>27
日本語変だよー　どうしたの？　まさか…？！

日本語変だよー　どうしたの？　まさか…？！

　　　　math かっ!!!!!!!!!!! でありmath.


正鵠であり　math



　　　normalなのに　日本語変?????

２ｃｈ　初邂逅　初夜　で　ありmath.


　　　　ｗｗｗで通用する　言語所望。

　日付けカワリ；

http://www.wolfram.com/index.ja.html


ｗｗｗ　で　非アブnormal な　言語で
\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
記述　在りmath.

眉唾か　懐疑　顛末　録　まで　長時間　を　要し math.

　新規購入された方　　　(研究者　非　コピ　ぺ）

　　　　　　　　　　　　　　は　　その後

　　　　　　　　　　オシエテ　無用で　あったか　を!
　

>>31
お前が無用だよ。まったく。

お前が無用だよ　って　貴方　何様。

眼には眼　歯には　歯　

眼＿＿i＿＿＞眼


http://www.google.com/search?sourceid=navclient&ie=UTF-8&q=%E7%84%A1%E7%94%A8%E3%81%AE%E7%94%A8

埴輪ハオ！

FIRフィルタ

eg;


http://216.239.63.104/search?q=cache:PFXSaxwx2EgJ:www.additive-net.de/software/mathematica/packages/mma.packages.digimageproc.shtml+FIR+filter++mathematica+nb&hl=ja

>>35
そういうのがあるのは知っているけれど、結局はCでの実装でしょ？
そのループ部分はCでやっているもん。
ループする部分をMathematicaの構文でやらないとダメですー。

実は微分積分忘れました！

実は結構な項書き換えに基づく
関数型言語だったりするんだよなーこれ。
でも高すぎ。誰かクローン作ってくれないものかな。

>>38
クローンとは違うけれど、

Risa/Asir
http://risa.cs.ehime-u.ac.jp/index.html

>>39
おお、日本製でこんなのがあったのですか。
情報ありがとうございます。
こういうのって言語機能自体よりライブラリとか作るほうが
大変なんだろうなあ。

Gröbner Bases　mathematica nb

で検索し　学界の先達　の

知恵を　有難く　拝借（*nb Get)

　　　　　　　　　知恵熱が∂まで

(Risa/Asir)

>>38
評価版で全機能が使えるわけだけど。

中級者向けの本を教えてください

>>43
「The Mathematica Book」 っていう分厚い本なかったっけ？
あれひととおり読めばだいたいのこと分かると思うけど。
http://www.amazon.co.jp/exec/obidos/ASIN/1579550223/249-8544534-2312305

ituneやテレ朝CIでやられとるの音声のヴィジュアライズって何をつかって作っておるのでしょうか？
ご存知の方をいらっしゃればお知恵を...;;

>>45
スレ違い

カレントディレクトリを簡単に変える方法はないのですか？

楕円のシュミレーションなんですけど、楕円の回転のさせ方おしえてください。

動画を作ってファイルに保存する方法を教えてください

おれはmaximaでことたりるな。無料だし。

>> 1
ついでにMapleのスレも立ててくれ。反応眺めるだけでおもろかったりする。
スレの定着っぷり見るとやっぱMATLABがメジャーなんかね。

mapleってつかえないよな

はなはだおこがましいとは思いますが、URL:
http://www.vector.org.uk/v201/cliff201.htm#note4
にある下記のソースは、Ｃのコードでいうとどういうことになるか
書き換えできないでしょうか。（入門のURLを見ましたが、'|','$','@'
などの表記が分からず、ここへ来ました）

↓行番号（書きこみに当たり追記）
１　NB. c barrel_undistort raw_image
２　barrel_undistort=: 4 : 0
３　wh=.|.hw=.2{.$y.
４　und=.wh sn_to_pix@(x. sn_undist)@(wh pix_to_sn)
５　y.=.(y.,{:a.),"_1 {:a.
６　if. x.<0 do.
７　'dw dh'=.|und 0 0
８　else.
９　'dw dh'=.-<./und -:wh*=i.2
10　end.

（>53 の続き）
11　z=.((hw+2*dh,dw),3)${:a.
12　dsn=. -{.-/hw pix_to_sn 0 0,:1 1
13　w=.j./"1 hw pix_to_sn (-dh),.(-dw)+i.(2*dw)+{:hw
14　for_k. (-dh)+i.({.hw)+2*dh do.
15　 t=.12 o. w
16　 r=.|10 o. w
17　 r=.x. R_to_RD r
18　 in=.;/_1 >. hw <."1 hw sn_to_pix +.r r. t
19　 z=.(in{y.)(dh+k)}z
20　 w=.w+dsn
21　 end.
22　NB. 'one'+1
23　z
24　)

Mathematicaでc並のスピードは出せますか？

mahtematicaまともに使える奴はここにはおらんか・・・

普通の人は数式処理なんてしませんから。

>>55
C並とまでは行きませんが、ある程度のスピードは出ます。

mathematicaとSのそれぞれの特徴(違い)を教えてください。

下記の一見ごく普通の二重フーリエ積分がMathematicaでどうしても実行できま
せん。

Integrate[(Cos[omega1 T1]+I Sin[omega1 T1]) (Cos[omega2 T2]+I Sin[omega2
T2]) (2.0 Exp[-0.5 (T1+2 T2)] Cos[2500. T1]+Exp[-0.5 T1-0.5 T2] Cos[2500.
(T1-T2)]),{T1,0,Infinity},{T2,0,Infinity},Assumptions->{omega1>0,omega2>0}]

Mathematica 5.0のWindows版で実行すると、何のエラーメッセージも無くMathKernel
そのものが落ちます。（それまでの計算履歴も全て消失します）

Mathematica 5.0のUnix版で実行すると、Segmentation fault (core dumped)という
エラーメッセージとともにMathematicaが落ちます。

ちなみに上記プログラムのひとつ目の指数関数の前の「2.0」の数値定数を「2」にす
るだけで上記積分は実行できるようになります。

なお、上記積分を実行中のワークステーションのメモリをチェックしていると、計算
が止まる前後でメモリはわずかに上昇するものの、ワークステーションのメモリ限界
値の半分も使用していない状態です。（実際、"No more memory available"などのエ
ラーは一切表示されません。）

こんな症状は初めてで原因が全く分からず戸惑っており、もしこのような症状に出く
わしたことがある方がおられたなら、対処法をご教授いただければ幸いです。

何卒よろしくお願いいたします。

>>60
そういうこと、よくあるよ。俺も簡単なDSolveが落ちるって苦情
出したことあるし。でも、「新しいバージョンでは直っていますので
買ってください」としか返ってこないけど(w

マゾ？

Mathematica 5.1では、DendrogramPlot[]が利用できるようになりましたが、
発現解析等で良く利用される2軸(?)のDendrogramを簡単に描くことは可能でしょ
うか。
なにかヒントがあればご教示下さい。

http://acykhm035074.adsl.ppp.infoweb.ne.jp/
うぇｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗっうぇっうぇおｋｗｗｗ
ｗｗｗっｗｗｗｗｗｗｗｗｗｗｗｗおｋｗｗｗ
うぇｗｗｗ
うぇｗｗｗおｋｗｗｗっうぇｗｗｗｗｗｗｗｗｗｗｗｗｗ

CUIで使ってるときに、1行に表示される文字数を変更する方法はありませんか？
例えば、

In[1]:= N[Pi,100]
Out[1]= 3.1415926535897932384626433832795028841971693993751058209749445923078\
> 16406286208998628034825342117068

としたときに、結果を改行せずに表示してほしいのですが。

http://218.33.212.67.eo.eaccess.ne.jp/
おｋｗｗｗおｋｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗ
ｗｗｗｗｗｗｗｗｗｗｗｗｗおｋｗｗｗうぇｗｗｗ
ｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗ

hoshu

wwｗwwｗｗｗっうぇｗｗｗｗｗｗｗｗｗｗｗｗっうぇｗｗｗｗｗｗｗｗｗｗっうぇｗｗ

入門のサイトを教えてください

>>69
正直、付いてくるチュートリアルを1つずつ順番に読んでいくのが
一番勉強になる。

age

Mathematicaについての質問なんですけど，
Mathematicaで作った***.nbのファイルをコマンドプロンプトから実行させたいんですけど，
どのようにコマンドラインを入力すればわからないのでどなたは教えてくれませんか？

どなたは教えてくれません。

ちね

ImageProcessing最高!!

ﾎｼｭｰﾘ

保守すんな

新スレたてました！

フーリエ変換が趣味です。
http://hobby8.2ch.net/test/read.cgi/hobby/1133005088/

6.0ﾏﾀﾞｰ?

ここは僕の待ち望んだスレです

ファイルのデータを読み込んで、その全てを足すコマンドってないのでしょうか？

作れ自分で

マスエムエルって凄いの？

MATLAB>>Mathematica>Mathcad

ここはMathematicaクローンを作ってフリーソフトにするスレですか？

そりゃ無理でしょ

ダウンロードしたファイルを表示させるには？

こんにちは、

下記をダウンロードして評価させたのですが、音がするだけで、結果が表示しません。
どうしたら、表示するのでしょうか？windows2000,
mathematica4.1を使用しております。

http://library.wolfram.com/infocenter/Articles/2446/#downloads

NuclideView.nb

多入力多ルールファジィ推論プログラム　が全くわかりません。
どなたか分かりますかね？

S = {2, 5, 9, 22, 89, 105, 200};
a = 105;
h = 1; t = Length[S];
k = Floor[(h + t)/2];
While[h <= t, If[S[[k]] == a, Break[],
If[S[[k]] > a, t = k - 1, h = k + 1]];
k = Floor[(h + t)/2]; ];
If[h > t, k = -1];
Print[a, " ", k];

なにを求めるプログラムなのかわかんない…
だれか教えてください。

>>89
配列Sから要素aを探索するプログラムと思われ
試しにaを別な値に書き換えて実行してみ
見つかった場合にはその位置インデクス
見つからない場合には-1が返る

>>90
ありがとうございます
やってみます

たびたびすみません
どういう構造になってるんですか？

長さlの配列Sが昇順に整列されているとする。

Sの中の要素aを探すには、
まず中央の要素S(m)とaを比べる。a＝S(m)ならそこで検索終了。
もしS(m)＞aならば、S(n)＝a (n≧m) を満たすnは存在しないので、S(n) (0≦n＜m) を探せばいいことになるので、
S'(n) = S(n) (n=1,2,3…m-1) とする。S(m)＜aの場合も同様にして、
S'(n) = S(m+n) (n=1,2,3…l-m) とする。
S'(n) についても上記のことを繰り返す。

>>92
Cで書き換えるとこんな感じ（インデント全角空白注意）
int main() {
　int S[7] = { 2, 5, 9, 22, 89, 105, 200 };
　int a = 105;
　int h = 0;
　int t = 6;
　int k = (h + t) / 2;
　while (h <= t) {
　　if(S[k] == a) break;
　　else if(S[k] > a) t = k - 1;
　　else h = k + 1;
　　k = (h + t) / 2;
　}
　if (h > t) k = -1;
　printf("%d %d¥n", a, k);
　return 0;
}

>>92
Javaならこうな（インデント全角空白注意）
public class Sample {
　public static void main(String[] args) {
　　int[] S = new int[] {2, 5, 9, 22, 89, 105, 200};
　　int a = 105;
　　int h = 0;
　　int t = S.length - 1;
　　int k = (h + t) / 2;

　　while (h <= t) {
　　　if (S[k] == a) break;
　　　else if(S[k] > a) t = k - 1;
　　　else h = k + 1;
　　　k = (h + t) / 2;
　　}
　　if (h > t) k = -1;
　　System.out.println(a + " " + k);
　}
}

(^_^;)

ありがとうございます

ちゅうかMathematicaなんだから、用意されてる関数を使おうよ。
Positionで一発。

>>98
そりゃそうだな
おれ処理系ないけど

エキスポートなどで作られるファイルを
自分の指定するフォルダにすることってできますか？

Export["c:/hogehoge/...
か、SetDirectory

データをTSV形式でExportすると、
599.8が599.7999999999999のようになってしまいます。
このようにならないようにすることはできないのでしょうか。

>>102
コンピュータってのは2進数で動いているので、
そういう風にできているんです。

ListPlotやPlotで図を表示させないようにしたいんですけど、どうすればいいでしょうか？

>>104
DisplayFunction -> Identity

>>105
できました。
ありがとうございます。

数式を特定の変数で整理するのはどうしたら良いのでしょうか？
よろしくお願いいたします

>>107
Collect[]

>>108
できました！ありがとうございました！

初めまして、


下記について、ご教示頂きましたら幸いです。

http://www.geocities.jp/dirac_equation/bohr.pdf
（注意；そのままクリックしても表示しませんので、お手数ですがコピーして貼り付けてください。）

の式（６A-1）を式（６A-21）に代入して計算したのですが、答えが一致しません。どこが悪いのでしょうか？
Clear[q];
Clear[f];
Clear[k];
Clear[Rp];
(*式(6A-1)*)
Rp[l_,m_]:=R0*(Sum[a[l,m]*(-1)^m*SphericalHarmonicY[l,-m,q,f],{l,0,l},{m,0,m}]);
l=0;
m=0;
(*式(6A-22)*)
f0=3/(4*Pi)*Z*e/y;
f0=1/r;
dr=3/(4*Pi)*Z*e/R0^2*(4*Pi)^(1/2)*a0*Integrate[(DiracDelta[r-R0]+1/2*R0*D[DiracDelta[r-R0],r])*f0*r^2,{r,R0,0}]//FullSimplify

クーロンエネルギーのmathematicaプログラム

原子核（液滴模型）のクーロンエネルギーをmathematicaで計算したいのですが、既に計算したプログラムは無いでしょうか？ありましたらご教示頂きましたら幸いです。

また、自分でプログラムを作る場合、下記のどちらを元に作る方が良いでしょうか？

http://www.geocities.jp/dirac_equation/bohr.pdf
の6A-1f Coulombエネルギー

もしくは
http://www.geocities.jp/dirac_equation/nogami.pdf
の変形による静電エネルギーの変化

Mathematicaプログラムを有料で作る計算屋さんについて

下記のプログラムを作りたいのですが、自分では出来ないので困っています。
そこで、プロに有料で頼みたいのですが、計算屋さんというか？そんなプロって日本に存在するのでしょうか？

お願いしたいプログラム

下記HPのプログラム

http://www.geocities.jp/dirac_equation/bohr.pdf
の6A-1eの表面エネルギーと6A-1f Coulombエネルギー

と

http://www.geocities.jp/dirac_equation/nogami.pdf
の表面エネルギーの増加と、変形による静電エネルギーの変化

>>112
原子核の液滴模型か

俺は一応、大学で量子力学を専攻したんだけど、
金額次第では作っても良いけど？

>>113


>俺は一応、大学で量子力学を専攻したんだけど、
>金額次第では作っても良いけど？

プログラムは、だんだん出来そうな感じになってきました。
すいません。もうしばらくして、やはり無理な場合は、考えます。

質問失礼します。

１から１０までの数字を重複せずランダムで１０個取り出すプログラムは
どう書けばよろしいでしょうか？
[3.2.7.5.6.4.10.9.1.8]
↑のように取れればいいんですが、必ず数字が重複してしまって・・・

よろしくお願いします。

>>115
超初心者ですが一応目的は達成できていると思います

Clear[a, b, c]; b = {};
For[i = 1, Length[b] < 10, i++, a[i] = Random[Integer, {1, 10}];
c = Position[b, a[i]];
If[c == {}, AppendTo[b, a[i]]];
]

bが１〜１０のランダムなリストになってるはずです

>>115
連続ですみません
ヘルプブラウザのRandomの説明の他の例の一番下にスマートなやり方が乗ってましたよ

>>116と同様の手法ですが，関数型言語っぽく，関数の再帰的定義を用いて
重複のない乱数列生成関数を作ってみました。

Clear[RandomList]

RandomList[start_Integer, end_Integer, list_:{}] :=
  If[Length[list] > end - start, list,
    RandomList[start, end,
      (If[MemberQ[#1, #2], #1, Append[#1, #2]] &)
         [list, Random[Integer, {start, end}]]
    ]]

RandomList[1, 10] (* 1〜10 の乱数順列生成 *)
RandomList[1, 10, {2, 5}] (* 初項が {2, 5} で始まる 1〜10 の乱数順列生成 *)

任意の位で切り上げる関数が見つかりません
私が見つけられないだけで、きちんと用意されているのでしょうか？
おしえてください

>>119
やりたいのはこういうこと？
ceiling[num_Real, order_Integer] := 10^order Ceiling[num 10^(-order)]
N[ceiling[Pi, -3]] (* π の値を 10^(-3) オーダーの桁を切り上げて表示 *)

>>120
まさしくその通りです
参考にさせていただきます
ありがとうございました

FindRootで表示桁数を上げる方法を教えてください

>>122
WorkingPrecision→100

＞大学で量子力学を専攻したんだけど、

え？

compute algebraについて参考になる日本語の本ってありますか？

Pari/GPという計算ソフトを用いて作られたプログラムをMathematicaに変換できる人いませんでしょうか？
先日尋ねられたのですが、わたしはどちらも詳しくないので教えて頂きたいと思います。
元のプログラムは、、、

{
n=6;
g=matrix(n,n);
gs=matrix(n,n);
m=matrix(n,n);
u=matrix(n,n);
h=matrix(n,n);
c=0;
"====================step1====================";
g[1,]=[1,0,0,0,0,-21];
g[2,]=[0,1,0,0,0,-5];
g[3,]=[0,0,1,0,0,-10];
g[4,]=[0,0,0,1,0,-20];
g[5,]=[0,0,0,0,1,-24];
g[6,]=[0,0,0,0,0,49];
print("--------------------step1--------------------");
print("G=");
for(s=1,n,print(g[s,]));
print("");
print("M=");
for(s=1,n,print(m[s,]));
print("");
print("GS");
for(i=1,n,print(gs[i,]));
print("M*GS=");
mgs=m*gs;
for(s=1,n,print(msg[s,]));
print("");print("");

"====================step2====================";
gs[1,]=g[1,];
u[1,1]=1;
for(i=2,n,gs[i,]=g[i,];
for(j=1,i-1,uu=g[i,]*gs[j,]~;us=gs[j,]*gs[j,]~;u[i,j]=uu/us;
gs[i,]=gs[i,]-u[i,j]*gs[j,]
);
u[i,i]=1;
);
m=u;

print("--------------------step2"--------------------");
print("G=");
for(s=1,n,print(g[s,]));
print("");
print("M=");
for(s=1,n,print(m[s,]));
print("");
print("GS=");
for(i=1,n,print(gs[i,]));
print("M*GS=");
mgs*m*gs;
for(s=1,n,print(mgs[s,]));
print("");

i=2;
print("i=",2);
print("");
print("");
u=m;

"====================step3====================";
while(i<=n,
j=i-1;
while(j>0,
g[i,]=g[i,]-round(u[i,j])*g[j,]);
m[i,]=m[i,]-round(u[i,j])*m[j,]);
u[i,]=m[i,];
j=j-1;
print("--------------------step3"--------------------");
print("G=");
for(s=1,n,print(g[s,]));
print("");
print("M=");
for(s=1,n,print(m[s,]));
print("");
print("GS=");
for(i=1,n,print(gs[i,]));
print("");
print("M*GS=");
mgs=m*gs;
for(s=1,n,print(mgs[s,]));
print("");print("");
);

if((i>1&&(gs[i-1,]*gs[i-1,]~>2*gs[i,]*gs[i,]~),
h[i,]=g[i,];g[i,]=g[i-1,];g[i-1,]=h[i,];
mi=i-1;
for(ii=mi,n,
gs[ii,]=g[ii,];
for(j=1,ii-1,
uu=g[ii,]*gs[j,]~;us=gs[j,]*gs[j,]~;u[ii,j]=uu/us;
gs[ii,]=gs[ii,]-u[ii,j]*gs[j,]
);
u[ii,ii]=1;
);
i=i-1,
i=i+1
);
m=u;
print("--------------------step3(b)--------------------");
print("G=");
for(s=1,n,print(g[s,]));
print("");
print("M=");
for(s=1,n,print(m[s,]));
print("");
print("GS=");

for(i=1,n,print(gs[i,]));
print("");
print("M*GS=");
mgs=m*gs;
for(s=1,n,print(mgs[s,]));
print("");
print("i=",i);
print("");print("");

);

print("====================Finish====================");
print("G=");
for(s=1,n,print(g[s,]));
print("==============================================");

}

ちょっと長いですがこういうプログラムです。
どうかお願い致します！

半径0.5の円で以下の図を書くプログラムは？

○●○●○●○●○
●○●○●○●○●
○●○●○●○●○
●○●○●○●○●
○●○●○●○●○
●○●○●○●○●
○●○●○●○●○
●○●○●○●○●
○●○●○●○●○

半径1の円の軌道上に等間隔で0.1の円を40個起き、置かれた円が
・・・○○○●○○○●○○○●・・・になる図は？



上は
Black = Disk[{i,j},0.5];
White = Circle[{i,j},0.5];
で||か&&を使う？？
下は完全に謎です

どなたかよろしくお願いします

ヒント：
sin(a) cos(a)

さ、サインコサインですか・・・

もう忘れちゃってるやorz
ちょっとこれから学校に行って再チャレンジしてきます


ダメそうだったらまた質問にきます

>>131
処理系ないけど過去レスから文法を判断するとこんな感じか？

Circle[{x, y}, 1];
a = Pi / 20;
For [i = 1, i <= 40, i++,
t = a * i;
x = Cos[t];
y = Sin[t];
If ((i % 4) == 0) Disk[{x, y}, 0.1];
Else Circle[{x, y}, 0.1]
];

mathematicaメーリングリストに入会したいので、HPに記載の通りやったのですが、下記メールが１ヶ月前に来てから何の返事も無く、入会できません。どうしたら、入会できるのでしょうか？また、クレーム等はどこに連絡すればよいのでしょうか？



Good. Your subscribe request is confirmed and forwarded to the maintainer.
He/She will subscribe you, so PLEASE WAIT A LITTLE.

あなたからのメーリングリスト登録意志を確認しました。
登録依頼はメーリングリストの管理者へフォワードされました。
少し待ってみて下さい。

131です
ヒントを元に>131の二個ある内の下のプログラムは組むことができました
色々試したけど上のプログラムは結局できませんでした

どなたか上のプログラムを教えてください、よろしくお願いします

ヒント：2重ループ
てか上の奴の方は簡単すぎだろ

今数種類の微分方程式の解の最大値を求めるプログラムを組んでいます
For[i = -5, i < 6, i = i + 1,
NDSolve[{y''[x] + y'[x] + y[x] ==
1 - Cos[( 2 Pi i )x], y[0] == 0, y'[0] == 0},
y, {x, 0, 1}, MaxSteps -> 10^10]
data[[i]] = Max[Table[y[x], {x, 0, 1, 1/10^6}]]

]
これでdataに最大値の近似値が格納されると思ったのですが…
どなたか解決法をご教授下さい

２次元プロットの質問です。

たとえば　f=a1*Sin[x]+a2*Cos[x]
の式で、a1がa2より大きい場合とa2がa1より大きい場合を調べたい時、
どうプロットすれば良いのですか？
a1,a2に適当な数値を当てはめるしかないですか？

なんかmathematicaがようやくversion 6になったんだって? ﾌﾟｯ、おくれてんなー
matlabが6になったの何年前だと思ってんだよ、おめでてーな(ﾜﾗ

おめでたいな

デフォルトの文字サイズ設定するのってどうやる？

>>142
デフォルトのスタイルシートの変更
http://support.wolfram.com/mathematica/interface/options/stylesheet.html

mathematicaでアボットの負荷曲線を描きたいのですが、どうしても描けません！
誰か教えてください。

mathimaticaで出力されるグラフの軸の目盛りを外向きにする方法を、
ご存知の方が居られましたらご教授いただけないでしょうか。

何卒よろしくお願いします。

>>145
オプションでTicksを一つ一つ設定するしかないかも。

Ticks ->
　{
　　{{1, 1, {0, .02}}}, Automatic
　}

みたいな感じで。
もっと簡単な方法があるのだろうか？

TicksあるいはFrameTicksに関数を指定できるよ。

Mathematica式の挿入について


こんにちは、
下記のように、ｙを微分して得られる解をまたｙに代入したいのですが


y = (2/5)*(1 - x)*a^2 - (4/105)*(1 + 2*x)*a^3;
y /. a -> Solve[D[y, a] == 0, a][[2]]


と計算してもうまく計算できません。
どうしたら良いでしょうか？

要するに、下記にように

Simplify[y /. a -> -((7*(-1 + x))/(1 + 2*x))]

計算して

-((98*(-1 + x)^3)/(15*(1 + 2*x)^2))

を得たいのです。

>>148

y /. Solve[D[y, a] == 0, a][[2]]

>>146,147

お答えありがとうございます。
やって見ます。

こんにちは、
下記は、mathematicaで亜鈴の図を、平面図と立体図で書いたものです。
平面図のように、立体図を横にする方法を教えてください。


<< "Graphics`ParametricPlot3D`"
ParametricPlot[{{r1 = 1 + 1*LegendreP[2, Cos[q]]; r1*Cos[q], r1*Sin[q]}}, {q, 0, 2*Pi}, AspectRatio -> Automatic];
r2[n_, thet_] := LegendreP[2, Cos[thet]]
Do[SphericalPlot3D[1 + r2[2, thet], {thet, 0, Pi, Pi/30}, {phi, 0, 2*Pi, Pi/15}]]

Mathematicaで強制振動を扱っているHP

こんにちは、

Mathematicaで強制振動を扱っているHPを探しています。一定の力を加えたら、どんどん振幅が大きくなるっていくような図はないでしょうか？

>>151
<<"Graphics`ParametricPlot3D`"
SphericalPlot3D[1 + LegendreP[2,Cos[thet]],
{thet, 0, Pi, Pi/30},{phi, 0, 2*Pi, Pi/15},
ViewPoint->{-3.382, 0.096, 0.076},
ViewVertical->{0,1,0}, ImageSize->300]

>>153
お返事ありがとうございます。
ずーと、出来なくて困っていました。
助かりました、今後ともよろしくお願い致します。

>>149

お返事ありがとうございます。

わかりました。

10年前にこのスレ欲しかったぜ

Mathematica何が悪いのでしょうか？

こんにちは、

下記は、HPをそのまま写したのですが、評価するとエラーになります。
どこが間違っているのでしょうか？

http://www.rs.kagu.tus.ac.jp/~phlabex/LabExercise/LCR/frame/lcr2.htm


L := 2.65/10^3;
r := 73.32;
c := 12.55/10^9;
e1 := 5;
i1 = y[x] /. NDSolve[{Derivative[2][y][x] == (-(r/L))*Derivative[1][y][x] - (1/(L*c))*y[x] + e1/L, y[0] == 0.0001, Derivative[1][y][0] == 0}, y, {x, 0, 0.0001}]
Plot[i, {x, 0, 0.0001}]
q1 = D[i, x]
q2 = plot[g1, {x, 0, 0.0001}]

追伸

このグラフは、ｘが増加するに連れて、減少してますが、増加させるにはどこを変更すれば良いでしょうか？

>>157
i1 じゃなくて i
g1 じゃなくて q1

＞このグラフは、ｘが増加するに連れて、減少してますが、増加させるにはどこを変更すれば良いでしょうか？

Plot[-i, {x, 0, 0.0001}]
とでもすれば？

>>158

お返事ありがとうございます。

＞i1 じゃなくて i
＞g1 じゃなくて q1

了解いたしました。


＞Plot[-i, {x, 0, 0.0001}]
＞とでもすれば？

振幅の符号が逆転しましたが、やはり減少しております。

こんにちは、

下記プログラムで、円から亜鈴に変形する図を作成したのですが、このプログラムを書き換えて、アニメのように、ある場所で固定させて、時間の経過とともに、円から亜鈴、亜鈴から円、円から亜鈴、亜鈴から円に、繰り返し変形するようにすることは可能でしょうか？


For[t = 0, t <= 5, t++,
r0 = 1;
a1 = 0.2*t;
ParametricPlot[{{r1 = r0 + a1*LegendreP[2, Cos[q]]; r1*Cos[q], r1*Sin[q]}}, {q, 0, 2*Pi}, AspectRatio -> Automatic];
];

＞このグラフは、ｘが増加するに連れて、減少してますが、増加させるにはどこを変更すれば良いでしょうか？

何のグラフを描きたいのかがわかりません

お返事ありがとうございます。

＞何のグラフを描きたいのかがわかりません


説明不足ですいません。下記を評価しますと、ｘが増加するにつれて、振幅が減少しますが、ｘが増加するにつれて、振幅が増加するグラフを描きたいです。要するに、左から右に行くにつれて、揺れが大きくなるようにしたいです。


L := 2.65/10^3;
r := 73.32;
c := 12.55/10^9;
e1 := 5;
i = y[x] /. NDSolve[{Derivative[2][y][x] == (-(r/L))*Derivative[1][y][x] - (1/(L*c))*y[x] + e1/L, y[0] == 0.0001, Derivative[1][y][0] == 0}, y, {x, 0, 0.0001}]
Plot[i, {x, 0, 0.0001}]
q1 = D[i, x]
q2 = Plot[q1, {x, 0, 0.0001}]

Plot[-i, {x, 0, 0.0001}]

それはわかるんだけど、だったら単に x*sin[x] でもいいのかね？

お返事ありがとうございます。

下記のように、((r/L))*Derivative[1][y][x]の符号を、マイナスからプラスにしましたら、振幅が増加するようになりました。これが、減衰振動の反対の現象だと考えてよろしいでしょうか？

L := 2.65/10^3;
r := 73.32;
c := 12.55/10^9;
e1 := 5;
i = y[x] /. NDSolve[{Derivative[2][y][x] == ((r/L))*Derivative[1][y][x] - (1/(L*c))*y[x] + e1/L, y[0] == 0.0001, Derivative[1][y][0] == 0}, y, {x, 0, 0.0001}]
Plot[i, {x, 0, 0.0001}]
q1 = D[i, x]
q2 = Plot[q1, {x, 0, 0.0001}]

>>160
Animate Selected Graphics (MacだとCommand-Y)?

方程式の答えを図に反映させたいのですが、上手くいきません。
x = 0.7;
a2 =. ;
y = a2 + (7*(-1 + x))/(1 + 2*x);
Solve[y == 0, a2]
答えは
a2=
となります。
これを、ParametricPlot[{{r1 = 1 + a2*LegendreP[2, Cos[q]]; r1*Cos[q], r1*Sin[q]}}, {q, 0, 2*Pi}, AspectRatio -> Automatic];のa2に代入して、図を書きたいのですが、どうしたら書けるでしょうか？


x = 0.7;
a2 =. ;
y = a2 + (7*(-1 + x))/(1 + 2*x);
Solve[y == 0, a2]
ParametricPlot[{{r1 = 1 + a2*LegendreP[2, Cos[q]]; r1*Cos[q], r1*Sin[q]}}, {q, 0, 2*Pi}, AspectRatio -> Automatic];

>>166
x = 0.7;
a2 =.;
y = a2 + (7*(-1 + x))/(1 + 2*x);
rule1 = Solve[y == 0, a2][[1]];
ParametricPlot[{{r1 = 1 + a2*LegendreP[2, Cos[q]] /. rule1; r1*Cos[q],
r1*Sin[q]}}, {q, 0, 2*Pi}, AspectRatio -> Automatic];

>>167

お返事ありがとうございます。

わかりました。

申し訳ございません。更に、下記でa2の数値を表示させるには、どうしたら良いでしょうか？


x = 0.5;
a2 =.;
y = a2 + (7*(-1 + x))/(1 + 2*x);
rule1 = Solve[y == 0, a2][[1]];
Print["x=",x];
Print["a2=",a2];

>>168
x = 0.5;
a2 =.;
y = a2 + (7*(-1 + x))/(1 + 2*x);
rule1 = Solve[y == 0, a2][[1]];
Print["x=", x];
Print["a2=", a2 /. rule1];

>>169
いつも、お返事ありがとうございます。
＞Print["a2=", a2 /. rule1];
わかりました。
すいませんが、更に下記を教えて下さい。

次のHPは、核分裂の臨界エネルギーを求めたものです。
http://space.geocities.jp/gamow_gold/nuc.pdf
（上記をコピーして、アドレスに貼り付けて、エンターしてください）

質問
P6のFig.4.5の図を得るため、 (30) 式の下の a2=7(1-x)/(1+2x)　という式から a2 を求め、後はこの値を（１７）式に代入し、mathematicaで作図したのが下記です。しかし、評価してみますと、Fig.4.5の図と形が異なります。
（１７）式を、(1 + a2*LegendreP[2, Cos[q]])としたのが、多分、間違っていると思われますが、（１７）式をどのように修正すれば、P6のFig.4.5の図が得られるでしょうか？




For[n=0,n<=10,n++,
x = 0.1*n;
a2 =.;
y = a2 + (7*(-1 + x))/(1 + 2*x);
rule1 = Solve[y == 0, a2][[1]];
Print["x=",x];
Print[rule1];
ParametricPlot[{{r1 = (1 + a2*LegendreP[2, Cos[q]]) /. rule1; r1*Cos[q],
r1*Sin[q]}}, {q, 0, 2*Pi}, AspectRatio -> Automatic];
];

こんにちは、

Mathematicaで作成した図を、LaTeXに貼り付ける方法および手順を教えてください。例えば、下記で作成した図について教えてください。

{r0,r1,n,qs}={1,0.1,8,Pi/32};
ParametricPlot[{{Cos[q],Sin[q]},{r=r0+r1 Sin[n q];r Cos[q-qs],r Sin[q-qs]}},{q,0,2Pi},AspectRatio®1];

抵抗とコンデンサからなる並列回路の複素インピーダンスとその絶対値を求めるプログラムを作成せよ。
プログラムは、抵抗の値、コンデンサの容量、周波数を入力すると、複素インピーダンスとその絶対値を求めその結果を表示するものである。
このプログラミングはどのようになりますか？おねがいします

おまえら2ちゃん見てる暇があるならせめてググれ
もしくはチュートリアルくらい読め

＞＞173
171了解。ところで、「ググれ」の意味は何でしょうか？

こんにちは、
以前このサイトでご教示頂きました下記の図形なのですが、横軸の数字が重なって気になるので、いろいろとViewPoint®{-3.382,0.096,0.076}の値を変えているのですが、一番綺麗に作図するには、どの値にすれば良いでしょうか？個人の好みもあると思うのですが、、

<<"Graphics`ParametricPlot3D`"
SphericalPlot3D[1+LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint®{-3.382,0.096,0.076},ViewVertical®{0,1,0},ImageSize®300]
SphericalPlot3D[1+LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint®{-1.8,0.5,0.1},ViewVertical®{0,1,0},ImageSize®300]
SphericalPlot3D[1+LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint®{-0.8,0.5,0.1},ViewVertical®{0,1,0},ImageSize®300]
SphericalPlot3D[1+LegendreP[2,Cos[thet]],{thet,0,Pi,Pi/30},{phi,0,2*Pi,Pi/15},ViewPoint®{-1.0,0.5,0.1},ViewVertical®{0,1,0},ImageSize®300]

>>175
すいません。文字化けしております。
元々のソースは「153」です。ViewPointの値だけを変えております。

>>174
「ググれ」の意味わかりました。以後「ググり」ます。

＞＞175
事情により追求を残念しましたので回答は結構です。すみません、お騒がせ致しました。

このスレの書き込む人って頭よさそうな書き込みするね

すみません！質問なんですが、
�@２変数関数ｚ＝（2x^2+y^2)exp(1-X^2-Y^2)のグラフを、-2≦X≦2,-2≦Y≦2の範囲でPlot3Dを使って描きなさい。

�A上記関数において∂z/∂x=0,∂z/∂y=0を連立して解いて全ての臨海点を求めなさい。
このとき、∂z/∂x=0,∂z/∂y=0は、どちらもFactorを用いてexp(1-X^2-Y^2)を括りだした後に
exp(1-X^2-Y^2)で割っておき、多項式の方程式にしてからSolveを使うのがよい。

�B各臨界点におけるｚの値を求め、またＨ＝、(∂^2)z/(∂^2)x・(∂^2)z/(∂^2)y　及び、(∂^2)z/(∂^2)x
の符号を調べて臨界点を分類しなさい。

この３つの問題が全然わからなく困っています。どなたか教えていただけるとありがたいです（＞＜）

本当に全然分からないなら君にはできないのであきらめる。
そうじゃないならどこまで考えてどこまで分かってるのか書く。

質問なのですが、グラフの数字のフォントを大きくしたいのですが、
どうすればよいのでしょうか？
わかる方がいましたらどうが教えてください。
（ちなみに現在は一度Illustratorで編集を行い数字のサイズの調節を行っています）

>>182
日本公式サイトのFAQにそのものずばりの項目がある.

ttp://www.cybernet.co.jp/matlab/support/techkwdb/search_view.php?id=1008

>>183
お前。。。

>>181
5ならPlotのヘルプに書いてあるような気がするが。
Mathematica Plot font sizeでググったらいくらでも出てくるし。
何をやるにしてもこれくらいのことは自分で調べる能力をまず身に付けろ。

俺はリンクを間違えない能力をまず身に付けるわ。。

上のは>>182の間違い。

>>183　回答ありがとうございます。
しかし、MathematicaではなくMATLABの説明のページでした。

>>185　申し訳ありません。おもに書籍を用いてMathematicaの勉強をしていたためヘルプの使い方が未熟でした。
もし、よろしければ
Plot[sin[x],{x,0,2}]
のグラフの数字のフォントを大きくした例を示していただけると非常に助かります。

ヘルプに書いてあるって言ってるのになんで見ないの?

>>188 ヘルプを見たところ「FontSize」を使えばなんとかなりそうなのですが、
Plot[sin[x],{x,0,2},FontSize->20]
と書いても、
StyleFont[Plot[sin[x],{x,0,2}],FontSize->20]
と書いてもうまくいかず、どのように記述すればよいのかわからないからです。
私のヘルプ活用能力、Mathematicaのスキルが劣っているためだと認識しております。
どうがお力を貸してください。

うちのヘルプにはそんな風には書いてないんだけど。
ヘルプに書いてある通りに本当にやってんのか?
で上でググれとも書いたんだけどググったのか?

>>190
ご迷惑をおかけしました。解決しました。
ヘルプについてですが使用しているMathematicaのバージョンが4.0のため若干異なっているのかもしれません。（検索ワードに問題があったのかもしれませんが・・・）
「ググれ」についてですが、私は「日本語のページを検索」にチェックを入れて検索を行う癖がありまして、それが原因でうまくいっていませんでした。
親切な返答に対しいい加減な調べ方しかしていなかった自分を恥ずかしいく思います。
本当にありがとうございました。助かりました。

旧バージョン使ってるなんて重要なことは最初に書く。
エラーが出たならエラーメッセージの内容を書く。
解決したなら何が問題だったのか、どうやって解決したのかを書く。

もういないんだろうけど。

>>192
今後は気をつけます。
エラーの原因は>>189で示したような間違った文法を用いていたからです。
どうやって解決したかは、グーグルを用いてMathematica Plot font sizeで検索を行い、
ttp://documents.wolfram.com/mathematica/Built-inFunctions/GraphicsAndSound/2DPlots/FurtherExamples/Plot.html
の文法の通りに修正を行いました。
掲示板を使用する上での私のマナー、常識が欠けていたために>>192の人の気分を害してしまったことを申し訳なく思っております。
>>192の人には掲示板を使用する上での常識、マナーまで教えていただき本当に感謝しております。
ありがとうございました。

そんなことで気分を害したりはしないけども。
掲示板を使用する上で、とかは全然関係ないと思うよ。
問題を解決する基本的なスキルの話だよ。何をするときでも一緒。

こんにちは、

青色の字で、アルファベット「abc」と印字したいのですが、
どのようなコマンドを使用すれば良いでしょうか？

Print[“abc”];

一番直感的な方法は
　　Print["abc"]
と入力した後、abc の部分を選択して、メニューの書式から文字の色を選んで青を選択する。
これじゃかっこ悪いと思うなら
　　Print[StyleForm["abc",FontColor->RGBColor[0,0,1]]]
でやったらいかが．

Mathematica 6.0使っている人に聞きたいんですけど、
5.2のヘルプのようなMathematica Bookがないって本当ですか？

Mathematica 6.0をLinuxで使ってます。
デフォルトのフォントサイズが小さいので、毎回大きく設定しなおして
使ってます。デフォルトのフォントサイズを変更する方法って
ないですか？

とりあえず、Formatから、Screen Environmentを選択して、
そのあと、Presentationを選択することにしました。

>>197
6.0のヘルプは確かに見づらい。検索は向上してるけど、
5.2にあるMathematica Bookのように一から読んでいけるようなものはない。
初めて使う人には使いづらいかも。

ヘルプに頭から読んでいくような内容なんて不要。
そういうもの「も」ついててもいいとは思うが。
基本は強力な検索とハイパーテキストと豊富な実例がヘルプの命。

6のヘルプの嫌なところは、built-in FunctionがAlphabetical Listingでしかlisting出来ないところだな
あと最初のページが見づらい

簡単な式変形で質問です。

z = y Sqrt[m] Sqrt[x^2+y^2]
とやると
z = y Sqrt[m] Sqrt[x^2+y^2]
になると思いますが

出てきた結果を根号で再度まとめることはできるんでしょうか?
具体的には
z = Sqrt[m (x^2 + y^2) y^2]
と変形して欲しいのです。

どうもやり方が分からず困ってます。
どなたか御存じでしたら教えて下さい。お願いします。

mixiと2chへのマルチポストって初めて見たｗ

>>203
どういう場面でそういう操作が必要になるの？

>>205
レス遅れてすみません。えーと、問題になってるのは、 複雑な式でSimplifyできないので困ってます。

Log[Sqrt[A] Sqrt[B]]
(A、Bは非常に複雑な式)
のような場合に上手くSimplifyしてくれないので困ってたんですが、これは
SqrtSimplify[sqrt_] := (sqrt//. Sqrt[x_] Sqrt[y_] -> Sqrt[x y])//Simplify
のような形でSimplifyをかけると簡略化に成功しました。
ですが、今度は、
(-2 m + r0) r[t]
Out[1]= (Sqrt[----------------]
m (r0 - r[t])
Sqrt[r[t]]
> ((Sqrt[-2 m + r0] (4 m + r0) ArcTan[---------------] -
Sqrt[r0 - r[t]]
3/2 Sqrt[-2 m + r0] Sqrt[r[t]]
> 4 Sqrt[2] m ArcTanh[-------------------------------])
Sqrt[2] Sqrt[m] Sqrt[r0 - r[t]]
> Sqrt[r0 - r[t]] + Sqrt[-2 m + r0] Sqrt[r[t]] (-r0 + r[t]))) /
> (Sqrt[2] Sqrt[-2 m + r0] Sqrt[r[t]])

のような複雑な式でどうにもFullSimplifyがききません。本当は分母がSqrt[2m - r0]とSqrt[2] m で簡略化されるはずなんですが・・・
ここで上手く分母を上記でまとめないと、r[t] -> r0 などと仮定して計算を行うので1/0の形でエラーで止まります。
どうも良く分かりません。。

↑式が見にくかったので再度書きます。


　　　　　　　　　(-2 m + r0) r[t]
Out[1]= (Sqrt[----------------]
　　　　　　　　　　m (r0 - r[t])
　　　　　　　　　　　　　　　　　　　　　　　　Sqrt[r[t]]
> ((Sqrt[-2 m + r0] (4 m + r0) ArcTan[---------------] -
　　　　　　　　　　　　　　　　　　　　　　　Sqrt[r0 - r[t]]
　　　　　　　　　　　　　　3/2 Sqrt[-2 m + r0] Sqrt[r[t]]
> 4 Sqrt[2] m ArcTanh[-------------------------------])
　　　　　　　　　　　　　　Sqrt[2] Sqrt[m] Sqrt[r0 - r[t]]
> Sqrt[r0 - r[t]] + Sqrt[-2 m + r0] Sqrt[r[t]] (-r0 + r[t]))) /
> (Sqrt[2] Sqrt[-2 m + r0] Sqrt[r[t]])

googleで検索できないので質問させてください。
数式 /. {条件}
のときの/.ってどういういみですか？

頭使え。グーグルなら Mathematica slash period で検索。
Mathematica持ってるなら?/.でヘルプが出る。

>>203、207
203の式はy>0じゃないと成り立たないよね
変数に条件与えてRefineでなんとかならんかね
少しは手で式を変形しとかんと上手くいかんだろうけど

こんにちは、
下記のような数値計算で、yが０に最も近づく
ｎを求めたいのですが、どのようなプログラム
になるのでしょうか？
但し、ｙは、０以上であり、
yが０に最も近づいたときの、nとｙのみを
知りたいです。
下記は、あくまでも例であり
N[4000/95]
により、
n=43
y=4085
なので、この例の数値を知りたい訳ではなく、
プログラムの仕方を知りたいです。



For[n=1,n£100,n++,
y=4000-95*n;
If[y==0,Print[n]];

];

>>211
どれくらいの精度でｎを求めたいかによって変わるんじゃない。
与えられた情報だけでは答えるのが難しい。

お返事ありがとうございます。
質問の仕方が悪かったです。

一応下記で出来ました。
すいませんでした。
For[n=1,n<100,n++,
y=4000-95*n;
If[y<0,Print[n]];
If[y<0,Print[y]];
If[y<0,Break[]];
];

こんにちは
Solve[4-x^2Š0,x]

を計算しますと
X=-2,
X=2

の解が得られますが、
X=2
だけを取り出して、
y=x^3+1
でしたら、
９
とするには、

{{x->-2},{x->2}}

から、どのように２だけを
取り出せばよいでしょうか？

x=.;y=.;
y=x^3+1/.Solve[4-x^2==0,x][[2]]

y=.;
は不要でした。

数式計算でCとかC++とかいってるやつってあふぉなのか？
普通フォートランだろ

>215

お返事ありがとうございます。

この方法もありますね。
下記で、x->は不要です。
２だけを取り出す仕方はどうでしたか？

Flatten[Solve[4-x^2==0,x]][[2]]

>>218
質問の意味が？？？

>>218
質問の真意はよく解らないが、単純に数値の2のみを取出したいのなら
　　x=.
　　x/.Flatten[Solve[4-x^2==0,x]][[2]]
でよいと思う。但しこの場合xには2が代入されていない。これを実行した後に
　　x
を実行すると
　　x
が返ってくるだけである。

もしxに2を代入した状態にしたいというのであれば
　　x=.
　　x=x/.Flatten[Solve[4-x^2==0,x]][[2]]
としなければならないように思う。これを実行した後に
　　x
を実行すると
　　2
が返ってくる。

>>220

お返事ありがとうございます。
よくわかりました。
x=x/.Flatten[Solve[4-x^2==0,x]][[2]]
で、希望のものを得られました。

ParametricPlot3D[{(2 + 0.5Cos[u])Cos[v], (2 + 0.5Cos[u])Sin[v], 0.5Sin[u]}, {u, 0, 2Pi}, {v, 0, 2Pi}]
を使用し３Dの図形を書いたのですが、それをどうすれば回転させることができますでしょうか？よろしくお願いします

Print["足し算テストです."]
i = 0; atari = 0; hazure = 0;
While[i < 10,
{com1 = Random[Integer, {0, 100}];
com2 = Random[Integer, {0, 100}];
com3 = com1 + com2;
Print[com1, "+", com2, "=?"];
you = Input[];
If[you == com3,
{Print["当たり"]; atari = atari + 1},
{Print["はずれ"]; hazure = hazure + 1}];
i = i + 1
}]
Print["ゲ-ム終了です."]
kekka = atari/(atari + hazure)*100;
Print["あなたの点数は", kekka, "点です."]
If[kekka > 80, Print["合格"],
Print["不可"]]

上記のような簡単な計算のゲームを、今度は足し算の問題と引き算の問題がランダムに入り交じって１０問出題されるようにするには、どうすればよいのでしょうか？
IfやRandomを使っても、僕の力ではよくわかりませんでした。
どなたかよろしくお願い致します。

ｻｲｸﾛｲﾄﾞの作り方教えてください

ParametricPlot[ {t+Sin[t], Cos[t]}, {t, 0, 2 Pi} ]
サイクロイドだよ

すみませんがこの問題誰か解けませんか？

９行９列の整数型配列変数uzu[ ] [ ]に1から1づつカウントアップ下値を渦巻き状に代入するプログラムを作る。整数型の配列変数はuzu[9 ] [9 ]とする。

こんな感じにしたいんですよ。

1 2 3 4 5 6 7 8 9
32 33 34 35 36 37 38 39 10
31 56 57 58 59 60 61 40 11
30 55 72 73 74 75 62 41 12
29 54 71 80 81 76 63 42 13
28 53 70 79 78 77 64 43 14
27 52 69 68 67 66 65 44 15
26 51 50 49 48 47 46 45 16
25 24 23 22 21 20 19 18 17

>>226
感覚的には９×９のマスを動くロボットのようなプログラムを作ればいい。
壁か自分が既に通ったマスにぶつかったら自分の進行方向から見て右に曲がる。
次のマスにつくたびにカウントを１プラスしてuzu[x] [y]にその値を書き込めばそんな感じになる。
はず。

学校で宿題がでたのですが、さっぱりわかりません。
期限が明日の１８：００までなので急いでいます。
以下問題文

ランダムウォーク

(1)
平方和が１になるようなｘ成分、ｙ成分を持つ乱数を発生させる関数randomXY[]を作成せよ
ヒント：単位円上に一様に値を持つ乱数を考えよ．
つまり，qを2 Piの乱数として，2つの要素x=Cos[q], y=Sin[q]を持つリストを生成する
乱数の発生については，組み込み関数のRandom[] のHelpを参照

(2)
(1)で作った乱数の1000回分のリストを2次元平面上にプロットし，
ランダムな動きを実現していることを確かめよ
（これがランダムウォークの動く方向と大きさを表示している）
ヒント：randomXY[]を1000回適用して得られるTableを作成せよ
（そのリスト・テーブルの名前をrandomWalkStepとする）

(3)
randomWalkStepをListPlotせよ．
FoldList，Plusなどの関数を用いてrandomWalkStep 値を順に加えた、
randomWalkというリストを生成せよ
randomWalkを2次元の図に表せ（これがランダムウォークの粒子の動きの軌跡となる）
HelpのListPlotの項目を参考に

長レスすみません。
お願いします。

先生にここ見られたらおしまいじゃん。
まあもう期限過ぎたしいいか。
randomXY[] := (q = 2Pi*Random[]; {Cos[q], Sin[q]})
randomWalkStep = Table[randomXY[], {i, 1, 1000}];
ListPlot[randomWalkStep]
randomWalk = Transpose[{FoldList[Plus, 0, Transpose[randomWalkStep][[1]]],
　　　　　FoldList[Plus, 0, Transpose[randomWalkStep][[2]]]}];
ListPlot[randomWalk]
で出来るはず。

宿題の答えを聞くバカ
教えるバカ

Windowsの[Alt+Tab]のように、
アクティブなノートブックを変えるショートカットキーってあります？

スレ違いかも知れないけど。

linuxにMathematica6.0がインストールできません。

MatheInstallerを実行すると、最後に

A valid password file entry already exists.

Installation complete.

となって成功したようにみえます。でも、起動しようとすると

# mathematica　&
Mathematica front end executable
/usr/local/Wolfram/Mathematica/6.0/SystemFiles/FrontEnd/Binaries/Linux/Mathematica
not found. Your Mathematica installation may be incomplete
or corrupted.

　OS は fedora7 です。

どなたかわかりますか？？

式の表示につきまして

こんにちは、

下記の式、ｙ１とｙ２があります。
ｙ１を評価させて、ｙ２の通り表示させるには、どうしたら良いでしょうか？



y1 = Expand[(a1*p1 + a2*p2 + a3*p3 + b*mc)^2]


y2 = p1^2*a1^2 + p2^2*a2^2 + p3^2*a3^2 + mc^2*b^2 + (a1*a2 + a2*a1)*p1*p2 + (a2*a3 + a3*a2)*p2*p3 + (a3*a1 + a1*a3)*p3*p1 + (a1*b + b*a1)*p1*mc + (a2*b + b*a2)*p2*mc + (a3*b + b*a3)*p3*mc



ちなみに、ｙ１をそのまま評価しますと、下記になり、ｙ２の通り表示しません。

p1^2*a1^2 + 2*p1*p2*a1*a2 + p2^2*a2^2 + 2*p1*p3*a1*a3 + 2*p2*p3*a2*a3 + p3^2*a3^2 + 2*mc*p1*a1*b + 2*mc*p2*a2*b + 2*mc*p3*a3*b + mc^2*b^2

どなたか教えてください。
mathematicaで非集計ロジットモデルのパラメータ推定できるのでしょうか。
Rでは出来るようですけど、mathematicaの方が慣れてるのでできればmathematicaでやりたいのです･･･

# export LANG=C
# mathematica &
とやってみそ。

student　ver持ってるんですが、日にちが過ぎたので
ウインドウズの右下の日付を少し戻したんですが使えないままです。
ウインドウズの日付というのは値渡しということなんでしょうか？

Mathematicaで
For[i = 1, i <= 500, i++, xdate1[i] = Random[Real, {-1, 1}]]
For[i = 1, i <= 500, i++, xdate2[i] = Random[Real, {-1, 1}]]
For[i = 1, i <= 500, i++, If[(xdate1[i])^2 + (xdate2[i])^2 > (2/Pi), ydate[i] = 1, ydate[i] = -1]]
ってとこまでやったんだけど
ListPlot[Table[{xdate1[i], xdate2[i]}, {i, 1, 500}]]]]
のグラフ書いたとき
ydate[i]=1の部分だけ色変えたいときってどうすればいい？

>>237
グラフのオプションはバージョン5以前と6以降で大きく違うから、
バージョンも書いたほうがいいよ

楽なのはydate[i]=１のときのグラフと-1のときのグラフを別々に作って、
Showで表示する方法だと思う

あとdateじゃなくてdataじゃないの？

>>237
For[i = 1, i <= 500, i++, xdate1[i] = Random[Real, {-1, 1}]]
For[i = 1, i <= 500, i++, xdate2[i] = Random[Real, {-1, 1}]]
tbl1 = {}; tbl2 = {};
For[i = 1, i <= 500, i++,
If[test[i] = (xdate1[i])^2 + (xdate2[i])^2 > (2/Pi),
AppendTo[tbl1, {xdate1[i], xdate2[i]}],
AppendTo[tbl2, {xdate1[i], xdate2[i]}]]]
ListPlot[{tbl1, tbl2}, PlotStyle -> {Blue, Red}]

こんな感じかな
ver6だけど、5以前でも変わらないと思う

おお、サンクス

While を使って1 から整数i までの間の素数の数を数えるプログラムを作
れ。整数i は最初の行にi=4; のように与えること。

これ教えてくださいお願いします

maximaと文法似てるの？

>>241
PrimeQ使えば至極簡単なんだけど、使ってもいいの？

i=4;n=2;c=0;While[n<=i,n=NextPrime[n];c++]

初歩的なことだと思いますが，
定義域ごとに関数を定義したいんですが
例えば，
For[i=0,i<100,i++,x[t_] := 0.5^i /; i<t<i+1]
とするとうまくいきません．
誰か教えてください．

>>245
例えが何がしたいのか良くわからないが、
定義域が引数に対してなら、こんな感じになる

x[t_ /; t >= 0] := Log[t]
x[t_ /; t < 0] := Log[-t]

関数xは引数tが0以上のときにLog[t]、0以下のときにLog[-t]を返す

質問の意図がよくわからんかったから見当はずれの答えかもしれん

説明不足ですいません。
x[t_ /; 0<t<1] := 0
x[t_ /; 1<t<2] :=1
とすれば、階段関数みたいのものを定義できますが、
これを繰り返し、定義域 i<t<i+1に対して
例えば 関数
x[t_] := i
みたいなもの定義したいのです。
もちろんこのケースは
x[t_] := Integer[t]
とすれば階段関数が定義できますが、
もっと複雑に i に依存して、
それぞれの区間 i<t<i+1 にたいして
Forループで関数を定義したいのですが、
For[i=0,i<100,i++,x[t_] := i /; i<t<i+1]
ではうまく定義されません。
どうしたらよいでしょうか？

>>247
反復演算子を使ってるのが良くないみたいですね
局所定数にしてしまえば、上手くいくんじゃないでしょうか

For[i = 0, i < 100, i++, With[{j = i}, x[t_] := j /; j < t < j + 1]]

>>246
大変ありがとうございました。
うまく行きました。
たすかりました。

mathematicaで微分演算子を含むハミルトニアンを独立に扱いたのですが，できるのでしょうか？

例えば，一次元調和振動子のハミルトニアンでしたら

H＝-d^2/dx^2+x^2

ですが，これを独立に扱って

(H-E)^2*ψ

などの計算をmathematicaで簡便に行いたいのですが（ここでψは簡単な指数関数，ψ＝exp(-a*x^2) など），このような計算のプログラミングはどのようにすればできるのでしょうか？

こんにちは、
下記ＨＰをコピーして、実行させたいのですが、そのままコピーして使おうとしますと文字しかコピーできません。手で文字を１字ずつ入力して実行させる以外に方法はないでしょうか？

http://www-antenna.ee.titech.ac.jp/~hira/hobby/edu/commu/ss/index.html

>>250
（・∀・）つ 「物理学のためのMathematica―古典力学から宇宙論まで」
　　　　　　　ロバート・L. ジンマーマン (著), フレデリック・I. オルネス (著)

マスマティカでwigner分布できるパッケージある？

マスマティカってwigner分布啓さんできる？

interpolationを使ったら
座標が構造化テンソル積グリッド上にありません
って言われた。

ヘルプにも乗ってないんだけど、ヘルプが古いのかな？

MathematicaにDistributeという関数があるのですが、この関数を使って自分で作った関数に分配則を適用できます。
ここで問題なのですが、
Distribute[f[a+a+b,c]]と入力をすると、f[a,c]+f[a,c]+f[b,c]と出力してほしいのに、
f[2a,c]+f[b,c]と出力されてしまいます。
f[a,c]+f[a,c]+f[b,c]と出力するには、
a+a+bの部分を勝手に計算してしまうMathematicaの悪い癖を正さないといけないと思うのですが、
何か良い方法はないでしょうか？

すいません。初歩的な質問なんですが
n乗の計算式を求めたいんですが・・
mathematicaでは自分で数値を入力するようなことはできないのですか？

>>256
n乗の計算式って何？
x^nと記述すれば、xのn乗が求まります。
具体的な数値を知りたければ
In:N[Pi^3]
Out:31.006277
となると思います。
知りたい事がよく分からないので、勝手に回答しちゃいました。

>>255
解決しました。
SetAttributes[f, HoldFirst]で属性を加えると
第一引数については計算しないで固定しておいてくれるみたいです。
しかし残念ながら、第二引数、第三引数については固定しておく方法が見つかりませんでした。
まぁ、問題は解決したので良かったのですが・・・。

Mathematicaは数値計算も出来るけど
本格的な数値計算ならMATLAB系がいいんだよね？

?-g1.
input=taro.
input=jiro.
input=2.
input=..
.
input=end.
yes
endが入力されるまで繰り返して、入力をするg1を作りたいんですけど、全くわかりません。知恵を貸して下さい。お願いします。

アルファベット(aからz)までのさいころをつくって1000回投げて、目の頻度をカウントするにはどうしたらいいんですか？

>>261

RandomChoice[CharacterRange["a", "z"], 1000]
Tally[%]

自動焼人 ★　＝　自動保守 ◆KAWORUKOFI　＝　自動保守#K9K?_D[L

名言集　その２
『お前が規制系キャップ取れるか審査してやるよ』

http://yutori7.2ch.net/test/read.cgi/news4vip/1249830540/　ID:PVAf+dux0　＝　自動焼人 ★

> 36 ：以下、名無しにかわりましてVIPがお送りします [sage] ：2009/08/10(月) 00:31:30.02 ID:PVAf+dux0
> >>33
> キャップとコテハンの違いは何？

> 46 ：以下、名無しにかわりましてVIPがお送りします [sage] ：2009/08/10(月) 00:38:05.34 ID:PVAf+dux0
> >>45
> その回答では落ちるなｗ
> 答えは教えないがｗ

> 50 ：以下、名無しにかわりましてVIPがお送りします [sage] ：2009/08/10(月) 00:41:29.96 ID:PVAf+dux0
> Q.キャップとコテハンの違いは何？
> A.２ちゃんねるのボランティアの登録制度

> それがお前の答えかｗ

> 52 ：以下、名無しにかわりましてVIPがお送りします [sage] ：2009/08/10(月) 00:43:10.06 ID:PVAf+dux0
> まぁ、どうせ正解が出るわけもないし、次の問題。
> 君が思う面白いスレはどんなの？
----------------------------------------------
この自動焼人 ★メールマガジンの配信停止をご希望される方は
http://qb5.2ch.net/test/read.cgi/sec2chd/1250169591/
にて自動焼人 ★までご連絡ください

BubbleChart[{{0.5,0.5,0.5}},PlotRange->{{0,1},{0,2}}]
とやると、バブルが楕円になってしまうのですが、
PlotRangeを使いつつ、バブルを円にするにはどうしたらいいですか？

ダミーデータ入れて、自己解決しました。
BubbleChart[{{0.5,0.5,0.5},{0,0,0},{1,2,0}},PlotRange->{{0,1},{0,2}}]

2変数の偏微分方程式(境界はすべてディリクレ型)を解こうとしています。
変数の範囲はx,yともに0〜1で、d^2u[x,y]/dx^2+d^2u[x,y]/dy^2=0
が支配方程式です。境界条件として、
u[x,0]=sin(2*PI*x), u[x,1]=-sin(2*PI*X),
u[0,y]=0,u[1,y]=0
としています。これをNDSolveで解かせようと思い、以下のコード
sol1 = NDSolve[ {
D[u[x, y], x, x] + D[u[x, y], y, y] == 0 ,
u[0, y] == 0,
u[1, y] == 0,
u[x, 0] == Sin[2 Pi x],
u[x, 1] == -Sin[2 Pi x]
},
u[x,y],
{x, 0, 1},
{y, 0, 1}
]
を実行しても
Boundary values may only be specified for one independent variable. \
Initial values may only be specified at one value of the other \
independent variable.
と言われて結果が出ません。境界条件の与え方が悪いのでしょうか？
どうしたらいいか教えてください。

すいません、教えてください。
たとえば、
P[l_, x_]:=D[(x^2-1)^l,{x,l}]
と定義して
P[1,t]とすると2tと正しい結果になるのですが
P[1,1]とするとエラーになります。
t=1のときの値を計算するにはどうすれば
よいのでしょうか。

>>267
>P[l_, x_]:=D[(x^2-1)^l,{x,l}]

純関数を使えばできます。

P[l_, x_] := Derivative[l][(#^2 - 1)^l &][x]

ありがとうございます！
おかげでルジャンドル関数が定義できました。
（mathematica calccenterを使ってるのですが
　特殊関数が入っていないのです・・・）

大学で
一様な電場中に導体球がおかれたときの電位・電場の様子を可視化せよ
という課題が出たんですが
特に境界条件とかどうしていいのかよくわかりません

誰かヒントや解いてくれる方いないでしょうか？
提出期限は、10月７日の18:00までです。お願いします。

マルチw

どなたか下記のプログラミング（mathematicaです）の誤りを
ご指摘頂けないでしょうか？
Labelが見つからないと表示されます。
教本２冊と奮闘しておりますが、わかりません。

x = -10;
d = 2;
Label[begin]
s = x^2 - (x + d)^2
If[d/2 < 0.01, Goto[end], Goto[hazime]]
Label[hazime]
If[x^2 - (x + d)^2 >= 0, {x = x + d; Goto[begin]}, {d = d/2;
Goto[begin]}]
Label[end]
Print[x]

宜しくお願いします。

なんでgotoなんてつかってるの？

>>272
やりたいことがよくわからないのだが、こう書いたほうがいいのでは?

x = -10;
d = 2;
While[d/2 >= 0.01,
If[x^2 - (x + d)^2 >= 0, x = x + d, d = d/2];
]
Print[x]

Goto は使わないほうがいい。

こんにちは、
下記HPの図１は、直列共振回路です。その波形は、図２です。更に、mathematicaで計算したのが、図５です。（下記にプログラムを示します。）
では、図７のようなインパルス列を入力するようなプログラムはどのように書けばよいでしょうか？波形は図８のようになるはずです。
http://www.geocities.jp/dirac_equation/index.htm

Clear[m, l, k, F0, w, x, t]
solution = DSolve[{m*Derivative[2][x][t] + l*Derivative[1][x][t] + k*x[t] == F0*Cos[w*t], x[0] == 0, Derivative[1][x][0] == 0}, x[t], t];
solutionNew = FullSimplify[solution];
x[t_] = x[t] /. solutionNew[[1]];
L = 300/10^6;
R = 10;
c = 47.59/10^12;
V = 1/10^3;
Print["(*振動数　f　が　固有振動数の場合*)"];
f = 1/(2*Pi*Sqrt[c*L])
f = 1.332*10^6;
w = 2*Pi*f;
Plot[Evaluate[x[t] /. {m -> L, l -> R, k -> 1/c, F0 -> V*w}], {t, 0, 200/10^6}]

インパルス列の条件です。
Pulse Amplitude 5V
Period 3us
Pulse Width 300nS
Rise Time 1nS
Fall Time 1nS

u(x)=x^(2/5)
の効用関数の(a,u(a))における接線をu(x)のグラフに書き加えるためのModuleを作成したいのですが
どうModuleを使ったらいいのかよくわかりません・・・

別々に計算してグラフに描くだけならできるのですが…(ヽ'A`)

>>276

Module を使うというところがよくわからないのですが、単に

u[x_] := x^(2/5)
a = 1;
Plot[{u[x], u'[a] (x - a) + u[a]}, {x, 0, 3}]

では駄目なのでしょうか?

>>277

u[x_] := x^(2/5)
Plot[u[x], {x, 0, 5}, AxesLabel -> {"x", "u(x)"}]

上記でかいたグラフに点(a,u(a))における接線を書き加えるModuleを作りたいんです
tangent[a_] := Module[u'[x]*(x - a) + u[a]]

としてtangent[1.5]などと入れるとそのaの値における接線が書き加えられたグラフを描くようにしたいのですが
うまくいきません

>>278

Module は使わなかったのですが、

u[x_] := x^(2/5)
g1 = Plot[u[x], {x, 0, 5}, AxesLabel -> {"x", "u(x)"}]
tangent[a_] := Show[g1, Plot[u'[a]*(x - a) + u[a], {x, 0, 5}]]

でやりたいことが実現されるでしょうか?

>>279

u[x_] := x^(2/5)
tangent[a_] := Module[a],
Plot[u[x], u'[a]*(x - a) + u[a], {x, 0, 5}, AxesLabel -> {"x", "u(x)"}]

こんな感じにしたいのですができませんでした・・・

>>280

こんな感じでいけます。

u[x_] := x^(2/5)
tangent[a_] :=
Plot[{u[x], u'[a] (x - a) + u[a]}, {x, 0, 5}, AxesLabel -> {"x", "u[x]"}]
tangent[1]

ご教授下さい！
4次元データのcontourを書きたいのです．
ControurPlot3Dを使えば良さそうなのは分かるのですが，使い方がわかりません．
陰関数f(x,y,z)=0の球面上にg(x,y,z)の大きさでcontourを描きたいんです．

http://tech.ckme.co.jp/img/sphere8cdata_matlab.png
こんな感じです

よろしくお願いします！

現在単一パーコレーション・クラスター・プログラムを作成しています。
とりあえず以下のプログラムになりました。

In[1] := Epidemic[n_, p_] :=
Module[{choices = {{1, 0}, {0, 1}, {-1, 0}, {0, -1}},
reject, pickAndChoose, select, newpers},
reject = {};
pickAndChoose :=
(select = #[[2, Random[Integer, {1, Length[#[[2]]]}]]];
If[Random[] <= p,
newPers =
Complement[Union[Map[Function[y, y + select],
choices], #[[2]]],
{select}, #[[1]], reject];
{Join[#[[1]], {select}], newPers},
reject = Join[reject, {select}];
{#[[1]], Comple, ent[#[[2]], {select}]}]) &;
FixedPoint[pickAndChoose, {{{0, 0}}, choices}, n,
SameTest -> (#2[[2]] == {} ||
Length[#2[[1]]] == n &)][[1]]
]

上を実行すると

SetDelayed::write : In[1]のタグInはProtectedです．

と、エラーが出るんですけど、どうしたら解決しますか？
初心者でヘルプを見てもよくわからなくて・・・わかる方いましたら教えてください。

maximaみたいにlispで拡張作ることできますか？

>>542
なるほど。ありがとう

１００ステップの2次元ランダムウオークを作れ。
ただし、東西南北をそれぞれ、｛１、０｝｛ー１、０｝｛０、ー１｝｛０、１｝で
表すことにする。
（ヒント：りすとdirectotions={{1,0},{-1,0},{0,-1},{0,1}}を作り、1、2、3、4の範囲の
乱数kを発生させ、directotionsのk番目の要素をとるとランダムな方角が得られる）

すいませんが練習だと思って解いてくれませんか？課題が出せないんでお願い

Graphics[Line[Accumulate[RandomChoice[{{-1, 0}, {1, 0}, {0, 1}, {0, -1}}, 100]]]]

整数nの各位の数字の合計を出力する関数f4[n]を書け。
例えば、f4[1234]=10.

わかる人います？教えてくれませんか？

>>288
ヒント：1234 をなんとかして {1,2,3,4} の形にして足す。

教えて

正整数nのケタ数を与える関数f2を書け。たとえば、f2[40172]＝５

↑
問題間違えました。

０より大きく１より小さい既約分数で分母がn以下のものを
小さい方から並べてできるリストを出力する関数f5[n]を書け。たとえば、

f5[n]={1/4,1/3,1/2,1/3,1/4}

（このような数列をnに対するFarey数列といい、いろいろおもしろい性質が知られている。
たとえば、隣り合う項をa/b,c/dとすると、必ずab-cd=-1となっている。）

すいませんがお願いします。

>>291

あまり美しくないが。

f5[n_] := DeleteDuplicates[Sort[Flatten[Table[Table[Which[0 <= a/b <= 1, a/b, True, {}], {a, 0, n}], {b, 1, n}]]]]

ありがとうございます

正整数nの各位の数字を逆転させた数をnの逆転数といい、r(n)と書く。
たとえば、f[1234]＝4321

r(n)=n/4となるような最小の正の整数を求めよ。

お願いします。意味わかりません。

まず

f[1234]＝4321

を作るとどうなる？

どうなるって？

まず、1234 を 4321 にするにはどうしたら良い？

わからん

5555-1234

それだと 2345 のときに 5432 にならない。
もっと頭を柔らかく、Mathematica 的に考える。
Mathematica が得意なのは何か？

ヒント：>>288

正の整数リストから元の順序を変えずに重複するものを取り除く関数を書け。
{3,1,1,5,2,1,2}み適応し、{3,1,5,2}になることを確かめよ。
（Unionを使うと{1,2,3,5}になるので注意）

なんか方法ある？

タクシーの番号を見てラマヌジャンが
「その数は2通りの仕方で2つの正の立方数（立方数とはn^3）の和で表される
自然数のうち最小のものです」さてその数は何か？

答えはたしか９１になればいいと思います。
僕にとってムズイので教えてください。

91 = 6＾3 + (-5)＾3 = 4＾3 + 3＾3

タクシーの番号を見てラマヌジャンが
「その数は2通りの仕方で2つの正の立方数（立方数とはn^3）の和で表される
自然数のうち最小のものです」さてその数は何か？

答えはたしか1729になればいいと思います。
僕にとってムズイので教えてください。

6より大きな偶数nは二つの素数p,qの和で書けると予想できている。（ゴールドバッハ予想）
6より大きな偶数nに対し、n=p+q(p<=q)となるような素数の｛p,q｝を
与える関数goldbach[n]を書け。
goldbach[1000]を実行せよ。

goldbach[n_] := Table[{i, n - i}, {i, 3, n/2}]
までやったんだけど、goldbach[10]をやると、
{{3, 7}, {4, 6}, {5, 5}}になるんだ。どうすりゃ｛4，6｝消える？
PrimeQはどこで入れるの？最初から違う

正整数nが対称数とは、nの逆転数とnが等しいものである。
例えば１２３２１は対称数である。
正整数nに対し、それが対称数ならnを出力し、
そうでなければ、nとその逆転数を加え、それが対称数かどうか調べる。
この操作を対称数が得られるまで繰り返し、
得られた対称数を出力するような関数f[n]を書け。
たとえば、
78→78+87=165→165+561=726→726+627=1353→1353+3531=4884
なのでf[78]=4884である。

わかんないです。助けてください。

マセマティカによりモンテカルロ法を用いて、何か平面図形の面積を求めてください。
そして、その数値と求積の公式や積分を利用して得られる理論値との違いを調べてください。
すなわち答えることは以下の4点です。
（１）どうゆう図形について考えましたか？
（２）その図形の面積をモンテカルロ法で求めるとどうゆう結果がでましたか？
（３）その図形の面積を数学で求めるとどうなりますか？
（４）あなたがモンテカルロ法で計算した結果は理論値に対して何％のずれがありますか？

ポーカー
トランプのカードを
{H,1},{H,2}.....,{H,13}（ハート）
{S,1},{S,2}.....,{S,13}（スペード）
{D,1},{D,2}.....,{D,13}（ダイヤ）
{C,1},{C,2}.....,{C,13}（クラブ）
で表す。
（１）52枚のカードからランダムに5枚選び出す関数を書け。
（２）5枚の手札の（ポーカーにおける）点を下の表に従って与える関数を書け。
（３）（１）を1000回繰り返したときの平均点数を求めよ。

（表）
ワンペア　5点　ツーペア　10点　スリーカード　15点
ストレート　20点　フラッシュ　25点　フルハウス　30点
フォーカード　40点　ストレート・フラッシュ　50点
ロイヤル・ストレート・フラッシュ　100点

行列計算のsolveはできないのか
solve(A.x=b,x)

すいません。よくわかりません。

ベクトルと行列の方程式を解く方法はないのかってことだよ

できないかも。
俺へたっぴだから。

例えば、Erf関数名は知ってるけど、数学的な誤差関数の定義を知りたいときに、使えるコマンドってありますか？
Hoge[Erf]ってすると、∫e^t2dtって感じの出力が得られるような関数です。

Hoge[exp]の出力は

>>313
なるほど。ありがとうございます。

なにこれ

mathematicaの
.nbファイルをもらったんだけど，
mathematicaでは簡単な計算をしたことはありますが，
プログラムをしたことがありません．

定義されている変数とか，関数とかを見るにはどうすればいいですか？

うーん

http://reference.wolfram.com/mathematica/ref/Names.html
によると
Names["List*"]
ってやるとListで始まるシンボルが列挙されるらしいぞ。

http://reference.wolfram.com/mathematica/ref/Information.html
によると
シンボルが分かったら
? シンボル
でいいらしいぞ。

オライリーのマスマティカ本欲しい！ほしゅぃいいいいいいい！

こんなのがあるよ。

「Mathematica Cookbook 日本語版」出版記念講演
http://www.hulinks.co.jp/software/mathematica/seminar_cookbook.html

4月27日(水)14:00〜17:00
参加費無料

Rubyバカにしてる子ってさ
変数に＄ついてる言語触ってるって事だよね

いちいちSHIFT＋４キーおして ＄ 打ちまくってる感触はどう？


はい次（笑）

∨∨∨∨∨
＜＜＜＜＜＜＜＜＜＜＜ こんなのがあるよ。 ＞＞＞＞＞＞＞＞＞＞＞(キリッ!キリッッッ！！ッ！！！
∨∨
＜＜＜＜＜＜＜＜＜＜＜＜ 「Mathematica Cookbook 日本語版」出版記念講演 ＞＞＞＞＞＞＞＞＞＞＞＞(キリッッッッ!キリッッッッ！！！！キリッッッッッ！
∨
＜＜＜＜＜＜＜ 4月27日(水)14:00〜17:00 ＞＞＞＞＞＞＞（←キリッ!!!きリッッ！！！キリ!!キリ!!ッッ！！
ゴミじゃねーか

２０１１年、Ruby,Perl,PHP,Pythonって並べたときにさ
ここで、Ruby以外を選ぶ奴ってマジでなんなんだろうな
ゴミだよなぁほんと

こんばんは。質問があります。
ペーペー野郎を助けて頂けませんか？
Mathematicaって、けっこうバグとかあったりするものです？

非線形方程式、非線形常微分方程式とか解かせているんですが、
確認のために、解を元の方程式に代入しても方程式が成り立た
なくなってしまうんです・・・

何か特別なことをしないといけないのですかねぇ・・・。

数学板のほうが活発みたいだよ
〓 Mathematica 5 〓
http://kamome.2ch.net/test/read.cgi/math/1285859504/

http://www.mathworks.co.jp/products/demos/shipping_ja/matlab/_graf2d2_06_thumbnail.png
みたいなプロットできる？

>>327
ListPointPlot3Dが近いかも

ある関数の変数の変化率を連続ではなく、１、２、３・・・・とΔ＝１で飛び飛びにとってプロットしたいのですがどうすればいいのでしょうか？

日本語における文字セットと符号化方式の関係についておしえてくださいい

Newton法を用いて与えられたαに対しての√αの近似値を求める
xn=xn-1- x^2n-1-α/2xn-1=1/2(xn-1+α/xn-1) (n=2,3…)
を使って計算するのって、そのままだと駄目なんですか？
最近マセマティカを勉強し始めたばかりなので、変な質問ですみませんがどなたか教えてください。

なんでこのスレこの板にあるんだ?

数学板はAA貼る人がいるからじゃ？

リストを２次元の散布図にプロットした後、誤ってリストを消してしまったのですが、
作成した散布図のデータからリストのデータを取り出す方法ってありますか？
散布図に対してFullFormすれば一応わかるのですが、散布図の数が多いのでリストデータで欲しいのですが
どなたか回答お願いします。

tabで区切られたデータファイルの読み込み方法を教えてください

状態密度ρ(∈)=Σδ(∈-E(k))
を表すグラフを書きたいのですが
どうしたらいいでしょうか？
∈での箱を自分で定めて、確率がだとかなんとか…

みんなホームエディション買っている？

ドキュメンテーションセンターが「Mathematica 9」になってる。

新バージョン ｸﾙ━━━━(ﾟ∀ﾟ)━━━━!!

ふと思ったこと。

2004 11/18 〜 2013 2/11

この先、何年続くの？

書き込みが1000になるのは2028年くらいかな

あと10年以上あるね。
それまでに Mathematica の使い方マスターしたいと思います。

一日一善とかいうのは削除されたのか？

For文の中の変数iを、回す文章の中のファイル名に入れたいのですが、どう書けばいいのでしょうか？
例えば、Import="フォルダ名//基本名i"という形で、iによって読み込むファイルを変えたり、基本名i=○○○という変数をFor文の中で定義したりしたい場合はどうすればいいのでしょうか。

こんにちわん。

多対一マッチングってmathematicaでどうすりゃええんや

RasberryPi、数式処理ソフト「Mathematica」を同梱。個人利用は無料
http://internet.watch.impress.co.jp/docs/news/20131122_624731.html

「一日一善」の書き込み消えてるやん。
まあ。いいけど。

1.98のような近似値を198/100のような厳密値に変換する組み込み関数はありますか？

あるよ

書きたい答えを349が全く同じに書いてくれた
だからサービス　Rationalize[1.98, 0]　てかマニュアル見たがはやいよ

何向けの言語だこれ

数式処理システムとそれを扱うための言語