function??[AbUDLx]=?zhongxincafen（abn）
%輸入時:主意?阿爾法和貝塔的值是要變化的??y‘‘?+p（x）*y‘+q（x）*y=f（x）
%?a?是初始值?y（a）=阿爾法，b是最終值?y（b）=貝塔??，n是分段數
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*h;
????p（i）=0;???%?變化?請更改成為?x?的函數?例子?:p（i）?=-x/（1+x^2）
????q（i）?=（-5/1458）;??%?變化?請更改成為?x?的函數
????r（i）=（x^2-5*x）/（2916）;???%?變化?請更改成為?x?的函數
????if?i?==?1
????????A（11）=-2?+?h^2?*?q（i）;
????????A（12）=?1?+?0.5?*?h?*?p（i）;
????????b（11）=h^2*?r（i）-（1-0.5*h*p（i））*0;?%??請更改??最后一個數字是?初值的阿爾法
????elseif?i<=n-2
????????A（ii-1）=1-0.5*h*p（i）;
????????A（ii）=-2+h^2*q（i）;
????????A（ii+1）=1+0.5*?h?*p（i）;
????????b（i1）=h^2?*r（i）;
????el