#Métodos Numéricos - 2023. IMERL, Fing
#Código de ejemplo del ejercicio 9 del práctico 1
clear all;
imax=15;

#Parte b
p=ones(1,imax);
for i=1:imax-1
  p(1,i+1)=e-(i+1)*p(1,i);
endfor
figure(1)

#Parte c
p18=ones(1,imax)./8;
p20=20.*ones(1,imax);
p100=100.*ones(1,imax);
for i=imax:-1:2
  p18(1,i-1)=(e-p18(1,i))/i;
  p20(1,i-1)=(e-p20(1,i))/i;
  p100(1,i-1)=(e-p100(1,i))/i;
endfor

#Graficar sucesión calculada hacia adelante
figure(1);
plot(1:imax,p,'Color',[0,0,1],'LineStyle','-','marker','*','LineWidth', 1);

#Comparación cálculo hacia adelante y hacia atrás
#Sucesión hacia atrás con p20=1/8
figure(2);
plot(1:imax,p,'Color',[0.5,0,0.5],'LineStyle','-','marker','*','LineWidth', 1);
hold on;
plot(1:imax,p18,'Color',[0,0.5,0],'LineStyle','-','marker','*','LineWidth', 1);
legend({"Hacia adelante","Hacia atrás"},'FontSize',50,'location','southwest');
title("Comparación cálculo hacia adelante y hacia atrás",'FontSize',50);

figure(3);
plot(1:imax,p18,'Color',[1,0,0],'LineStyle','-','marker','*','LineWidth', 1);
hold on;
plot(1:imax,p20,'Color',[0,1,0],'LineStyle','-','marker','*','LineWidth', 1);
plot(1:imax,p100,'Color',[0,0,1],'LineStyle','-','marker','*','LineWidth', 1);
legend ({"0.125", "20","100"},'FontSize',50,'location','northwest');
title("Comparación diferentes p20",'FontSize',50)