ジュリア：エルロとは何ですか？ InexactError（） つつ：ラグランジュ補間は

Question

function lg(X,Y,Xint) 

    n = length(X) 
    L = ones(1:n) 

    for k = collect(1:n) 

    L[k] = 1 
    for c = collect(1:n)    
     if c!=k 
      L[k] = (L[k]*(Xint - X[c]))/(X[k] - X[c]) 
     end 
    end   
    end 
    return sum(Y.*L) 
end 

=========================

実行さ

とLoadError機能しますローディング[76]で、ライン1から始まる式で

In [74]：11？範囲の

Answer 1

onesがInt64の配列を作成：

julia> o = ones(1:3) 
3-element Array{Int64,1}: 
1 
1 
1 

julia> o[1] = 3.5 
ERROR: InexactError() 
in setindex!(::Array{Int64,1}, ::Float64, ::Int64) at ./array.jl:339 
in eval(::Module, ::Any) at ./boot.jl:226 

あなたはInt64の配列へのfloat64を割り当てることはできませんが（このエラーが発生します）。

あなただけのfloat64の配列を取得するためにones(n)を使用したい：

julia> ones(3) 
3-element Array{Float64,1}: 
1.0 
1.0 
1.0 

---

サイドノート：あなたは範囲を反復処理する前collectする必要はありません。

for k = collect(1:n) 

代わりに、範囲を反復するだけです：

for k = 1:n