# Generated by macrovsa_experiments.C do NOT edit
with(Statistics):with(plots):
X := Vector([100, 400, 1024, 2500, 4096, 10000, 1000000, 25000000], datatype=float):
Y := Vector([0.004445, 0.0157525, 0.0639475, 0.236397, 0.495985, 1.87314, 2089.46, 433039], datatype=float):
Xu := Vector([0.004445, 6.95311e-310, 0.0639475, 6.95311e-310, 0.495985, 1.87314, 6.95311e-310], datatype=float):
Yu := Vector([0.00095, 0.00079, 0.001645, 0.003435, 0.00507, 0.02295, 1.78942], datatype=float):
fit1 := Fit(a+b*d^c, X, Y, d, initialvalues = [a = 0.05, b = 0.00003, c = 1.35], output = [leastsquaresfunction, residualstandarddeviation]);
f1 := unapply(fit1[1], d):
fit2 := Fit(a+b*d^1.5, X, Y, d, output=[leastsquaresfunction, residualstandarddeviation], summarize=true);
f2 := unapply(fit2[1], d):
plotsetup(jpeg, plotoutput="../public/macrovsa_experiments/durations_interpolation.jpg", plotoptions="width=600,height=600"):
display(plot(f1(d), d=100..10000), pointplot(X, Y), title="binding computation time in msec");
fit3 := Fit(a+b*d, Xu, Yu, d, initialvalues = [a = 0.00, b = 0.00003], output = [leastsquaresfunction, residualstandarddeviation]);
f3 := unapply(fit3[1], d):
plotsetup(jpeg, plotoutput="../public/macrovsa_experiments/durations_interpolation_2.jpg", plotoptions="width=600,height=600"):
display(plot(f3(d), d=100..10000), pointplot(Xu, Yu), title="bundling computation time in msec");
quit: