script correction

hi, would somebody pleas help me to correct this script?

clear
clc;
BHv= [1,2,3,4,5,6,7,8];
fiv= [0,10,20,30];
h = [ 1 1.02  1.11  1.21  1.3   1.4  1.59  1.78; ...    
        1    1.11  1.35  1.62  1.95  2.33 3.34  4.77; ... 
        1.01 1.39  2.12  3.29  5.17  8.29  22    61; ...
        1.13  2.5  6.36  17.5   50   150  1400  14800];

[fim BHm] = meshgrid(fiv,BHv);
[XI,YI] = meshgrid(linspace(0,30),linspace(1,8));
F = TriScatteredInterp(fim(:),BHm(:),h(:));
figure,surf(F(XI,YI))

ylabel('y-axis,(B/H)')

xlabel('x-axis,phi')

zlabel('depth factor,hc')

title('depth factor(hc) measurements,mandel & salencon 1969')

[milad babaei]

after run the code.it doesnt show exact data limitation on graph,have a problem introducing my data,