%Este programa pretende mostrar cómo quedaría la comparación entre
%felicidad y bienestar, aplicando la expresión propuesta.

close all;clear all

%Definimos el primer caso, con una ganacia de bienestar.
xb=[-10:.1:10];
bienestar=[ones(1,numel(xb))];
bienestar(round(numel(xb)/2+1):end)=2;
xf=linspace(min(xb),max(xb),numel(xb)-1);

felicidad=ones(size(xf))*NaN;
for ii=1:numel(xb)-1
    if xb(ii)<=0
        felicidad(ii)=diff(bienestar(ii:ii+1));%Antes de producirse cualquier cambio de bienestar
    else
        felicidad(ii)=diff(bienestar(ii:ii+1))+exp(-xf(ii));%Una vez se produce el cambio de bienestar
    end
end

figure
pos=get(gcf,'position');
pos(4)=pos(4)*2;
pos(3)=pos(3)*2;
pos(1)=pos(1)-100;
pos(2)=pos(2)-400;
set(gcf,'position',pos,'paperpositionmode','auto')
h1=plot(xb,bienestar,'linewidth',2);
hold on
h2=plot(xf,felicidad,'r','linewidth',2);
axis([-10.5 10.5 -2.5 2.5])
xlabel('Tiempo','fontsize',18)
set(gca,'ygrid','on','fontsize',16)
legend([h1 h2],'Bienestar (\beta)','Felicidad (\phi)',4)
set(gca,'ygrid','on')
fich_sal='C:\MisDocumentos\AdHoc\Docs\InvestigaBlog\bienestar_felicidad_positivo';
eval(['export_fig ',fich_sal,' -png -transparent -a4 -opengl'])

%Definimos el segundo caso, con una pérdida de bienestar.
xb=[-10:.1:10];
bienestar=[ones(1,numel(xb))];
bienestar(1:round(numel(xb)/2))=2;
xf=linspace(min(xb),max(xb),numel(xb)-1);

felicidad=ones(size(xf))*NaN;
for ii=1:numel(xb)-1
    if xb(ii)<=0
        felicidad(ii)=diff(bienestar(ii:ii+1));%Antes de producirse cualquier cambio de bienestar
    else
        felicidad(ii)=diff(bienestar(ii:ii+1))-exp(-xf(ii));%Una vez se produce el cambio de bienestar
    end
end

figure
pos=get(gcf,'position');
pos(4)=pos(4)*2;
pos(3)=pos(3)*2;
pos(1)=pos(1)-100;
pos(2)=pos(2)-400;
set(gcf,'position',pos,'paperpositionmode','auto')
h1=plot(xb,bienestar,'linewidth',2);
hold on
h2=plot(xf,felicidad,'r','linewidth',2);
axis([-10.5 10.5 -2.5 2.5])
xlabel('Tiempo','fontsize',18)
set(gca,'ygrid','on','fontsize',16)
legend([h1 h2],'Bienestar (\beta)','Felicidad (\phi)',4)

fich_sal='C:\MisDocumentos\AdHoc\Docs\InvestigaBlog\bienestar_felicidad_negativo';
eval(['export_fig ',fich_sal,' -png -transparent -a4 -opengl'])