[yuanwanwu]

phi=0:2*pi/90:2*pi;
f=@(phi)sqrt(1-K1.^2*sin(phi*180/pi)*sin(phi*180/pi));
Q1=quad(f,0,pi/2);



为什么程序运行总说这里错了，搞不懂