In [2]:

    

include("./lagrange_spec.jl")
include("../src/HPFEM.jl")


    

    
Out[2]:





HPFEM

In [3]:

    

using PyPlot


    

In [4]:

    

lagrange_spec(12,13,2,fun)


    

    
Out[4]:





0.9924242424242424

In [6]:

    

fun(x) = sin(2*pi*x)
Marray = 3:50
j = 2

erros = Float64[lagrange_spec(M,M,1,fun) for M in Marray]
Marray2 = Float64[Marray;]
loglog(Marray2,erros,"r",label="Q = M")

for j in 2
    erros2  = Float64[lagrange_spec(M,M+j,1,fun) for M in Marray]
    loglog(Marray2,erros2,label = "Q = M+$j")
end    
legend()
title("loglog")


    

    













    
Out[6]:





PyObject <matplotlib.text.Text object at 0x7fdc2121a810>

In [ ]: