%   Copyright (c) 2001-2009 by Ales Cerny

%************************************************************************%
% chapter3sect5b.m - supplementary program to                            %
% Ales Cerny (2009) Mathematical Techniques in Finance (2nd ed.)         % 
% Princeton University Press http://press.princeton.edu/titles/9079.html %
%************************************************************************%
    
% This code is provided 'as-is', without any express or implied warranty.  
% 
% Permission is granted to anyone to use this code for any purpose, 
% subject to the following restrictions:
% 
% 1. The origin of this code must not be misrepresented; you must not
%    claim that you wrote the original code.
% 2. Modified code versions must be plainly marked as such, and must not 
%    be misrepresented as being the original code.
% 3. This notice may not be removed from any source distribution.

% NOTICE TO STUDENTS: To avoid accusations of plagiarism, if you use this 
% code or its modifications in assessed work you should prepend it with a 
% note stating: 
%   "This is the original/modified version of the code chapter2sect4a.m by 
%    Ales Cerny (2009), Mathematical Techniques in Finance (2nd ed.),  
%    Princeton University Press. The original version is available from 
%    http://www.martingales.info/mtfweb2".
% A similar acknowledgement should appear prominently inside your written 
% report.

clear;
clc

gamavect=[0.5 1.00001 2 5 15 10^8];
X=[0.18 -0.12];
XDistr=[0.5 0.5];
IPTolerance = 10^(-7);

%****************************%
% main body of the programme %
%****************************%

IP=zeros(size(gamavect));
beta=zeros(size(gamavect));

for ii = 1:size(gamavect,2);
    [IP(ii),beta(ii)]=HARAmax(X,XDistr,gamavect(ii),IPTolerance);
end;
gamavect(6)=inf;

disp('______________________________________________________________________________________________________________________________')                                          
disp(' ')  
disp(sprintf('%s',['gamma                            ' num2str(gamavect,'%12.1f')]))
disp(' ')  
disp(sprintf('%s',['normalized optimal investment   ' num2str(beta,'%12.3f')]))
disp(' ')
disp(' ') 
disp(sprintf('%s',['investment potential            ' num2str(IP,'%12.4f')]))
disp('______________________________________________________________________________________________________________________________')  
disp(' ') 
disp(' ') 
disp(' ')