Need dsolve code to run for Symbolic solution
1 次查看(过去 30 天)
显示 更早的评论
syms U(y) T(y) C(y) M Kp Pr phi Gr Gc Sc Kc a
% % % Gr = 1; Gc = 1; M =1; Kp = 0.1; Pr = 1; phi = 0.1; Sc = 0.22; Kc = 1;
xa = 0; xb = Inf;
Cond = [U(0) == a/s^2; T(0) == 1/s^2; C(0) == 1/(s-1); U(Inf) == 0; T(Inf) == 0; C(Inf) == 0];
S = dsolve( [diff(U,2) + Gc*C + Gr*T == U*(Kp + M + s), diff(T,2) == Pr*T*(phi + s), diff(C,2) == Sc*(Kc + s)*C], Cond);
U = S.U; T = S.T; C = S.C;
%% Code didn't run for INFINITY boundary condition
%% Can there be any other way to put INFINITY boundary condition after getting the general solution U, T, and C
%% I need to find inverse LAPLACE TRANSFORM of U, T, and C
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!