你好!
clc;clear all
format long
%要么
Y1=quad2d(@(x,y)(y.*(y-0.042))./(y.*(2+pi*sin(x))-0.084),pi/60,(4.2724/180)*pi,@(x)(0.084./(2-sin(x))),@(x)(0.0905./(2+sin(x))))
%或者
Y2=integral(@(x) arrayfun(@(x) integral(@(y)(y.*(y-0.042))./(y.*(2+pi*sin(x))-0.084),(0.084./(2-sin(x))),(0.0905./(2+sin(x)))),x),pi/60,(4.2724/180)*pi)
试一试吧!!!
Y1 =
6.636967271338606e-08
Y2 =
6.636967272022144e-08