models/johnpye/testconopt.a4c

REQUIRE "atoms.a4l";

MODEL testconopt;
	Profit, Inp, Out IS_A solver_var;
	P, L, K IS_A solver_var;

	Al, Ak, Ainp, Rho IS_A real_constant;
	Pinp, P0, D IS_A real_constant;
	W, L0, K0 IS_A real_constant;

	Out = ( Al*L^(-Rho) + Ak*K^(-Rho) + Ainp*Inp^(-Rho) )^(-1/Rho);

	Out = ( P0 - P ) / D;

	L0 < L;
	K = K0;
	0 <= P;
	0 < Out;
	0 < Inp;

	W :== 1.0;
	L0 :== 0.1;
	Pinp :== 1.0;
	Al :== 0.16;
	Ak :== 2.00;
	Ainp :== 0.16;
	Rho :== 1.0;
	K0 :== 4.0;
	P0 :== 2.0;
	D :== 0.5;

	Profit = P*Out - (L-L0)*W - Inp*Pinp;

METHODS

METHOD default_self;
END default_self;

METHOD specify;
END specify;

METHOD values;
END values;

METHOD on_load;
	RUN default_self;
	RUN specify;
	RUN values;
END on_load;

METHOD self_test;
	ASSERT L = 0.295686;
	ASSERT Inp = 0.295686;
	ASSERT Out = 0.632019;
	ASSERT P = 1.683990;	
END self_test;

END testconopt;
1	johnpye	786	REQUIRE "atoms.a4l";
2
3			MODEL testconopt;
4			Profit, Inp, Out IS_A solver_var;
5			P, L, K IS_A solver_var;
6
7			Al, Ak, Ainp, Rho IS_A real_constant;
8			Pinp, P0, D IS_A real_constant;
9			W, L0, K0 IS_A real_constant;
10
11			Out = ( AlL^(-Rho) + AkK^(-Rho) + Ainp*Inp^(-Rho) )^(-1/Rho);
12
13			Out = ( P0 - P ) / D;
14
15			L0 < L;
16			K = K0;
17			0 <= P;
18			0 < Out;
19			0 < Inp;
20
21			W :== 1.0;
22			L0 :== 0.1;
23			Pinp :== 1.0;
24			Al :== 0.16;
25			Ak :== 2.00;
26			Ainp :== 0.16;
27			Rho :== 1.0;
28			K0 :== 4.0;
29			P0 :== 2.0;
30			D :== 0.5;
31
32			Profit = POut - (L-L0)W - Inp*Pinp;
33
34			METHODS
35
36			METHOD default_self;
37			END default_self;
38
39			METHOD specify;
40			END specify;
41
42			METHOD values;
43			END values;
44
45			METHOD on_load;
46			RUN default_self;
47			RUN specify;
48			RUN values;
49			END on_load;
50
51			METHOD self_test;
52			ASSERT L = 0.295686;
53			ASSERT Inp = 0.295686;
54			ASSERT Out = 0.632019;
55			ASSERT P = 1.683990;
56			END self_test;
57
58			END testconopt;