function??[AUDLrbx]=?xianghoucafen（abn）
%輸入時:主意?阿爾法和貝塔的值是要變化的??y‘‘?+p（x）*y‘+q（x）*y=f（x）
%?a?是初始值?y（a）=阿爾法，b是最終值?y（b）=貝塔??，n是分段數(shù)
format?long
h?=?0.5;%??有些時候?需要輸入?例如?等于?0
disp（h）;
A=zeros（n-1n-1）;
b=zeros（n-11）;
syms?x
q?=[];
p=[];
r=[];
for?i?=?1:?n-1?
????x=0+（i+1）?*?h;
????p（i）=0;???%?變化?請更改成為?x?的函數(shù)?例子?:p（i）?=-x/（1+x^2）
????q（i）?=（-5/1458）;??%?變化?請更改成為?x?的函數(shù)
????r（i）=（x^2-5*x）/（2916）;?%?變化?請更改成為?x?的函數(shù)
????if?i?==1
????????A（11）=-2?+h*p（1）;
????????A（12）=?（1?+h*p（1）+（h^2）*?q（1）?）?;
????????b（11）=（h^2）*?r（i）-0;?%??請更改??最后一個數(shù)字是?初值的阿爾法
????elseif?i?????????A（ii-1）=1?;
????????A（ii）=-2?-?h?*?p（i）;
????????A（ii+1）=1+h*p（i）+（h^2）*q（i）?;
????????b（i1）=（h^2）?*?r?（i）;
????el