%
% Aproximacao cubica dos minimos quadrados e interpolacao com splines
%
clear;
clc;
% dados
t=0:1:10;
y=sin(t);
% valores onde se pretende calcular o polinomio aproximado
t1=0:0.1:10;
% opoees
r=input('Aproximacao global dos minimos quadrados (1) ou segmentada com splines (2) ? ');
switch r
case 1
   coef=polyfit(t,y,3);
   y1=polyval(coef,t1);
   plot(t,y,'-o',t1,y1,'r'), title('Aproximacao cubica global'), legend('dados','P_3(t)')
case 2
   y1=spline(t,y,t1);
   plot(t,y,'-o',t1,y1,'r'), title('Interpolacao cubica segmentada'), legend('dados','S_3(t)')
otherwise
   disp('Erro')
end
