我有这个功能:
psi = @(x,y)(-(11*x^(6/5)*y^(1/5) - 19*x + 9)/(10*x^(1/10)*y^(1/10)))
解决x只能隐式完成:
R_sp = @(y)solve(psi(x,y), x);
我需要用R_sp作为下限评估以下双积分:
int(s * int(s*t, t, R_sp(s), 1), s, .3, .4)
评估时,matlab说:MuPAD命令出错:需要列表或条目表。 [(DOM ::矩阵(DOM :: ExpressionField()))::创建]
R_sp的值在[.3,.4]中有明确定义,可以很容易地计算出来:
R_sp(.3) = 0.845993307852830503940553192774
R_sp(.4) = 0.89641082924669996315283636152907
所以我假设matlab的int函数不能处理该隐式函数中的连续值,对吧?当然在这种情况下,我可以使用trapz函数通过计算for循环中被积函数的离散值对来近似上面的积分。但除了使用陷阱近似(使用连续值)之外,还有其他任何方法可以在matlab中评估这个积分吗?