MATLAB的答案

我想解决常微分方程采用数值有限体积方法造成的。但是我错误的代码。请帮助纠正错误。

2视图(30天)
(sanjeev kumar jha
(sanjeev kumar jha 2023年4月9日
评论道: Torsten 2023年4月9日
Lx = 0.1;
Ly = 0.1;
dx = 0.01;
dy = 0.01;
nx = Lx / dx;
纽约= Ly / dy;
锡= 500;
x = dx / 2: dx: Lx + (dx / 2);
y = dy / 2: dy: Ly + (dy / 2);
Nx = Nx + 1;
纽约=纽约+ 1;
M = 0 (Nx、纽约);
Tini = 500;
T = Tini * (1 + M);
告诉= T;
Counter = 0;
因为我= 1:Nx
j = 1:纽约
如果(i = = 1) & & (j = = 1)
dvdt = @ (t, W)[3 *锡+ 4 * W (i, j) 3 * W (i, j + 1)];
elseif (j = = 1)
dvdt = @ (t, W) [2 * W(张,j) + 4 * W (i + 1, j + 1) 3 * W (i, j)];
elseif (i = = 1)
dvdt = @ (t, W) [3 * 2 * W锡(i, j + 1) + 4 * 3 W (i, j - 1) * W (i, j)];
elseif (i = = 1) & & (j = =纽约)
dvdt = @ (t, W) [3 * Tin-2 * W (i, j) + 4 * W (i, j - 1)];
elseif (j = =纽约)
dvdt = @ (t, W) [2 * W(张,j) + 4 * W (i, j - 1) 3 * W (i, j)];
其他的
dvdt = @ (t, W) [W(张,j) + 2 * 4 * 3 W (i, j + 1) * W (i, j - 1)];
结束
[t W] =数值(0.4 [0]dvdt告诉)
结束
结束
y = W (: 1);
% z = W (:, 2);
情节(t, W)
2的评论
显示1年长的评论

登录置评。

在回答这个问题。

答案(0)

类别

找到更多的在常微分方程在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!