beam with 2 polarization in ecample

1 files changed, 18 insertions, 14 deletions

diff --git a/example_beam_trace.m b/example_beam_trace.m
index 95b5ff3..de980b7 100644
--- a/example_beam_trace.m
+++ b/example_beam_trace.m
@@ -1,8 +1,12 @@
-beam.k=[-1,0];
-beam.origin=[2,0.5];
-beam.face=NA;
-beam.intensity = 1;
-beam.polarization = 1; % 1 for s  and 2 for p
+beam1.k=[-1,0];
+beam1.origin=[2,0.5];
+beam1.face=NA;
+beam1.intensity = 1;
+beam1.polarization = 1; % 1 for s  and 2 for p
+
+beam2=beam1;
+beam1.polarization = 2; % 1 for s  and 2 for p
+
 
 face1.vertex1=[-1,0];
 face1.vertex2=[1,0];
@@ -24,26 +28,26 @@ faces={face1,face2,face3};
 
 border1.vertex1=[-2,-2];
 border1.vertex2=[2,-2];
-border1.n_right = 1.;
-border1.n_left= 1;
+border1.n_right = [1, 1];
+border1.n_left= [1, 1];
 
 border2.vertex1=[2,-2];
 border2.vertex2=[2,2];
-border2.n_right = 1.;
-border2.n_left= 1;
+border2.n_right = [1, 1];
+border2.n_left= [1, 1];
 
 border3.vertex1=[2,2];
 border3.vertex2=[-2,2];
-border3.n_right = 1.;
-border3.n_left= 1;
+border3.n_right = [1, 1];
+border3.n_left= [1, 1];
 
 border4.vertex1=[-2,2];
 border4.vertex2=[-2,-2];
-border4.n_right = 1.;
-border4.n_left= 1;
+border4.n_right = [1, 1];
+border4.n_left= [1, 1];
 
 borders={border1,border2, border3, border4};
-beams={beam};
+beams={beam1,beam2};
 
 
 img=zeros(1000,1000);