I am a starter for Phased array calculation, having this toolbox is much convenient and useful for my study.

Recently, I want to import a NEC output pattern for the calculation and I tried the example named 'Ex3b.m' in the 'Examples' folder. And I found it will raise a error while executing the statements as:

Then error message from Matlab commander as follows:
Calculating Directivity : d(Th)=10.00 d(Phi)=10.00
| |
...??? Error using ==> interp at 37
Input parameters are out of range.

Error in ==> SUMCODE at 33
case 11,[EleAmp,CPflag]=interp(thloc,philoc);

Error in ==> fieldsum at 100
[EleAmp,CPflag]=sumcode(eltype,thloc,philoc); % Sumcode calls the relevant
element model 'eltype'

Error in ==> calc_directivity at 67
Emultiple=fieldsum(range_config,theta,phi); % E(theta,phi)

Error in ==> Test_Ex3b at 1339
calc_directivity(10,10); % Calc directivity using 10deg
theta and phi steps

------------------------------------

If I changed 'interp' to 'user1' or 'patchr' something in function 'place element', the program running successfully without errors.
Therefore, I doubted that there are some mistakes using "interp" function when calculating the directivity. My matlab version is Matlab r2009a (V7.8). Could you comment it and give me some advices about how to solve the problem?

Thank you so much for your kindest attention and help.