怎么用 matlab画出这个公式的图，求编程过程

写了个m file，运行得到你的图像，不过有点小小的bug： clear figure1=figure('color',[1 1 1]); syms x; k=0:0.2:10; n=zeros(1,length(k)); for i=2:length(k) y=k(i)/2*int(1/((x+k(i)^2)^(3/2)*(x+1)),x,0,inf); n(i-1)=vpa(y,10); end plot(k,n