Integration of scattered data

Hello,

I'm having the following problem: I'm using TriScatteredInterp to interpolate discrete points.

   vq=TriScatteredInterp([x(:),y(:),z(:)],p_abs(:),'natural')

i.e. the value at (x1,y1,z1) can be evaulated with vq(x1,y1,z1).

I would like to integrate this value over the thickness, e.g. at x=0, y=0.

   f=@(z)vq(0,0,z); 
   quad(f,0,1)

I get the following error:

   Error using TriScatteredInterp/subsref 
   The input data has inconsistent size. 

   Error in @(z)vq(0,0,z) 

   Error in quad (line 76) 
   y = f(x, varargin{:}); 

   Error in fit_press_test_3d (line 30) 
   quad(f,0,1)

Integration with trapz gives the same error.

Answer 1

[ChristianW]

vq(scalar,scalar,vector) doesnt work.

f=@(z)vq(zeros(size(z)),zeros(size(z)),z);

[PJ]

Thank you !