syms rho rho0
integrand = rho^3 / ((rho^2 - rho0^2) / (1 - rho0^2) - rho^6)^(1/2);
zeta = int(integrand, rho, rho, 1);
disp(zeta);
rho0 값에 따라 결과가 바뀔텐데, 이 값이 0.1인 경우는 아래와 같이 계산시킬 수 있습니다.
subs(zeta, rho0, 0.1)

Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!