| 查看: 478 | 回复: 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 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有265人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有87人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
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
