数值计算时，怎么在某一区间取更多的点？我下边的程序能通过Table实现吗？

Question

Clear["Global`*"]
W[x_] := a1 (x^(-4) + 2 x^2 - 3)
p2[x_] := 
 Integrate[
  D[W[x], x]/(1 - x^3)/
   a1, {x, ((\[Delta] + 1)^3*(xb^3 - 1) + 1)^(1/3), xb}]
yy = Table[
   NIntegrate[
    D[W[x], x]/(1 - x^3)/
     a1, {x, ((\[Delta] + 1)^3*(xb^3 - 1) + 1)^(1/3), xb}], {\[Delta],
     0.5, 1, 0.5}, {xb, 1, 4, 0.1}];
xx = Range[1, 4, 0.1];
ListLinePlot[Transpose[{xx, yy[[2]]}], Mesh -> All, 
 AxesLabel -> {xb, p2}]

图形结果如下

程序中1至4步长都是0.1，但我想在1至1.5取步长0.01,1.5之后取步长0.01，该怎么做？

Answer 1

可以的，可以用Table的第五种用法，Table[expr,{i,{i1,i2,...}}].或者事先定义i的列表，如

j=Range[5];Table[i^2,{i,j}]

对于你这个，

xblist = Join[Range[1, 1.5, 0.01], Range[1.6, 4.0, 0.1]];

ylist = Table[NIntegrate[D[W[x], x]/(1 - x^3)/a1, {x, 
         ((\[Delta] + 1)^3*(xb^3 - 1) + 1)^(1/3), xb}], {\[Delta],0.5, 1, 0.5},
         {xb, xblist}];
ListLinePlot[Transpose[{xblist, ylist[[2]]}], Mesh -> All,AxesLabel -> {xb, p2}]

Comments

没有好好理解Table基本用法，再次感谢耐心解答！