MATLAB求一道定积分方程已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/28 03:32:09
MATLAB求一道定积分方程已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.

MATLAB求一道定积分方程已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.
MATLAB求一道定积分方程

已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.

MATLAB求一道定积分方程已知 τ = 0.65,要求出 θ 的值.最好能给出代码和结果.
clc;clear
format long
f=@(t) t./(exp(t)-1);
t=@(x) 1+4./x.*(1./x.*quad(f,0,x)-1)-0.65;
x1=1;
x2=10;
tol=1e-6;
while abs(x1-x2)>tol
    if feval(t,x1)*feval(t,(x1+x2)/2)>0
        x1=(x1+x2)/2;
    elseif feval(t,x1)*feval(t,(x1+x2)/2)<0
        x2=(x1+x2)/2;
    else
        break
    end
end
 x=(x1+x2)/2;
disp('精确到小数点后6位的解为')
disp(x)
t=feval(t,x)+0.65;
disp('将x带入式子结果为')
disp(t)
结果为
精确到小数点后6位的解为
   9.437618821859360

将x带入式子结果为
   0.649999994395076