close all
clear all

g = 9.81;
lambda = 100;

x0 = 0:1:lambda;
x = 0:lambda/8:lambda;
z = 0:-lambda/16:-lambda/2;
z = z';

k = 2*pi/lambda;

w = sqrt(k*g);

T = 2*pi/w;

A0 = 4;
A = 1;

t = 0:T/100:3*T;

u = w*A*exp(k*z)*cos(k*x-w*t(1));
v = w*A*exp(k*z)*sin(k*x-w*t(1));

figure
set(gca,'FontSize',16);
p1 = plot(x0,A0*cos(k*x0-w*t(1)),x0,0,'b');
ylim([-lambda/2 2*A0]);
hold on;
q1 = quiver(x,z,u,v,.45,'MaxHeadSize',1.5);
set(q1,'LineWidth',2);
%axis equal;