1 files changed, 35 insertions, 0 deletions
diff --git a/macros/lar2rc.sci b/macros/lar2rc.sci
new file mode 100644
index 0000000..15fa4ee
--- /dev/null
+++ b/macros/lar2rc.sci
@@ -0,0 +1,35 @@
+function k=lar2rc(g)
+    
+//lar2rc convert log area ratios to reflection coefficients.
+// Calling Sequence
+// k = lar2rc(g)
+// Parameters
+// g: define log area ratios.
+// k: returns the reflection coefficients.
+// Examples
+//X = [7 6 5 8 3 6 8 7 5 2 4 7 4 3 2 5 4 9 5 3 5 7 3 9 4 1 2 0 5 4 8 6 4 6 5 3];
+// k = lar2rc(X)
+// or t=[2 5 6; 8 6 5; 8 9 4]
+// k = lar2rc(t)
+//
+// See also
+//
+// Author
+// Jitendra Singh
+//     
+//Modified to match MATLAB o/p when i/p is of type char and is a string by Debdeep Dey
+  if or(type(g)==10) then
+    [r,c]=size(g);
+    if r==1 & c==1 then
+        k=ones(1,length(g));
+    else
+        k=ones(size(g,1), size(g,2))
+    end
+    
+else
+    if ~isreal(g) then
+        error('Log area ratios must be real.')
+    end
+    k=-(tanh(-g/2));
+    end
+endfunction