3 files changed, 29 insertions, 0 deletions
diff --git a/tests/unit_tests/test_jacobi/scilabcode/main.sci b/tests/unit_tests/test_jacobi/scilabcode/main.sci
new file mode 100644
index 0000000..ee229eb
--- /dev/null
+++ b/tests/unit_tests/test_jacobi/scilabcode/main.sci
@@ -0,0 +1,18 @@
+// Demo function to test %sn and %k functions
+// %sn is the Jacobi's elliptic function
+// %k is the Jacobi's complete integral
+function main()
+	u=[1+2*%i 2+3*%i 3+1*%i];  // the corresponding complex vector
+        k=0.6;                     // parameter of elliptic integral
+        [out]=%sn(u,k);
+        disp(out);
+	m=[0.1 0.2 0.3];	   // parametrt of Jacobi complete integral
+        s=%k(m);
+        disp(s);
+endfunction
+
+//Output
+//for %sn
+// 1.5641568 - 0.2037920i    1.0580372 + 0.0100213i    1.1386855 - 0.4631921i  
+// for %k
+// 1.6124413    1.6596236    1.7138894 
diff --git a/tests/unit_tests/test_jacobi/scilabcode/modktest.sci b/tests/unit_tests/test_jacobi/scilabcode/modktest.sci
new file mode 100644
index 0000000..613cd33
--- /dev/null
+++ b/tests/unit_tests/test_jacobi/scilabcode/modktest.sci
@@ -0,0 +1,5 @@
+function modktest()
+    m=[0.1 0.2 0.3];
+    s=%k(m);
+    disp(s);
+endfunction
diff --git a/tests/unit_tests/test_jacobi/scilabcode/modsntest.sci b/tests/unit_tests/test_jacobi/scilabcode/modsntest.sci
new file mode 100644
index 0000000..dc5d435
--- /dev/null
+++ b/tests/unit_tests/test_jacobi/scilabcode/modsntest.sci
@@ -0,0 +1,6 @@
+function modsntest()
+    u=[1+2*%i 2+3*%i 3+1*%i];  // the corresponding complex vector
+    k=0.6;                     // parameter of elliptic integral
+    [out]=%sn(u,k);
+    disp(out);
+endfunction