Matlab - Problema punto mas cercanos

clear all

clc

close all

syms x

n=10;

f=x^2+1;

X=linspace(0,2,n);

F=subs(f, x, X);

plot(X,F,'-*')

set(gca,'Xtick', X);

set(gca,'Ytick', double(F));

grid on

axis([X(1)-0.5 X(end)+0.5 min(double(F))-0.5 max(double(F))+0.5 ])