lab2

Tn=0;
Tv=2*pi/3;
t=Tn:0.005:Tv;
x1=t.*t+2*cos(t)+sin(3*t);
figure(1);
plot(t, x1, 'k');
hold on;
grid on;
t0=Tn;
dt=pi/6;
for j=1:1:5
f(j)=dt*(j-1).*dt*(j-1)+2*cos(dt*(j-1))+sin(3*dt*(j-1));
plot(pi.*(j-1)/6, f(j), 'gh');
end
df11=f(2)-f(1);
df12=f(3)-f(2);
df13=f(4)-f(3);
df21=df12-df11;
df22=df13-df12;
df3=df22-df21;
p=