function?p=fcatch（Pmdn）
%函數(shù)的求取過(guò)程較為復(fù)雜，應(yīng)該可以簡(jiǎn)化，可將矩陣中的元素提取出來(lái)，進(jìn)行相關(guān)的運(yùn)算
L=zeros（nn）;%定義下三角矩陣
U=zeros（nn）;%定義上三角矩陣
L（11）=Pm（11）;
U（11）=1;
for?i=2:n
????U（ii）=1;
????U（i-1i）=Pm（i-1i）/L（i-1i-1）;
????L（ii）=Pm（ii）-Pm（ii-1）*U（i-1i）;
????L（ii-1）=Pm（ii-1）;
end
%----------------------追的過(guò)程————————————————————%
y=zeros（n1）;
y（1）=d（1）/L（11）;
for?i=2:n
????y（i）=（d（i）-L（ii-1）*y（i-1））/L（ii）;
end
%—————————————趕得過(guò)程，求壓力—————————————————%
p=zeros（n1）;
p（n）=y（n）;
for?i=n-1:-1:1
????p（i）=y（i）-U（ii+1）*p（i+1）;
end

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件????????2347??2018-05-05?21:47??一維徑向流模擬\Main.m
?????文件????????1244??2014-06-28?23:21??一維徑向流模擬\Solve.m
?????文件?????????651??2014-05-18?16:18??一維徑向流模擬\fcatch.m
?????目錄???????????0??2018-05-05?21:36??一維徑向流模擬\