lab 4 / numer_diff

function numer_diff(x,y)
p=x(1):0.1:x(end);
f7=interp1(x,y,'spline');
f=diff(f7);
f=[f f(end)];
h=p(2)-p(1);
f=f/h;
fg=gradient(f7,h);
f8=@(p)interp1(x,y,p,'spline');
for i=1:length(p)
Integ(i)=integral(f8,p(1),p(1));
end
figure;
yyaxis left
plot(p,f,p,fg)
xlabel('X');
ylabel('Df(x)');
yyaxis right
plot(p,Integ)
title('Numerical gradient and intergration');
ylabel('Intergral(x)');
legend({'diff', 'gradient','integral'},'Location','best');
end