function[key]=wave（）
beta1=[input（‘請輸入遭遇角：‘）];
V1=[input（‘請輸入航速（節）：‘）];
beta=beta1*pi/180;
V=V1*1852/3600;
g=9.8;
A=0.0081;
B=10;
N=25;

for?t=1:1:200
????a1=0;
????for?w=0.4:0.05:1.6
????????kt=1-0.0833*w+0.0032*w^2-0.0000044*w^6;
????????kb=1+0.0049*