Hello:
    I  want to get the variational value of omega when KR[N]>0.2, but  the omega of my procedure remain the same,  how do I change my program to get the variational value of omega when KR[N]>0.2.

Thank you !

My procedure:
restart:
> Digits:=30;
>
> h0:=0.156;
> d:=0.32*h0;
> l:=2;
> h1:=h0-d;
> h2:=h0+d;
> h3:=0.6*h0;
> g:=9.8;
> d1:=1;

Dear Sir,  I want to know why the procedure  appear the error : (in U1) unable to execute add.

My procedure:

restart:
>
> h0:=0.156;
> d:=0.32*h0;
> l:=1;
> h1:=h0-d;
> h2:=h0+d;
> h5:=evalf(2/3*h0);
> h7:=h0;
> s1:=evalf(2/3);
> s2:=evalf(2/3);
> g:=9.8;
> d1:=1;
> Term:=10;
> Num:=5:
> n:=1:
>
>
> for N from 1 to Num do
> k0:=evalf(0.5*Pi+2*(N-1)*Pi/(Num-1)):...

    Why the zeta1 does not plug in numerical calculation.
 My program:

restart:
>
> h0:=0.156;
> d:=0.32*h0;
> l:=1;
> h1:=h0-d;
> h2:=h0+d;
> h5:=evalf(2/3*h0);
> h7:=h0;
> s1:=evalf(2/3);
> s2:=evalf(2/3);
> g:=9.8;
> d1:=1;
> Term:=10;
> Num:=5:
> n:=1:
>
>
> for N from 1 to Num do
> k0:=evalf(0.5*Pi+2*(N-1)*Pi/(Num-1)):

restart:
> Digits:=30:
> with(linalg):
>
> h[0]:=0.156;:
> d:=0.32*h[0]:
> l:=4:
> h[1]:=h[0]-d:
> h[2]:=h[0]+d:
> g:=9.8:
> Term:=12:
> Num:=3:
> n:=4:
>
> for N from 1 to Num do
> lambda:=2*n*Pi/l:## N1 wei sha ba tiao shu
> epsilon:=evalf(0.5+2*(N-1)/(Num-1)):
> k[0]:=evalf(0.5*Pi+2*(N-1)*Pi/(Num-1)):
> tau[0]:=evalf(k[0]*h[0]):
> omega:=evalf((g*k[0]*tanh(k[0]*h[0]))^(1/2)):