In [ ]:

    

using NMfE


    

In [5]:

    

Julia()


    

    
Out[5]:




true

In [6]:

    

setsymata(:xin, List(0.0,0.5,1.0))
setsymata(:yin, List(0.0,0.6,1.0))


    

    
Out[6]:




158

In [7]:

    

f = @sym Compile(Evaluate(LagrangePolynomial(xin,yin)));


    

In [8]:

    

println("\nf(0.75) = ", f(0.75))


    

    




f(0.75) = 0.8250000000000001

In [9]:

    

function lagrangefunc(x1,y1)
    setsymata(:xin, List(x1...))
    setsymata(:yin, List(y1...))
    symeval(parse("Compile(Evaluate(LagrangePolynomial(xin,yin)))"))
end


    

    
Out[9]:





lagrangefunc (generic function with 1 method)

In [11]:

    

f1 = lagrangefunc([0.0,0.5,1.0], [0.0,0.6,1.0])


    

    
Out[11]:





(::#7) (generic function with 1 method)

In [12]:

    

println("f1(0.75) = ", f1(0.75))


    

    




f1(0.75) = 0.8250000000000001

In [13]:

    

xl1 = [0, 1/3, 2/3, 1]
yl1 = [0, 0.3, 0.6, 1]
f2 = lagrangefunc(xl1, yl1)


    

    
Out[13]:





(::#9) (generic function with 1 method)

In [14]:

    

println("f2(0.75) = ", f1(0.75))


    

    




f2(0.75) = 0.8250000000000001

In [ ]: