Boundary conditions are not satisfied after application of an interpolating function

I have a little question about boundaries conditions. I solved a differential equation and these are x-dependent functions. They are classified using two indices (i,j) like:

funcs[i][j]

For example if i=0,j=0 we have a particular function:

funcs

and the output is: Every function are defined in domain [0,a+Lin+i incrementoL+win+ j incrementow], thet depends on values of indices i and j, where 'a' is a constant, 'incrementoL' is a constant and 'incrementow' is a constant. At the start and at the end of domain the function is zero.

I can plot this function: In the domain defined, the function is zero at the boundary.

Now I nedd to crate an interpolation function, to vary the indices, so I create the following table:

matrixwave = Flatten[Table[{Lin + incrementoL i , win + incrementow j ,
funcs[i][j]}, {i, 0, passiL}, {j, 0, passiw}], 1]

where Lin,incrementoL,incrementow,passiL,passiw are real constants. Then I try to interpolate the functions:

intfuncs=Interpolation[matrixwave]

Then I try to plot several slices of the interpolation functions:

Plot[intfuncs[Lin + incrementoL,win + incrementow],{x,0,a+Lin+win+ incrementoL+incrementow}

and I obtain: and the function is not correct.It isn't zero a the boundaries. Any tips, helps ?