# Returns integral of a polynomial function given by array 'p' within layer #layer

function [ val ] = integral( layer, p )

	x = polyint(p);

#layer
#x

	val = dz(layer)*(polyval(x,1.0)-polyval(x,0.0));

endfunction