matlab色散方程

pi=3.1415926,g=9.8;

T=input('输入周期(s)T:\n');

h=input('输入水深(m)h:\n');

x0=1.0e-8;

x=(4*pi*pi*h)/(g*T*T*tanh(x0));

while abs(x-x0)>1.0e-8;

x0=x

x=(4*pi*pi*h)/(g*T*T*tanh(x0));

end

L=2*pi*h/x;

k=2*pi/L;

c=L/T;

str1=strcat('你输入的周期：',num2str(T), ' s ');

str2=strcat('你输入的水深：',num2str(h), ' m ')

str3=strcat('计算所得的波长：',num2str(L), ' m ') str4=strcat('计算所得的波速：',num2str(c), ' m /s') str5=strcat('计算所得的波数：',num2str(k), ' 1 /m') disp(str1), disp(str2) ,disp(str3), disp(str4) ,disp(str5)