![]() |
请教符号积分时遇到的问题
想解一个二重积分,看网上有人说可以先做符号积分int,再做数值积分quadv。可是我int那里就出问题了,不知道怎么解决。
clear; syms t w L=0.08; a=0.1;b=0.08;c=-0.15;uz=0.2;vz=0.2;wz=0.1; f1=[(-3/20+w)/((2/25-t)*(2/25-conj(t))+(-3/20+w)*(-3/20+conj(w)))^(3/2)+(3/20-w)/(1/25+(2/25-t)*(2/25-conj(t))+(-3/20+w)*(-3/20+conj(w)))^(3/2)+(-3/20-w)/((2/25-t)*(2/25-conj(t))+(-3/20-w)*(-3/20-conj(w)))^(3/2)+(3/20+w)/(1/25+(2/25-t)*(2/25-conj(t))+(-3/20-w)*(-3/20-conj(w)))^(3/2), 0, 1/5/(1/25+(2/25-t)*(2/25-conj(t))+(-3/20+w)*(-3/20+conj(w)))^(3/2)+1/5/(1/25+(2/25-t)*(2/25-conj(t))+(-3/20-w)*(-3/20-conj(w)))^(3/2)]; f2=int(f1,t,-0.1,0.1); %%%%%出错的地方 f2=inline(f2/L); result=quadv(f2,w,-wz,-wz-L); f1有点复杂,是个行向量,带有两个变量t和w;int是想对t求积分,quadv想对w求积分。 ??? Attempted to access e(1); index out of bounds because numel(e)=0. Error in ==> sym.sym>char2sym at 518 sk = x(s(k):e(k)); Error in ==> sym.sym at 95 S = char2sym(x); Error in ==> sym.maple at 92 result = sym(result,'keepijalias'); Error in ==> sym.int at 51 r = reshape(maple('map','int',f(:),[x.s '=(' a.s ')..(' b.s ')']),size(f)); Error in ==> tmpy at 10 f2=int(f1,t,-0.1,0.1); 错误信息如上,不知道怎么解决。请各位大虾帮帮忙,谢谢! |
所有时间均为北京时间。现在的时间是 05:32。 |
Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.