| 查看: 381 | 回复: 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 |
回复此楼» 猜你喜欢
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有213人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有1人回复
-
投稿chemical physical letters不送审？
已经有2人回复
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
已经有22人回复
-
杂志要求提供official language-editing certificate
已经有4人回复
-
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有16人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
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