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 ...

2楼: Originally posted by 1592203609 at 2015-08-26 23:01:16
确定可以用Origin作图？

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 ...

5楼: Originally posted by luo06009228 at 2015-08-27 18:59:42
不确定，所以今天装了MAtlab，不过楼下竟然画出来了，不知道是不是origin画的。...

7楼: Originally posted by 月只蓝 at 2015-08-27 19:43:30
MATLAB做的。
代码中已经给出了 k、φ0 的具体数值，输入到Origin画图也可。
Origin直接画图貌似很难实现。...

9楼: Originally posted by lfq6785676 at 2015-09-01 21:17:37
先用matlab算出数据，在用origin画图就行