4个回答
>> syms x y
epsilonr1=int(int((1/150)*(1/150)*exp((-i)*(0.0419*x+0.0419*y)),x,0,75),y,0, 75)
double(epsilonr1)
epsilonr1 =
(4*i*((10000*sin(1257/400))/419 - (20000*i*sin(1257/800)^2)/419))/(3771*exp((1257*i)/400)) - (80000*sin(1257/800)^2)/1580049 - (40000*i*sin(1257/400))/1580049
ans =
-0.1013 + 0.0001i
>>
>> int(int('y*exp(x)/(x+y)','y','2*y'),1,2) warning: explicit integral could not be found. ans = int(-(y*(ei(2*y) - ei(3*y)))/exp(y), y = 1..2)
ei是maple的函数，需要调用maple内核计算
图中是用maple计算的 i是虚数单位
要用MATLAB画一个表达式中含二重积分的曲线，但是不会写程序。看到您之前回答最简单的方法是用数值积分，但是上下限一改差别非常大，说明你这个函数的振荡
>> syms x y;
>> epsilonr1=int(int((1/150)*(1/150)*exp((-i)*(0.0419*x+0.0419*y)),x,0,75),y,0, 75)
epsilonr1 =
-40000/1580049*((-1)^(-1257/400/pi))^2+80000/1580049*(-1)^(-1257/400/pi)-40000/1580049
以上是我运行的结果，没有问题啊
你缺少M文件，未对int定义