using Polynomials
x=df2[:RozmiarMacierzy]
y=df2[:avg_Naive_t]

y_Naive_std=df2[:std_Naive_t]
yNaive=df2[:avg_Naive_t]

y_Better_std=df2[:std_Better_t]
yBetter=df2[:avg_Better_t]

y_BLAS_std=df2[:std_BLAS_t]
yBLAS=df2[:avg_BLAS_t]

fitNaive=polyfit(x,yNaive,3)
fitBetter=polyfit(x,yBetter,3)
fitGsl=polyfit(x,yBLAS,3)
xD=100:100:500

scatter(x,y,yerr=y_Naive_std)
plot!(xD, polyval(fitNaive, xD))
plot!(xD, polyval(fitBetter,xD))
plot!(xD,polyval(fitGsl,xD))