k12 = 0.01;

delB = 0;

w = sqrt(k12^2+delB^2);

a0 = 1;



z = linspace(0,1000,1000);

a1 = a0*cos(w*z)-j*delB/w*sin(w*z);

a2 = -j*k12/w*sin(w*z);



subplot(311);

plot(z,abs(a1).^2); hold on;

plot(z,abs(a2).^2,'--r');

legend('P_1(z)','P_2(z)');

xlabel('Distance [{\mu}m]');

ylabel('Power [W]');

title('Case 1: {\Delta}{\beta} = 0');

