I think I know what the problem is. For some reason DarkBasic doesn't automatically calculate the normals for the matrices, so they don't render with any shadows/highlights, just flat colour. Here's some code that manually calculates the normal data for your matrix (You need only run this code once after making changes to your matrix):
for z=1 to 32
for x=1 to 32
rem Get matrix heights
h8#=get matrix height(currentmatrix,x,z-1)
h4#=get matrix height(currentmatrix,x-1,z)
h#=get matrix height(currentmatrix,x,z)
h2#=get matrix height(currentmatrix,x,z)
rem Calculate projected angle X using heights
x1#=(x-1)*25.0 : y1#=h#
x2#=(x+0)*25.0 : y2#=h4#
dx#=x2#-x1#
dy#=y2#-y1#
ax#=atanfull(dx#,dy#)
ax#=wrapvalue(90-ax#)
rem Calculate projected angle Z using heights
z1#=(z-1)*25.0 : y1#=h2#
z2#=(z+0)*25.0 : y2#=h8#
dz#=z2#-z1#
dy#=y2#-y1#
az#=atanfull(dz#,dy#)
az#=wrapvalue(90-az#)
rem Make normal from projected angle
nx#=sin(ax#)
ny#=cos(ax#)
nz#=sin(az#)
rem Setting matrix normal for smoothness
set matrix normal currentmatrix,x,z,nx#,ny#,nz#
next x
next z
update matrix currentmatrix
Change the '32' in the 'for x' and 'for z' loops for however many x-segments and z-segments your matrix has.