Question

怎么用matlab计算这个二重积分

Answer 1

题主给出的二重积分，用matlab的int()函数无法得到其解析值。如l=0.5，R=0.39，则可以通过下列代码，求得其数值解

l=0.5,R=0.39

syms x y

x0=sqrt(R^2-l^2/2);

x1=l/2;

y0=sqrt(R^2-l^2*x^2/4);

y1=l/2;

fun=R./sqrt(R^2-x.^2-y.^2);

I=int(int(fun,y,y0,y1),x,x0,x1);

I=vpa(I)

追问

但是文章里说是可以解的啊 l等于1的情况下R在这个区间变化求得的积分值变化如图第三段