% harjutus 4 - vedru venitamine

clear
x=[0 0.01 0.02 0.03 0.04 0.05 0.06]; % x(cm) -> x(m) tasakaalupunkti kaugusest
F=[0 5.1 10.2 16.1 25.4 40.2 61.9] % F(N)rakendatav jõud

%teeme võrgu, üle mille abil näeme graafikut
xvork=0:0.01:0.06; %NB! semikoolon lõppu!
Fkuup=interp1(x,F,xvork,'spline'); %NB semikoolon

%joonis
plot(x,F,'r*',xvork,Fkuup)
xlabel('kaugus x [m]')
ylabel('jõud F [N]')
legend('mõõdupunktid','kuupsplain')


vastustrapets=trapz(x,F) %trapetsvalemiga

%simpsoni omaga

f=@(t)interp1(x,F,t,'spline');, %defineerime nii matlebi jaoks kuupsplineõi
vastussimpson=quadv(f,0,0.06)

