| 查看: 374 | 回复: 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 |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有111人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
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