这是弄的一个牛顿迭代的程序,之前都没出现问题,这一下就出现了无法转换的问题,有点懵,麻烦各位给点意见,谢谢。
x是符号,不是值,所以根据x计算出的T、c也是一个表达式而不是一个具体的值,所以没法和一个具体的值1e-6来比较,要在使用前先把x转换为求解的值xk并计算才能到的一个具体的值。另外,注意到程序中在while里并没有更新c的值,我也添加了一下。
clear;
clc;
syms x;
y = (exp((-x^2)/2))*((2*pi)^(-1/2));
xk=0.5;
i=0;
T=int(y,x,0,xk);
Tk = eval(subs(T,x,xk));
c=abs(Tk-0.45);
while c>=1e-6
xk=xk-(T-0.45)/((exp((-xk^2)/2))*((2*pi)^(-1/2)));
i=i+1;
T=int(y,x,0,xk);
Tk = eval(subs(T,x,xk));
c=abs(Tk-0.45);
end