Script:
- Codice:
g=@(t)(t.*(sin(t)).^2.*exp(-t));
m=10;
a=-2;
b=1;
F=@(x)(x(end)-(x(:,end)+1)/6/m.*( g(x(:,1))+2*sum(g(x(:,3:2:2*m)),2)+4*sum(g(x(:,2:2:2*m)),2)+g(x(:,2*m+1))) );
k=0;
t=linspace(a,b,100);
for tt=t
k=k+1;
x(k,:)=linspace(-1,tt,2*m+1);
end
figure(1)
hold on
plot(t,F(x),'k')
pause
n=3;
A=0;
x = linspace(a,b,m+1);
for j=1:m
x_in=linspace(x(j),x(j+1),n+1)';
for k=1:n+1
xF=linspace(-1,x_in(k),2*m+1);
y(k)=F(xF);
end
a = get_polyn(x_in,y');
t = linspace(x(j),x(j+1),100);
yp= polyval(a,t);
plot(t,yp,'r');
A=A+Cavalieri_Simpson_polinomi(x(j),x(j+1),m,a);
pause
end
title('Grafico')
legend('F','polinomio composito')
xlabel('asse x')
ylabel('asse y')
però il polinomio interpolante non interpola per niente $F$ nei nodi, anche se è molto simile a $F$, questi sono i due grafici:
non capisco cosa sia sbagliato, sembra tipo traslato il polinomio interpolante rispetto a $F$...