MATLAB 线性方程求解函数 - 扩展区间功能

function u = xianxing(h, v1, v2) format long f = @(t)((gamma(4)/gamma(3.6))t^2.6+3t^3+6); a = @(n,k)((n-k+1)^0.6-(n-k)^0.6); u = zeros(1,1/h+1); u(1) = 2; n = (v2-v1)/h; u(2) = (h^(-0.4)a(1,1)u(1)+gamma(1.6)f(v1+h))/(h^(-0.4)+3gamma(1.6)); for j=2:n an = 0; for i=1:j-1 an = an+(a(j,i+1)-a(j,i))u(i+1); end u(j+1) = (h^(-0.4)(an +a(j,1)u(1))+gamma(1.6)(f(v1+jh)))/(h^(-0.4)+3gamma(1.6)); end end