Added files for Lagrange Polynomials

TestLagrangePolynomial.m

@@ -0,0 +1,22 @@
+%Test LagrangePolynomial function
+
+close all
+clear all
+
+number_points = 10;
+rng('default')  %Reset the pseudo-random number generator
+
+for(m = 2:number_points)
+    x_data = 10 * rand(1, m);
+    n = m - 1;
+    disp(['Degree = ' num2str(n)])
+    A = zeros(m);
+    for(k = 0:(m - 1))
+        A(k + 1, :) = LagrangePolynomial(x_data, k, x_data);
+    end
+    disp(A)
+    
+    if(~isequal(A, eye(m)))
+        error('Values do not match')
+    end
+end