symbolic integration needed in the present code
1 回表示 (過去 30 日間)
古いコメントを表示
syms h q x f(x) f0 f1(x) f2(x) f3(x)
f = f0 + f1*q + f2*q.^2 + f3*q.^3; Eq0 = diff(f,x) - diff(f0,x); Eq1 = diff(f,x) - f^2;
A = (1-q)*(Eq0) == h*q*Eq1; B = diff(A,q); C = subs(B,[q,f0],[0,1]);
D = int(C,x)
%%% D: f1 = h*x + Constant of integration; This answer should come when I run this code but ...
0 件のコメント
回答 (1 件)
Paul
2022 年 4 月 26 日
syms h q x f(x) f0 f1(x) f2(x) f3(x)
f = f0 + f1*q + f2*q.^2 + f3*q.^3;
Eq0 = diff(f,x) - diff(f0,x);
Eq1 = diff(f,x) - f^2;
A = (1-q)*(Eq0) == h*q*Eq1;
B = diff(A,q);
C = subs(B,[q,f0],[0,1])
As stated on its doc page, int() does not return a constant of integration in its result, so it has to be added explicitly:
syms C1
D = int(lhs(C),x) == int(rhs(C),x) + C1
参考
カテゴリ
Help Center および File Exchange で Symbolic Math Toolbox についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!