技術標籤：MATLAB學習筆記matlab

自定義函式

function [ R, q, err0, err1, h ] = e( f, a, b, n, tol0)
M = 1;
h = b-a;
y = integral(f,a,b);
err0 = 1;
err1 = 1;
J = 0;
R = zeros(n,n);
R(1,1) = h*(feval(f,a)+feval(f,b))/2;
while (err0>tol0)&&(J<n)
    J = J+1;
    h = h/2;
    s = 0;
    for p = 1:M
        x =
 
 a+h*(2*p-1);
        s = s+feval(f,x);
    end
    R(J+1, 1) = R(J, 1)/2+h*s;
    M = 2*M;
    for K=1:J
        R(J+1, K+1) = R(J+1, K)+(R(J+1,K)-R(J,K))/(4^K-1);
    end
    err1 = abs(R(J,J)-R(J+1,K+1));
    err0 = abs(R(J+1,K+1)-y);
end
q = R(J+1,J+1);
end

將該函式作為引數賦值給新的函式中

eps =1e-12;
a = 0;
b = 1
 
;
n = 8;
f1=@f;
[ R, q, err0, err1, h ] = e(f1, a, b, n, tol0)

多引數的函式

function [ output ] = fun(a,b,c )
output = a + b + c;
end

function [ output ] = fun1(f,a,b,c,d)
output = f(a,b,c)*d;
end

f=@fun;
fun1(f,1,2,3,4)