function mu = myMoment(n,name) mu = zeros(n+1,1); for i=0:n if strcmp(name,'legendre') if mod(i,2) == 0 mu(i+1) = 2/(i+1); end elseif strcmp(name,'log') mu(i+1) = 1/(i+1)^2; elseif strcmp(name,'tscheb') if mod(i,2) == 0 if i == 0 mu(i+1) = pi; else v1 = 1:2:i-1; v2 = 2:2:i; mu(i+1) = pi*prod(v1)/prod(v2); end end elseif strcmp(name,'laguerre') mu(i+1) = factorial(i); elseif strcmp(name,'hermite') if mod(i,2) == 0 mu(i+1) = gamma((i+1)/2); end end end