% Factor de descuento dada tasa de interes r
d = @(r) (1+r)^(-1) ;
% Rendimiento a partir de t = 2
w2 = @(r,x) (d(r)^2)*(1-d(r))^(-1)*x ;
% Rendimiento a partir de t = 1
w1 = @(r,x) d(r) * x + 0.5 * w2(r,x+10) + 0.5 * w2(r,x-10);
% Rendimiento a partir de t = 0
w0 = @(r) 30 + 0.5 * w1(r,40) + 0.5 * w1(r,20) ;
% Rendimiento neto si invertimos en t = 0
v0 = @(r) -300 + w0(r);
% Rendimiento neto si invertimos en t = 1 con X1 = x
v1 = @(r,x) -300*d(r) + w1(r,x);
% Rendimiento neto si invertimos en t = 2, con X2 = x
v2 = @(r,x) -300*(d(r)^2) + w2(r,x);
fprintf('Rendimiento neto inversion en t = 0: %f\n',v0(.1));
fprintf('Rendimiento neto inversion en t = 1, x1 = 20: %f\n',v1(.1,20));
fprintf('Rendimiento neto inversion en t = 1, x1 = 40: %f\n',v1(.1,40));
fprintf('Rendimiento neto inversion en t = 2, x2 = 10: %f\n',v2(.1,10));
fprintf('Rendimiento neto inversion en t = 2, x2 = 30: %f\n',v2(.1,30));
fprintf('Rendimiento neto inversion en t = 2, x2 = 50: %f\n',v2(.1,50));