1 files changed, 9 insertions, 0 deletions
diff --git a/gaussian_focus.m b/gaussian_focus.m
new file mode 100644
index 0000000..ff444ee
--- /dev/null
+++ b/gaussian_focus.m
@@ -0,0 +1,9 @@
+function [ w, s ] = gaussian_focus( w0, s0, f, lambda )
+%GAUSSIAN_FOCUS Summary of this function goes here
+%   Detailed explanation goes here
+
+zR = pi*w0^2/lambda;
+s = f*(1+(s0/f-1)/((s0/f-1)^2+(zR/f)^2));
+w = w0/sqrt((1-s0/f)^2+(zR/f)^2);
+end
+