Multidimensional integration with Normal inverse cumulative distribution function

I am trying to perform multidimensional integration, using the following algorithm:

q = @(x,y,m,n)abs(sin(norminv(n,0,1))+sin(norminv(m,0,1))+sin(norminv(y,0,1))).*exp(-norminv(x,0,1));
q1 = int(q,  @n, 0,1);
q2 = int(q1, @m, 0,1);
q3 = int(q2, @y, 0,1);
q4 = int(q3, @x, 0,1);
vpa(q4)

but this gives the following error:

Undefined function 'int' for input arguments of type 'function_handle'.

Error in untitled2 (line 8)
q1 = int(q,  @n, 0,1);

I am quite new in MATLAB and I think that the problem can be in non-symbolic usage of variables.

Can somebody please help me understand and fix my mistake?