11楼: Originally posted by meranl at 2015-09-02 11:38:27
本来是要推荐mathematica的，想想算了，没必要折腾新软件

function plot_figure
clear all;clc
global lower_limit
n=2;
lspan=linspace(0,pi/2,100);
for i=1:length(lspan)-1   
lower_limit=lspan(i);

int1=  quadgk(@(x)  (tan(x/2)).^(2*n)./(sin(x)).^3,lower_limit,pi/2);

int2=  quadgk(@(x)  (n^2+1+ (sin(x)).^2 + 2*n*cos(x)   ).*sin(x).* (tan(x/2)).^(2*n)       ,lower_limit,pi/2);

k(i)=int1/int2;
end
[lspan(1:end-1); k]'

figure(1)
plot(lspan(1:end-1),k,'--o')

3楼: Originally posted by 月只蓝 at 2015-08-27 11:42:09
仅供参考，附图1是 k~φ0 的曲线关系，其中φ0是以弧度制表示的。
function plot_figure
clear all;clc
global lower_limit
n=2;
lspan=linspace(0,pi/2,100);
for i=1:length(lspan)-1   
lower_limit=lsp ...