1 files changed, 42 insertions, 9 deletions

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