| 查看: 414 | 回复: 1 | ||
am99lx银虫 (小有名气)
|
[求助]
求大神帮我看看我这代码为什么出不了图，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 |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有20人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有118人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
am99lx
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 318.1
- 帖子: 59
- 在线: 47.1小时
- 虫号: 4105406
- 注册: 2015-09-27
- 性别: GG
- 专业: 电工材料特性及其应用
|
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)}] 改成这样 还是不出图 求大神指点 |
2楼2015-10-22 09:16:11
