#************************************************************* # Test Problem 1 from Chapter 11, section 3.3 # Simple First Order Irreversible Reaction #************************************************************* OPTION {{ FTOL = 1e-10; OTOL = 1e-10; }} DECLARATION {{ INDEX {m,i,p}; SET M = |1:10|; # Number of Data sets SET I = |1:2|; # Number of Variables SET P = |1:2|; # Number of Parameters XVAR {x(M,I), # Fitted data variables a(P) # Parameters }; ZVAR {z(I)}; ISPE {z(I)}; ZINC {1,0}; TIME {0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0}; # Observed data PARA xe(M,I) = {0.606,0.373,0.368,0.564,0.223,0.647,0.135,0.669,0.082,0.656, 0.050,0.624,0.030,0.583,0.018,0.539,0.011,0.494,0.007,0.451}; # Variable Bounds LBDS x(M,I) = < m E M | < i E I | 0 > >; UBDS x(M,I) = < m E M | < i E I | 1 > >; LBDS a(P) = < p E P | 0 >; UBDS a(P) = < p E P | 10 >; }} MODEL {{ MIN: <> >>; dae1: z'(1) =e= -a(1)*z(1); dae2: z'(2) =e= a(1)*z(1) - a(2)*z(2); err(m E M,i E I)[m]: z(i) - x(m,i) =e= 0; }}