%??Procedure?for?calculating?the?spatial?impulse?response
%??for?an?aperture.
%
%??Calling:??[h?start_time]?=?calc_h（Thpoints）;?
%
%??Parameters:??Th?????-?Pointer?to?the?transducer?aperture.??
%???????????????points?-?Field?points.?Vector?with?three?columns?（xyz）?
%????????????????????????and?one?row?for?each?field?point.?
%
%??Return:??????h??????????-?Spatial?impulse?response?in?m/s.
%???????????????start_time?-?The?time?for?the?first?sample?in?h.
%
%??Version?1.01?October?4?1996?by?Joergen?Arendt?Jensen

function?[h?start_time]?=?calc_h?（Thpoints）

%??Check?the?point?array

??[mn]=size（points）;
??if?（n?~=?3）
????error?（‘Points?array?must?have?three?columns‘）;
????end

%??Call?the?C-part?of?the?program?to?show?aperture

??[h?start_time]?=?Mat_field?（4001Thpoints）;