#Métodos Numéricos - 2023. IMERL, Fing
#Código de ejemplo para el ejercicio 7 del práctico 2
clear all;

xi=0;
xf=5;
alfa=0;
beta=sin(5);

N=50;
h=(xf-xi)/N;
i=1:N-1;
H=(xi+i*h).';

A=zeros(N-1,N-1);
b=zeros(N-1,1);

%Cargar A y b
%Cargo los bordes
b(1,1)=h^2*sin(H(1,1))*(e^(H(1,1))-1)-alfa;
b(N-1,1)=h^2*sin(H(N-1,1))*(e^(H(N-1,1))-1)-beta;

x=H(1,1);
A(1,1)=(e^x)*h^2-2;
A(1,2)=1;
x=H(N-1,1);
A(N-1,N-1)=(e^x)*h^2-2;
A(N-1,N-1-1)=1;


for j=2:N-2
  x=H(j,1);
  A(j,j-1)=1;
  A(j,j)=(e^x)*h^2-2;
  A(j,j+1)=2;
  b(j,1)=h^2*sin(x)*(e^x-1);
endfor

%Resolver con bslash
y1=A\b;
%Agregarle bordes
y1=[alfa,y1.',beta];

%Resolver con tridisolve
a3s=[1,diag(A,-1).',1];
b3s=[1,diag(A).',1];
c3s=a3s;
d3s=[alfa,b.',beta];
y2=tridisolve(a3s,b3s,c3s,d3s);

%Graficar
clear figure;
H=[xi,H.',xf];
hold on;
plot(H,sin(H),'LineWidth',3);
plot(H,y2,'LineWidth',3);
#plot(H,y1,'LineWidth',3);
legend('Solución exacta', 'Tridisolve','FontSize',28)