Hi Altin,
Indefinite integration can be done in matlab using the int() function. For the above function the codes are
syms x;
int(x^3+9/x,x);
You integrate numerically if limits are provided using Trapezoidal method. http://in.mathworks.com/help/matlab/ref/trapz.html#bua4lsr
