From 0b992f0195b8d0abe9f74ac3d9294d1b1c14c44b Mon Sep 17 00:00:00 2001
From: nutricato
Date: Tue, 7 Jul 2009 14:41:59 +0000
Subject: 

---
 .../scilabcode/mainfunction.sci                    | 51 ++++++++++++++++++----
 1 file changed, 42 insertions(+), 9 deletions(-)

(limited to 'tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci')

diff --git a/tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci b/tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci
index afea60c6..5e2a7728 100644
--- a/tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci
+++ b/tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci
@@ -1,16 +1,49 @@
+//SCI2C: DEFAULT_PRECISION= DOUBLE
 function mainfunction()
 
 // ---------------------------------------
-// --- Initialization of the operands. ---
+// --- Initialization of the parameters. ---
 // ---------------------------------------
+
+
+// ---------------------------------------
+// --- Set number of point for quadrature .
+// --- interval for quadrature [a,b]
+// ---------------------------------------
+
+
+
+Np =10 // number of point for quadrature 
+a=1
+b=5
+ 
+
+// ---------------------------------------
+// --- Compute Legendre-Gauss-Lobatto nodes, weights 
+//. ---In standard [-1 1] interval
+// ---------------------------------------
+
+[x,w]=lglnodes(Np);
+disp('x');
+disp(x);
+disp('w');
+disp(w);
+// ---------------------------------------
+// --- Compute the nodes  
+//. ---In [a,b]  pick a t points in the 
+// --- desired interval
+// -------------------------------------
+
+t=((b-a)./2)*x+(b+a)./2;
+
+
+
+Iab=((b-a)./2)*mysum(w.*Fint(t),11);
+
+disp('Number of nodes for quadrature ');
+disp(Np)
+disp('Integral value in [a b] with Np points');
+disp(Iab)
 
-v1 = zeros(1,5);
-disp('v1, before initialization.')
-disp(v1)
-for cnt = 1:5
-   v1(cnt) = cnt;
-end
-disp('v1, after initialization.')
-disp(v1)
 
 endfunction
-- 
cgit