matlab向后差分

function  [AUDLrbx]= xianghoucafen（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+1） * 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*p（1）;
        A（12）= （1 +h*p（1）+（h^2）* q（1） ） ;
        b（11）=（h^2）* r（i）-0; %  请更改  最后一个数字是 初值的阿尔法
    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