[d_spline,gam,h,ssM,sM,dM] = smu_spline(time,d_data,0.000075);


d_res3=d_data-d_spline;


figure
plot(time,d_data);
%axis([0.0 0.5 -0.00006 0.00006]);
axis tight;
hold on;
plot(time,d_spline,'g--');
hold off;
legend('experimental data','spline model');
xlabel('Time','FontSize', 18);
ylabel('Displacement','FontSize',18);


figure
clf;
plot(time,d_res3, 'g.');
axis tight;
%axis([0.0 0.5 -0.00006 0.00006]);
ylabel('Residual', 'FontSize', 18);
xlabel('Time', 'FontSize', 18);
grid on;

figure;
clf;
qqplot(d_res2);