clear;clcf=@(x,a)det([x^2+a sqrt(x^2+4*a);exp(x) 2*x+a^2])y=@(x)quadl(@(a)arrayfun(@(a)f(x,a),a),0,2)fplot(@(x)arrayfun(y,x),[0,5])