-
Notifications
You must be signed in to change notification settings - Fork 1
/
MLIAM11.m
36 lines (25 loc) · 1013 Bytes
/
MLIAM11.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
%MLIAM08
clear all
ctax0 = 0;
citer0 = [0 0 0];
NPV0 = -fMLIAM11(0.0, citer0, 1);
options = optimset('MaxFunEvals',1000,'MaxIter',1000);
ctax1 = fminsearch(@(ctax) fMLIAM11(ctax,citer0,1), ctax0, options);
ctax2 = fminsearch(@(ctax) fMLIAM11(ctax,citer0,2), ctax0, options);
ctax3 = fminsearch(@(ctax) fMLIAM11(ctax,citer0,3), ctax0, options);
citer1(1) = ctax1;
citer1(2) = ctax2;
citer1(3) = ctax3;
ctax1 = fminsearch(@(ctax) fMLIAM11(ctax,citer1,1), ctax0, options);
ctax2 = fminsearch(@(ctax) fMLIAM11(ctax,citer1,2), ctax0, options);
ctax3 = fminsearch(@(ctax) fMLIAM11(ctax,citer1,3), ctax0, options);
citer2(1) = ctax1;
citer2(2) = ctax2;
citer2(3) = ctax3;
ctax1 = fminsearch(@(ctax) fMLIAM11(ctax,citer2,1), ctax0, options);
ctax2 = fminsearch(@(ctax) fMLIAM11(ctax,citer2,2), ctax0, options);
ctax3 = fminsearch(@(ctax) fMLIAM11(ctax,citer2,3), ctax0, options);
citer3(1) = ctax1;
citer3(2) = ctax2;
citer3(3) = ctax3;
%NPV1 = -fMLIAM08(ctax1,ctax0,1);