your y2 function return Array for every single number as it's input.
syms x
y1 = linspace(0,1,10);
y2 = @(x) x.*y1;
y2(1)
according to documentation of integral function if fun has Array output, you should use option ArrayValued:
q = integral(y2,0,1,'ArrayValued',true)
