公告:1)网站程序升级:Q2A升级到1.8.6,Wordpress升级到5.7.2
2)修复了头像加载慢与提交问题反应慢等问题
2021-06-16

欢迎来到 Mathematica 问答社区

提问时请贴上文本代码

语法高亮:在编辑器中点击

被禁止的话题:广告破解

请阅读:《提问的智慧》

备用域名:mma.ooo

支持LaTex数学公式
行内公式标识符:\$ 或“$\backslash ($”+“$\backslash )$”,
行间公式标识符:\$\$ 或 “$\backslash [$”+“$\backslash ]$”

社区建议QQ群:365716997

分类

0 投票
287 浏览
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,该怎么做?

分类:列表操作 | 用户: keanhy (361 分)

1个回答

+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}]

用户: 随堂测验 (606 分)
采纳于 用户:keanhy
没有好好理解Table基本用法,再次感谢耐心解答!
...