f1='y*(x+y)/4'
f2=int(f1,'x','1','y*y')
f3=int(f2,'y','1','z')
p1=sym2poly(f3)
p2=conv(p1,[1 0])
f4=poly2sym(p2,'z')
f5=int(f4,'z','1','3')
numeric(f5)
f1 =
y*(x+y)/4
f2 =
1/8*y*(y^4-1)+1/4*y^2*(y^2-1)
f3 =
1/48*z^6+3/40+1/20*z^5-1/12*z^3-1/16*z^2
p1 =
0.0208 0.0500 0 -0.0833 -0.0625 0 0.0750
p2 =
0.0208 0.0500 0 -0.0833 -0.0625 0 0.0750 0
f4 =
1/48*z^7+3/40*z+1/20*z^6-1/12*z^4-1/16*z^3
f5 =
194/7
ans =
27.7143
|