beam.k=[1,-.5];
beam.origin=[0,0];
beam.face=NA;
beam.intensity=1;

face.vertex1=[1.1,-1];
face.vertex2=[2.7,1];
face.n_right = 1;
face.n_left= 2.5;

face2.vertex1=[-.1,-3];
face2.vertex2=[-.1,3];
face2.n_right = 1;
face2.n_left= 1;

face3.vertex1=[2.1,-3];
face3.vertex2=[2.1,3];
face3.n_right = 1.5;
face3.n_left= 1;

faces={face,face2,face3};

[is_face_hit, hit_position, hit_distance, new_beams] =  face_beam_interaction(beam, faces);
beam_reflected=new_beams{1};
beam_refracted=new_beams{2};


[is_face_hit, hit_position, hit_distance, new_beams] =  face_beam_interaction(beam_reflected, faces);

[is_face_hit, hit_position, hit_distance, new_beams] =  face_beam_interaction(beam_refracted, faces);