 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 406 | 回復(fù): 1 | |||
am99lx銀蟲(chóng) (小有名氣)
|
[求助]
求大神幫我看看我這代碼為什么出不了圖，mathematics
|
|
a = 1.2*10^7; b = 1.23*10^11; c = 2472960; d = 1.53*10^(-8); e = 6.42*10^4; f = 1.4*10^(-13); g = 2.6*10^16; h = 4*Log[2]; i = 4700; NDSolve[{70*M[z, t]*D[M[z, t], t, t] == -D[F[z, t], z, z] - g*(M[z, t] - N[z, t]) + Sqrt[h/3.1415]*0.07*i/f/d*E^(-z/d - h*((t - 2*f)/f)^2), 1/(a*M[z, t]^2 + b*N[z, t])*D[F[z, t], t, t] + F[z, t] == (-353*((M[z, t]/e)^2 + 0.16)^(4/5)*((M[z, t]/e)^2 + 0.44) (M[z, t]/e)^2)/((M[z, t]/e)^2 + 0.092)^(1/ 2)*((M[z, t]/e)^2 + 0.16*N[z, t]/e)* D[M[z, t], z, z], (c - 70*M[z, t])*D[M[z, t], t, t] == g*(M[z, t] - N[z, t]), M[z, 0] == 300, N[z, 0] == 300, F[z, 0] == 0, F[0, t] == 0, F[d, t] == 0}, {M[z, t], N[z, t], F[z, t]}, {z, 0, d}, {t, 0, 3*10^(-11)}] Plot3D[M[z, t], {z, 0, 1.53*10^(-8)}, {t, 0, 3*10^(-11)}] X9{{NRI}KP]6X[PHY0ZP}0Y.png |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有196人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有68人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
am99lx
銀蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 318.1
- 帖子: 59
- 在線: 47.1小時(shí)
- 蟲(chóng)號(hào): 4105406
- 注冊(cè): 2015-09-27
- 性別: GG
- 專(zhuān)業(yè): 電工材料特性及其應(yīng)用
|
a = 1.2*10^7; b = 1.23*10^11; c = 2472960; d = 1.53*10^(-8); e = 6.42*10^4; f = 1.4*10^(-13); g = 2.6*10^16; h = 4*Log[2]; i = 4700; DSolve[{70*M[z, t]*D[M[z, t], t, t] == -D[F[z, t], z, z] - g*(M[z, t] - Q[z, t]) + Sqrt[h/3.1415]*0.07*i/f/d*E^(-z/d - h*((t - 2*f)/f)^2), 1/(a*M[z, t]^2 + b*Q[z, t])*D[F[z, t], t, t] + F[z, t] == (-353*((M[z, t]/e)^2 + 0.16)^(4/5)*((M[z, t]/e)^2 + 0.44) (M[z, t]/e)^2)/((M[z, t]/e)^2 + 0.092)^(1/ 2)*((M[z, t]/e)^2 + 0.16*Q[z, t]/e)* D[M[z, t], z, z], (c - 70*M[z, t])*D[Q[z, t], t, t] == g*(M[z, t] - Q[z, t]), M[z, 0] == 300, N[z, 0] == 300, F[z, 0] == 0, F[0, t] == 0, F[d, t] == 0}, {M[z, t], Q[z, t], F[z, t]}, {z, 0, d}, {t, 0, 3*10^(-11)}] Plot3D[M[z, t], z -> 0, {t, 0, 3*10^(-11)}] 改成這樣 還是不出圖 求大神指點(diǎn) |
2樓2015-10-22 09:16:11
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 am99lx 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
