szldh2005
2013-04-17, 13:58
现有一广义积分,2010a版可以得到结果,而2008a版则出错,不知原因。特此求救!
% 广义积分计算
% 20130327-20130416
clc;clear;
syms x y u b
u=10^-34
%%%%%%% B 符号积分->数值计算 %%%%%%
% 一、原无限积分
disp('%%%%%% B 符号积分->数值计算 %%%%%%%%')
disp('一、原无限积分')
b=-log(10^-15);
f1=@(x) exp(-x)./x;
w1_k=@(u) vpa(int(f1,x,u,inf),20); % 符号积分-数值计算(原无限积分)
w1d=@(u) vpa(int(f1,x,u,b),20); % 符号积分-数值计算(原无限截断积分)
w1_k=w1_k(u)
w1d=w1d(u)
% 广义积分计算
% 20130327-20130416
clc;clear;
syms x y u b
u=10^-34
%%%%%%% B 符号积分->数值计算 %%%%%%
% 一、原无限积分
disp('%%%%%% B 符号积分->数值计算 %%%%%%%%')
disp('一、原无限积分')
b=-log(10^-15);
f1=@(x) exp(-x)./x;
w1_k=@(u) vpa(int(f1,x,u,inf),20); % 符号积分-数值计算(原无限积分)
w1d=@(u) vpa(int(f1,x,u,b),20); % 符号积分-数值计算(原无限截断积分)
w1_k=w1_k(u)
w1d=w1d(u)