Hello, users of forum

i need of help with o scilab, i am doing a computational physics project related a musical instruments

For that, i have what solve a differential partial equation of fifth order with a term partial of fifth order.

but it displays the following error:

Sorry,I still do not speak English perfectly

xilofone(0.1,0.1,1)
!--error 144
Undefined operation for data operands.
Check or set the function %s_a_l for overloading.
at line 5 of function xilofone called by :
xilofone(0.1,0.1,1)

below is the code of the program:

function[V]=G(x)
V=cos((2/10)*%pi-3)*x;
endfunction

function[o,f]=xilofone(ht,hx,L)
X=zeros(L+5,L+5);
for i=3:L+5
for j=3:L+5
X(i,j)=X(i,j)-(1\ht)*(X(i+1,j)-X(i-1,j))-((1\ht)^2)*(X(i+1,j)-2*X(i,j)+X(i-1,j))+((1\hx^4))(X(i+2,j)-4*X(i+1,j)+6*X(i,j)-4*X(i-1,j)+X(i-2,j))-(X(i+2,j)-4*X(i+1,j)+6*X(i,j)-4*X(i-1,j)+X(i-2,j)-(X(i+2,j-1)-4*X(i+1,j-1)+6*X(i,j-1)-4*X(i-1,j-1)+X(i-2,j-1))\ht*(hx^4));
end
end
o=matrix((1:L),L,1)
u=linsolve(X,o);
plot(u,o)
endfunction