假設有耦合PDE的一個系統,諸如 第一PDE在F(A,B) 第二PDE在F(A,B) 以下代碼能夠分別解決每個PDE: import numpy as np
import sympy as sp
# definition of variables
a, b = sp.symbols('a b')
f = sp.Function('f')
F = f(a, b)
Fda = F.dif
我試圖用OpenModelica來數值求解一個非常簡單的PDE du/dx = du/dt,邊界條件爲u(0,t)= t^2和u_x(0,t) = 0。我已經寫了下面的代碼: model pdetest_1
parameter Real L=1;
parameter Integer N=100;
parameter Real dx=L/(N-1);
par