Need some change in the present Analytical solution code which I'm unable
이전 댓글 표시
syms x F1(x) G1(x) H1(x) n fw h
Nt = 0.5; Nb = 0.5; Pr = 1; Le = 2; DF = diff(F1,x);
eq1 = diff(F1,x,3) + diff(F1,x,2) == 0; eq2 = diff(G1,x,1) + diff(G1,x,2) == 0; eq3 = diff(H1,x,1) + diff(H1,x,2) == 0;
EQ = [eq1, eq2, eq3]; S = dsolve(EQ);f1 = S.F1;
F(1) = fw + 1 - exp(-x);
for m = 2:5
for k = 1:m-1
K = F(k)*diff(F(m-k),2) - (2*n/(n+1))*diff(F(k))*diff(F(m-k));
F(m) = f1 + F(m-1) + h*( diff(F(m-1),3) + int(int(int(K))) + int(int(K)) );
end
end
F = F(1)+F(2);F0 = subs(F,x,0); dF0 = subs(diff(F,x),x,0); dFI = subs(diff(F,x),x,Inf);
sv = solve(F, {F0 == fw, dF0 == 0,dFI == 0})
%% I want to simplify 'F' by putting the conditions {F0 == fw, dF0 == 0,dFI == 0}, so that the values of the Constants present in 'F' can be find out.
댓글 수: 3
P.K. Pattnaik
2022년 5월 10일
P.K. Pattnaik
2023년 11월 4일
Walter Roberson
2023년 11월 4일
use subs() to do the replacements before solving.
subs(F, {list of things to replace}, {list of things to replace with})
채택된 답변
추가 답변 (0개)
카테고리
도움말 센터 및 File Exchange에서 Assumptions에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
