2010年8月15日 星期日

[自動轉寄] 引圖疑問

作者: MIZUYAMA (致命祈願)
標題: 引圖疑問
時間: Mon Aug 16 00:41:29 2010

首先感謝您昨晚的幫忙!
才沒有讓我今天下午要給老闆的報告開天窗
畢竟老闆急著要投稿....囧

這是我第一次使用mac打tex系列
之前則都是使用win
在沒遇到您我可能會用bootcamp開win打了 囧

--
昨晚您幫我改成eps離線後
我檢視您幫我轉成的程式碼

一開始以為是轉圖的格式
可是後來實驗套用其他用R語言畫出的圖發現都可以跑出來
相比之下我無法跑出圖的程式碼與您所轉eps的程式碼幾乎沒有差別

想問您方便撥冗說明一下嗎?

--
「如果昨天到了明天會變成兩天前,試問後天應該是今天或是昨天?」

∼時間操縱緒論期末考

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 59.115.201.126

[問題] mac texshop上的引用圖形

作者: MIZUYAMA (致命祈願) 看板: LaTeX
標題: [問題] mac texshop上的引用圖形
時間: Sun Aug 15 23:02:40 2010

爬了文跟翻了手邊的教學書還是搞不定

目前手邊的工作是用R語言畫出圖後儲存(pdf檔)
在引用pdf圖檔

所以我用了
\usepackage[dvipdfm]{graphicx}
.. ...
.. ...
\begin{figure}
\centering

\includegraphics[width=0.9\textwidth]{R圖檔名.pdf}
\end{figure}


但是就是沒辦法跑出圖來 囧
想請問高手該怎麼做....

--
大學院研究的主題是時間∼直到他們將時間用盡

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 59.115.201.126
推 chungyuandye:menu=>typeset=>pdftex 不過還是建議轉成eps  08/15 23:34

2010年8月6日 星期五

(af2486) Re: [程式] 有關mathematica求解的一個問題

作者: chungyuandye (養花種魚數月亮賞星星)
標題: Re: [程式] 有關mathematica求解的一個問題
時間: Thu Aug 5 23:03:49 2010

※ 引述《af2486 (我喜歡旅行)》之銘言:
: 不好意思...找不到部落格要在哪邊留言直接回信囉
: 關於您新寫的程式碼我跑過確定可以用沒問題囉,非常感謝!
: 然後為求方便我把資料輸出成notepad檔(方便整個複製)
: NestList[{#[[1]] + 0.001, #[[3]], Re@y} /.
: FindRoot[
: f[y*(2 Pi/(#[[1]] + 0.01)), opt[#[[1]] + 0.01]] ==
: 0, {y, #[[3]]}] &, {0.699, 1.6, 1.4435}, 901][[
: 2 ;; -1]] >>> temp17
: 像這樣子。然後跑出來的檔案形式是像這樣(列舉部分資料)
: {{0.7, 1.4435, 1.4435871120952857}, {0.701, 1.4435871120952857,
: 1.443560146928692}, {0.702, 1.443560146928692, 1.4435045255446237},
: {0.703, 1.4435045255446237, 1.4434489234482561},
: {0.704, 1.4434489234482561, 1.4433933402197179}, ......
: 可是我只想複製每個最後算出來的結果貼到excel(標紅色的數字)
: 像目前這樣的形式就不能讓我一次複製一大筆資料而要一個一個貼
: 請問要怎麼改讓跑出來的資料就只有我想要的那些資料呢?
: 謝謝!

temp17=NestList[{#[[1]] + 0.001, #[[3]], Re@y} /.
FindRoot[
f[y*(2 Pi/(#[[1]] + 0.01)), opt[#[[1]] + 0.01]] ==
0, {y, #[[3]]}] &, {0.699, 1.6, 1.4435}, 901][[
2 ;; -1]]

Export["~/temp17.csv",temp17[[All,3]]]

~/temp17.csv是你的絕對路徑,Win下面你可能要自己改一下
接著用Excel開就好了

: ※ 引述《chungyuandye (養花種魚數月亮賞星星)》之銘言:
: : 有問題可以到我blog留言討論
: : xiLimit = 15.;
: : Hankel1[m_, x_] :=
: : If[Abs[x] < xiLimit, BesselJ[m, x] + I BesselY[m, x],
: : Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
: : BesselJp[m_, x_] :=
: : If[Abs[x] < xiLimit,
: : D[BesselJ[m, x1],
: : x1] /. {x1 -> x}, -Sqrt[2/(Pi*x)] Sin[(x - m/2*Pi - Pi/4)]];
: : BesselYp[m_, x_] :=
: : If[Abs[x] < xiLimit, D[BesselY[m, x1], x1] /. {x1 -> x},
: : Sqrt[2/(Pi*x)] Cos[(x - m/2*Pi - Pi/4)]];
: : Hankel1p[m_, x_] :=
: : If[Abs[x] < xiLimit, D[Hankel1[m, x1], x1] /. {x1 -> x},
: : I*Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
: : kr[kz_, k_] := Sqrt[k^2 - kz^2];
: : \[Alpha]r[kz_, k_] := Sqrt[kz^2 - k^2];
: : f[kz_, opt_] :=
: : Module[{sub}, k0 = \[Omega] Sqrt[\[Epsilon]0 \[Mu]0] /. opt;
: : k1 = \[Omega] Sqrt[\[Epsilon]1 \[Mu]1] /. opt;
: : k2 = \[Omega] Sqrt[\[Epsilon]2 \[Mu]2] /. opt;
: : r0 = r0 /. opt;
: : r1 = r1 /. opt;
: : m = m /. opt;
: : kr0 = kr[kz, k0];
: : kr1 = kr[kz, k1];
: : \[Alpha]r2 = \[Alpha]r[kz, k2];
: : BJ00 = BesselJ[m, kr[kz, k0] r0];
: : BJ10 = BesselJ[m, kr[kz, k1] r0];
: : BJ11 = BesselJ[m, kr[kz, k1] r1];
: : BY10 = BesselY[m, kr[kz, k1] r0];
: : BY11 = BesselY[m, kr[kz, k1] r1];
: : BJp00 = BesselJp[m, kr[kz, k0] r0];
: : BJp10 = BesselJp[m, kr[kz, k1] r0];
: : BJp11 = BesselJp[m, kr[kz, k1] r1];
: : BYp10 = BesselYp[m, kr[kz, k1] r0];
: : BYp11 = BesselYp[m, kr[kz, k1] r1];
: : HA21 = Hankel1[m, I \[Alpha]r[kz, k2] r1];
: : HAp21 = Hankel1p[m, I \[Alpha]r[kz, k2] r1];
: : M = {{BJ00, 0, -BJ10, -BY10, 0, 0, 0, 0}, {0, BJ00, 0,
: : 0, -BJ10, -BY10, 0,
: : 0}, {(\[Omega] \[Epsilon]0)/(kr0 r0) BJp00, (m kz)/(kr0^2 r0^2) \
: : BJ00, -((\[Omega] \[Epsilon]1)/(kr1 r0)) BJp10, -((\[Omega] \
: : \[Epsilon]1)/(kr1 r0)) BYp10, -((m kz)/(kr1^2 r0^2)) BJ10, -((m \
: : kz)/(kr1^2 r0^2)) BY10, 0,
: : 0}, {(m kz)/(kr0^2 r0^2) BJ00, (\[Omega] \[Mu]0)/(kr0 r0) \
: : BJp00, -((m kz)/(kr1^2 r0^2)) BJ10, -((m kz)/(kr1^2 r0^2)) BY10, -((\
: : \[Omega] \[Mu]1)/(kr1 r0)) BJp10, -((\[Omega] \[Mu]1)/(kr1 r0)) BYp10,
: : 0, 0}, {0, 0, BJ11, BY11, 0, 0, -HA21, 0}, {0, 0, 0, 0, BJ11,
: : BY11, 0, -HA21}, {0,
: : 0, (\[Omega] \[Epsilon]1)/(kr1 r1) BJp11, (\[Omega] \
: : \[Epsilon]1)/(kr1 r1) BYp11, (m kz)/(kr1^2 r1^2) BJ11, (m kz)/(kr1^2 \
: : r1^2) BY11, (I \[Omega] \[Epsilon]2)/(\[Alpha]r2 r1) HAp21, (m kz)/(\
: : \[Alpha]r2^2 r1^2) HA21}, {0,
: : 0, (m kz)/(kr1^2 r1^2) BJ11, (m kz)/(kr1^2 r1^2) BY11, \
: : (\[Omega] \[Mu]1)/(kr1 r1) BJp11, (\[Omega] \[Mu]1)/(kr1 r1) BYp11, \
: : (m kz)/(\[Alpha]r2^2 r1^2) HA21, (I \[Omega] \[Mu]2)/(\[Alpha]r2 r1) \
: : HAp21}} /. opt;
: : equ = Det[M];
: : Return[equ];]
: : n0[lambda_] :=
: : Sqrt[1 + (0.68671749*lambda^2)/(lambda^2 -
: : 0.072675189^2) + (0.43481505*lambda^2)/(lambda^2 -
: : 0.11514351^2) + (0.89656582*lambda^2)/(lambda^2 - 10.002398^2)];
: : n1[lambda_] :=
: : Sqrt[1 + (0.6961663*lambda^2)/(lambda^2 - 0.0684043^2) + (0.4079426*
: : lambda^2)/(lambda^2 - 0.1162414^2) + (0.8974794*
: : lambda^2)/(lambda^2 - 9.896161^2)];
: : n2[lambda_] := 1.0;
: : opt[lambda_] := {m -> 1, \[Omega] -> 2 Pi/lambda, r0 -> 4.1,
: : r1 -> 62.5, \[Epsilon]0 -> n0[lambda]^2, \[Mu]0 ->
: : 1., \[Epsilon]1 -> n1[lambda]^2, \[Mu]1 -> 1., \[Epsilon]2 ->
: : n2[lambda]^2, \[Mu]2 -> 1., kz -> y};
: : NestList[{#[[1]] + 0.01, #[[3]], Re@y} /.
: : FindRoot[
: : f[y*(2 Pi/(#[[1]] + 0.01)), opt[#[[1]] + 0.01]] ==
: : 0, {y, #[[3]]}] &, {1.53, 1.4435, 1.4435},
: : 20][[2 ;; -1]] // TableForm

--
養花種魚數月亮賞星星

http://cydye1069.blogspot.com

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.129.12

2010年8月5日 星期四

[自動轉寄] (af2486) Re: [程式] 有關mathematica求解的一個問題

作者: chungyuandye (養花種魚數月亮賞星星)
標題: Re: [程式] 有關mathematica求解的一個問題
時間: Thu Aug 5 23:03:49 2010

※ 引述《af2486 (我喜歡旅行)》之銘言:
: 不好意思...找不到部落格要在哪邊留言直接回信囉
: 關於您新寫的程式碼我跑過確定可以用沒問題囉,非常感謝!
: 然後為求方便我把資料輸出成notepad檔(方便整個複製)
: NestList[{#[[1]] + 0.001, #[[3]], Re@y} /.
: FindRoot[
: f[y*(2 Pi/(#[[1]] + 0.01)), opt[#[[1]] + 0.01]] ==
: 0, {y, #[[3]]}] &, {0.699, 1.6, 1.4435}, 901][[
: 2 ;; -1]] >>> temp17
: 像這樣子。然後跑出來的檔案形式是像這樣(列舉部分資料)
: {{0.7, 1.4435, 1.4435871120952857}, {0.701, 1.4435871120952857,
: 1.443560146928692}, {0.702, 1.443560146928692, 1.4435045255446237},
: {0.703, 1.4435045255446237, 1.4434489234482561},
: {0.704, 1.4434489234482561, 1.4433933402197179}, ......
: 可是我只想複製每個最後算出來的結果貼到excel(標紅色的數字)
: 像目前這樣的形式就不能讓我一次複製一大筆資料而要一個一個貼
: 請問要怎麼改讓跑出來的資料就只有我想要的那些資料呢?
: 謝謝!

temp17=NestList[{#[[1]] + 0.001, #[[3]], Re@y} /.
FindRoot[
f[y*(2 Pi/(#[[1]] + 0.01)), opt[#[[1]] + 0.01]] ==
0, {y, #[[3]]}] &, {0.699, 1.6, 1.4435}, 901][[
2 ;; -1]]

Export["~/temp17.csv",temp17[[All,3]]]

~/temp17.csv是你的絕對路徑,Win下面你可能要自己改一下
接著用Excel開就好了

: ※ 引述《chungyuandye (養花種魚數月亮賞星星)》之銘言:
: : 有問題可以到我blog留言討論
: : xiLimit = 15.;
: : Hankel1[m_, x_] :=
: : If[Abs[x] < xiLimit, BesselJ[m, x] + I BesselY[m, x],
: : Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
: : BesselJp[m_, x_] :=
: : If[Abs[x] < xiLimit,
: : D[BesselJ[m, x1],
: : x1] /. {x1 -> x}, -Sqrt[2/(Pi*x)] Sin[(x - m/2*Pi - Pi/4)]];
: : BesselYp[m_, x_] :=
: : If[Abs[x] < xiLimit, D[BesselY[m, x1], x1] /. {x1 -> x},
: : Sqrt[2/(Pi*x)] Cos[(x - m/2*Pi - Pi/4)]];
: : Hankel1p[m_, x_] :=
: : If[Abs[x] < xiLimit, D[Hankel1[m, x1], x1] /. {x1 -> x},
: : I*Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
: : kr[kz_, k_] := Sqrt[k^2 - kz^2];
: : \[Alpha]r[kz_, k_] := Sqrt[kz^2 - k^2];
: : f[kz_, opt_] :=
: : Module[{sub}, k0 = \[Omega] Sqrt[\[Epsilon]0 \[Mu]0] /. opt;
: : k1 = \[Omega] Sqrt[\[Epsilon]1 \[Mu]1] /. opt;
: : k2 = \[Omega] Sqrt[\[Epsilon]2 \[Mu]2] /. opt;
: : r0 = r0 /. opt;
: : r1 = r1 /. opt;
: : m = m /. opt;
: : kr0 = kr[kz, k0];
: : kr1 = kr[kz, k1];
: : \[Alpha]r2 = \[Alpha]r[kz, k2];
: : BJ00 = BesselJ[m, kr[kz, k0] r0];
: : BJ10 = BesselJ[m, kr[kz, k1] r0];
: : BJ11 = BesselJ[m, kr[kz, k1] r1];
: : BY10 = BesselY[m, kr[kz, k1] r0];
: : BY11 = BesselY[m, kr[kz, k1] r1];
: : BJp00 = BesselJp[m, kr[kz, k0] r0];
: : BJp10 = BesselJp[m, kr[kz, k1] r0];
: : BJp11 = BesselJp[m, kr[kz, k1] r1];
: : BYp10 = BesselYp[m, kr[kz, k1] r0];
: : BYp11 = BesselYp[m, kr[kz, k1] r1];
: : HA21 = Hankel1[m, I \[Alpha]r[kz, k2] r1];
: : HAp21 = Hankel1p[m, I \[Alpha]r[kz, k2] r1];
: : M = {{BJ00, 0, -BJ10, -BY10, 0, 0, 0, 0}, {0, BJ00, 0,
: : 0, -BJ10, -BY10, 0,
: : 0}, {(\[Omega] \[Epsilon]0)/(kr0 r0) BJp00, (m kz)/(kr0^2 r0^2) \
: : BJ00, -((\[Omega] \[Epsilon]1)/(kr1 r0)) BJp10, -((\[Omega] \
: : \[Epsilon]1)/(kr1 r0)) BYp10, -((m kz)/(kr1^2 r0^2)) BJ10, -((m \
: : kz)/(kr1^2 r0^2)) BY10, 0,
: : 0}, {(m kz)/(kr0^2 r0^2) BJ00, (\[Omega] \[Mu]0)/(kr0 r0) \
: : BJp00, -((m kz)/(kr1^2 r0^2)) BJ10, -((m kz)/(kr1^2 r0^2)) BY10, -((\
: : \[Omega] \[Mu]1)/(kr1 r0)) BJp10, -((\[Omega] \[Mu]1)/(kr1 r0)) BYp10,
: : 0, 0}, {0, 0, BJ11, BY11, 0, 0, -HA21, 0}, {0, 0, 0, 0, BJ11,
: : BY11, 0, -HA21}, {0,
: : 0, (\[Omega] \[Epsilon]1)/(kr1 r1) BJp11, (\[Omega] \
: : \[Epsilon]1)/(kr1 r1) BYp11, (m kz)/(kr1^2 r1^2) BJ11, (m kz)/(kr1^2 \
: : r1^2) BY11, (I \[Omega] \[Epsilon]2)/(\[Alpha]r2 r1) HAp21, (m kz)/(\
: : \[Alpha]r2^2 r1^2) HA21}, {0,
: : 0, (m kz)/(kr1^2 r1^2) BJ11, (m kz)/(kr1^2 r1^2) BY11, \
: : (\[Omega] \[Mu]1)/(kr1 r1) BJp11, (\[Omega] \[Mu]1)/(kr1 r1) BYp11, \
: : (m kz)/(\[Alpha]r2^2 r1^2) HA21, (I \[Omega] \[Mu]2)/(\[Alpha]r2 r1) \
: : HAp21}} /. opt;
: : equ = Det[M];
: : Return[equ];]
: : n0[lambda_] :=
: : Sqrt[1 + (0.68671749*lambda^2)/(lambda^2 -
: : 0.072675189^2) + (0.43481505*lambda^2)/(lambda^2 -
: : 0.11514351^2) + (0.89656582*lambda^2)/(lambda^2 - 10.002398^2)];
: : n1[lambda_] :=
: : Sqrt[1 + (0.6961663*lambda^2)/(lambda^2 - 0.0684043^2) + (0.4079426*
: : lambda^2)/(lambda^2 - 0.1162414^2) + (0.8974794*
: : lambda^2)/(lambda^2 - 9.896161^2)];
: : n2[lambda_] := 1.0;
: : opt[lambda_] := {m -> 1, \[Omega] -> 2 Pi/lambda, r0 -> 4.1,
: : r1 -> 62.5, \[Epsilon]0 -> n0[lambda]^2, \[Mu]0 ->
: : 1., \[Epsilon]1 -> n1[lambda]^2, \[Mu]1 -> 1., \[Epsilon]2 ->
: : n2[lambda]^2, \[Mu]2 -> 1., kz -> y};
: : NestList[{#[[1]] + 0.01, #[[3]], Re@y} /.
: : FindRoot[
: : f[y*(2 Pi/(#[[1]] + 0.01)), opt[#[[1]] + 0.01]] ==
: : 0, {y, #[[3]]}] &, {1.53, 1.4435, 1.4435},
: : 20][[2 ;; -1]] // TableForm

--
養花種魚數月亮賞星星

http://cydye1069.blogspot.com

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.129.12

[自動轉寄] Re: [程式] 有關mathematica求解的一個問題

作者: af2486 (我喜歡旅行)
標題: Re: [程式] 有關mathematica求解的一個問題
時間: Thu Aug 5 22:01:05 2010

不好意思...找不到部落格要在哪邊留言直接回信囉
關於您新寫的程式碼我跑過確定可以用沒問題囉,非常感謝!
然後為求方便我把資料輸出成notepad檔(方便整個複製)
NestList[{#[[1]] + 0.001, #[[3]], Re@y} /.
FindRoot[
f[y*(2 Pi/(#[[1]] + 0.01)), opt[#[[1]] + 0.01]] ==
0, {y, #[[3]]}] &, {0.699, 1.6, 1.4435}, 901][[
2 ;; -1]] >>> temp17

像這樣子。然後跑出來的檔案形式是像這樣(列舉部分資料)
{{0.7, 1.4435, 1.4435871120952857}, {0.701, 1.4435871120952857,
1.443560146928692}, {0.702, 1.443560146928692, 1.4435045255446237},
{0.703, 1.4435045255446237, 1.4434489234482561},
{0.704, 1.4434489234482561, 1.4433933402197179}, ......

可是我只想複製每個最後算出來的結果貼到excel(標紅色的數字)
像目前這樣的形式就不能讓我一次複製一大筆資料而要一個一個貼
請問要怎麼改讓跑出來的資料就只有我想要的那些資料呢?
謝謝!


※ 引述《chungyuandye (養花種魚數月亮賞星星)》之銘言:
: 有問題可以到我blog留言討論
: xiLimit = 15.;
: Hankel1[m_, x_] :=
: If[Abs[x] < xiLimit, BesselJ[m, x] + I BesselY[m, x],
: Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
: BesselJp[m_, x_] :=
: If[Abs[x] < xiLimit,
: D[BesselJ[m, x1],
: x1] /. {x1 -> x}, -Sqrt[2/(Pi*x)] Sin[(x - m/2*Pi - Pi/4)]];
: BesselYp[m_, x_] :=
: If[Abs[x] < xiLimit, D[BesselY[m, x1], x1] /. {x1 -> x},
: Sqrt[2/(Pi*x)] Cos[(x - m/2*Pi - Pi/4)]];
: Hankel1p[m_, x_] :=
: If[Abs[x] < xiLimit, D[Hankel1[m, x1], x1] /. {x1 -> x},
: I*Sqrt[2/(Pi*x)] Exp[I*(x - m/2*Pi - Pi/4)]];
: kr[kz_, k_] := Sqrt[k^2 - kz^2];
: \[Alpha]r[kz_, k_] := Sqrt[kz^2 - k^2];
: f[kz_, opt_] :=
: Module[{sub}, k0 = \[Omega] Sqrt[\[Epsilon]0 \[Mu]0] /. opt;
: k1 = \[Omega] Sqrt[\[Epsilon]1 \[Mu]1] /. opt;
: k2 = \[Omega] Sqrt[\[Epsilon]2 \[Mu]2] /. opt;
: r0 = r0 /. opt;
: r1 = r1 /. opt;
: m = m /. opt;
: kr0 = kr[kz, k0];
: kr1 = kr[kz, k1];
: \[Alpha]r2 = \[Alpha]r[kz, k2];
: BJ00 = BesselJ[m, kr[kz, k0] r0];
: BJ10 = BesselJ[m, kr[kz, k1] r0];
: BJ11 = BesselJ[m, kr[kz, k1] r1];
: BY10 = BesselY[m, kr[kz, k1] r0];
: BY11 = BesselY[m, kr[kz, k1] r1];
: BJp00 = BesselJp[m, kr[kz, k0] r0];
: BJp10 = BesselJp[m, kr[kz, k1] r0];
: BJp11 = BesselJp[m, kr[kz, k1] r1];
: BYp10 = BesselYp[m, kr[kz, k1] r0];
: BYp11 = BesselYp[m, kr[kz, k1] r1];
: HA21 = Hankel1[m, I \[Alpha]r[kz, k2] r1];
: HAp21 = Hankel1p[m, I \[Alpha]r[kz, k2] r1];
: M = {{BJ00, 0, -BJ10, -BY10, 0, 0, 0, 0}, {0, BJ00, 0,
: 0, -BJ10, -BY10, 0,
: 0}, {(\[Omega] \[Epsilon]0)/(kr0 r0) BJp00, (m kz)/(kr0^2 r0^2) \
: BJ00, -((\[Omega] \[Epsilon]1)/(kr1 r0)) BJp10, -((\[Omega] \
: \[Epsilon]1)/(kr1 r0)) BYp10, -((m kz)/(kr1^2 r0^2)) BJ10, -((m \
: kz)/(kr1^2 r0^2)) BY10, 0,
: 0}, {(m kz)/(kr0^2 r0^2) BJ00, (\[Omega] \[Mu]0)/(kr0 r0) \
: BJp00, -((m kz)/(kr1^2 r0^2)) BJ10, -((m kz)/(kr1^2 r0^2)) BY10, -((\
: \[Omega] \[Mu]1)/(kr1 r0)) BJp10, -((\[Omega] \[Mu]1)/(kr1 r0)) BYp10,
: 0, 0}, {0, 0, BJ11, BY11, 0, 0, -HA21, 0}, {0, 0, 0, 0, BJ11,
: BY11, 0, -HA21}, {0,
: 0, (\[Omega] \[Epsilon]1)/(kr1 r1) BJp11, (\[Omega] \
: \[Epsilon]1)/(kr1 r1) BYp11, (m kz)/(kr1^2 r1^2) BJ11, (m kz)/(kr1^2 \
: r1^2) BY11, (I \[Omega] \[Epsilon]2)/(\[Alpha]r2 r1) HAp21, (m kz)/(\
: \[Alpha]r2^2 r1^2) HA21}, {0,
: 0, (m kz)/(kr1^2 r1^2) BJ11, (m kz)/(kr1^2 r1^2) BY11, \
: (\[Omega] \[Mu]1)/(kr1 r1) BJp11, (\[Omega] \[Mu]1)/(kr1 r1) BYp11, \
: (m kz)/(\[Alpha]r2^2 r1^2) HA21, (I \[Omega] \[Mu]2)/(\[Alpha]r2 r1) \
: HAp21}} /. opt;
: equ = Det[M];
: Return[equ];]
: n0[lambda_] :=
: Sqrt[1 + (0.68671749*lambda^2)/(lambda^2 -
: 0.072675189^2) + (0.43481505*lambda^2)/(lambda^2 -
: 0.11514351^2) + (0.89656582*lambda^2)/(lambda^2 - 10.002398^2)];
: n1[lambda_] :=
: Sqrt[1 + (0.6961663*lambda^2)/(lambda^2 - 0.0684043^2) + (0.4079426*
: lambda^2)/(lambda^2 - 0.1162414^2) + (0.8974794*
: lambda^2)/(lambda^2 - 9.896161^2)];
: n2[lambda_] := 1.0;
: opt[lambda_] := {m -> 1, \[Omega] -> 2 Pi/lambda, r0 -> 4.1,
: r1 -> 62.5, \[Epsilon]0 -> n0[lambda]^2, \[Mu]0 ->
: 1., \[Epsilon]1 -> n1[lambda]^2, \[Mu]1 -> 1., \[Epsilon]2 ->
: n2[lambda]^2, \[Mu]2 -> 1., kz -> y};
: NestList[{#[[1]] + 0.01, #[[3]], Re@y} /.
: FindRoot[
: f[y*(2 Pi/(#[[1]] + 0.01)), opt[#[[1]] + 0.01]] ==
: 0, {y, #[[3]]}] &, {1.53, 1.4435, 1.4435},
: 20][[2 ;; -1]] // TableForm
: ※ 引述《af2486 (我喜歡旅行)》之銘言:
: : 您好~不好意思我想請問您的程式碼是要從哪邊開始貼呢?
: : 或是說取代哪邊的程式碼...我取代"Do"這個指令後面的
: : 要執行的時候好像有點小問題(出在函數的引數[]那邊,有個錯誤訊息)
: : 不好意思才剛開始學連基本的debug都不太會...這小問題可能還要麻煩你了謝謝!

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.42.126