#using "types";

FlowSheet Alumina_Flow
	VARIABLES
	y as Real;
	x1 as Real;
	x2 as Real;

	PARAMETERS
	b0 as Real(Default=1);
	b1 as Real(Default=1);
	b2 as Real(Default=1);

	SPECIFY
	x1 = 0.647;
	x2 = 77.1;
	
	EQUATIONS
	y = b0 + b1*x1 + b2*x2;
	
	OPTIONS
	Dynamic = false;
	#NLASolver(AbsoluteAccuracy=1e-13);
end

Estimation Alumina_Est as Alumina_Flow

	ESTIMATE
	# PAR   START  LOWER  UPPER  UNIT
	b0     -1.00    -5  5;
	b1      2.00    -5  5;
	b2     -0.00    -5  5;

	EXPERIMENTS
	# FILE                WEIGTH	TYPE
# This sequence works correctly, but the simulation results go all to the first level of variables (fix me)
	"alumina_est.dat"		1		"fit";
	"alumina2_est.dat"		1		"test";

# This sequence does not show all data (fix me)
#	"alumina2_est.dat"		1		"fit";
#	"alumina_est.dat"      1		"test";

# This sequence does not compute statistics (fix me)
#	"alumina2_est.dat"		1		"test";
#	"alumina_est.dat"		1		"fit";

# This sequence has both problems above (fix me)
#	"alumina_est.dat"		1		"test";
#	"alumina2_est.dat"		1		"fit";

	OPTIONS
	Statistics(
		Fit=true,
		Parameter=true,
		Prediction=false
	);
	Significance = 0.95;
	BiLateral = true;
	NumJac = true;
	
	NLPSolver(
		AbsoluteAccuracy=1e-9,
		#File = "complex"
		File = "ipopt_emso"
	);
	Dynamic = false;
end