Matlab求的数值积分,为什么结果全是NaN?(Matlab2010b)
发布网友
发布时间：2024-10-23 22:20
共1个回答
热心网友
时间：2024-11-02 21:24
clear
%x不能等于0，因为是分母，eps是matlab里最小的正数，步长可以取的再小点，更准确
x=eps:0.01:1;
c=1:11;
P=zeros(1,11);%最好先给定P的大小，可以加快执行速度，这是一个好的习惯
for a=1:11
y=(x.^(-3)+c(a)^3-1).^(-1/3);%此处c应该只取一个值，故用c(a)
h=(y./x.^2-y.^7./x.^8)*2;
P(a)=trapz(x,h);
end
P
