function  [integral_outvalue, fun_outvalue]=integral_iota(v)    
%return integral as first argument
integral_outvalue=quad(@iota,v,5);
%return derivative of the integral, i.e. the function itself as the second
%argument
fun_outvalue=iota(v);
end