Use of Evalin in matlab

for ii=1:1:5
  xa(ii)=(2*ii+1)/4/pi*int(int(sin(thea)*dot([sin(thea)*cos(phia) sin(thea)*sin(phia) cos(thea)],vam)*evalin(symengine,'orthpoly::legendre(ii, cos(thea))'),thea,0,pi),phia,0,2*pi)
end

Here come the error Warning: Explicit integral could not be found. I think that is the function 'orthpoly::legendre' contains 'ii' which is the variable in 'for'. Is there any way to deal with 'ii' in evalin statement. Thanks

Answer 1

[Azzi Abdelmalek]

evalin(symengine,sprintf('orthpoly::legendre(%d, cos(thea))',ii),thea,0,pi),phia,0,2*pi)

[Xin]

Sorry to bother you. It doesn't work. Any other method? Thanks

[Azzi Abdelmalek]

What is your error messages

[Xin]

Error using mupadengine/evalin

Too many input arguments.

Thanks

[Walter Roberson]

 feval(symengine, 'orthpoly::legendre', ii, cos(thea))