Cilj ovog rada bio je izraditi MATLAB program za procjenu parametara (omjera kopolimerizacijskih rektivnosti) za slučajeve terpolimerizacija i izvršiti procjenu parametara za dva terpolimerizacijska sustava. Za izradu programa kao model procesa korištene su modificirane Alfrey-Goldfingerove jednadžbe dok je optimizacijski algoritam napravljen korištenjem modela grešaka u varijablama i genetičkog algoritma. Na oba sustava primijenjen je dizajn eksperimenata pri procjeni parametara a model je provjeren na podacima koji nisu korišteni prilikom procjene parametara. Prvi terpolimerizacijski sustav sastoji se od dimetilaminoetil-metakrilata (DMAEM), metil-metakrilata (MMA) i dodecil-metakrilata (DDMA). Eksperimentalni podaci za prvi sustav preuzeti su iz literature a procjena parametara za taj sustav izvršena je pomoću vlastito izrađenog programa. Rezultati procjene parametara za prvi sustav su r12= 0,73, r21= 0,82, r13= 0,91, r31= 0,89, r23= 1,32, r32= 1,34 što je blisko parametrima koji su navedeni u literaturi. Usprkos dobroj procjeni parametara, literaturni podaci o parametrima određenim u zasebnim binarnim eksperimentima daju bolje rezultate. Za prvi sustav simulacijom su pronađena šira pseudoazeotropna područja a postojanje jednog područja je i eksperimentalno potvrđeno. Drugi sustav sastoji se od metil-metakrilata (MMA), 2-(tert-butilamino)etil-metakrilata (TBAEMA) i oktadecil-metakrilata (ODMA). Za drugi sustav eksperimentalni podaci su prikupljeni provođenjem vlastitih polimerizacija u toluenu i analizom sastava kopolimera 1H NMR spektroskopijom. Rezultati procjene parametara za drugi sustav su r12= 0,84, r21= 0,90, r13= 1,75, r31= 0,62, r23= 1,10, r32= 0,77. Vlastito procijenjeni parametri bolje opisuju eksperimentalne podatke nego literaturno dostupni parametri iz zasebnih binarnih eksperimenata. Za drugi sustav simulacijom nisu pronađena pseudoazeotropna područja od većeg značaja. The aim of this work was to develop a MATLAB program for parameter estimation (ratios of copolymerization reactivities) for cases of terpolymerizations and to perform parameter estimation for two terpolymerization systems. Modified Alfrey-Goldfinger equations were used as process model, while the optimization algorithm was created using the error in variables model and the genetic algorithm. On both systems, the design of the experiments was applied when estimating the parameters and the model was tested on data that were not used when estimating the parameters. The first terpolymerization system consists of dimethylaminoethyl- methacrylate (DMAEM), methyl- methacrylate (MMA) and dodecyl- methacrylate (DDMA). Experimental data for the first system were taken from the literature and the estimation of parameters for this system was performed using a self-developed program. The results of parameter estimation for the first system are r12= 0.73, r21= 0.82, r13= 0.91, r31= 0.89, r23= 1.32, r32= 1.34 which is close to the parameters listed in the literature. Despite a good estimate of the parameters, the literature data on the parameters determined in separate binary experiments give better results. For the first system, wider pseudoazeotropic areas were found by simulation, and the existence of one area was experimentally confirmed. The second system consists of methylmethacrylate (MMA), 2-(tert-butylamino)ethyl- methacrylate (TBAEMA) and octadecyl- methacrylate (ODMA). For the second system, experimental data were collected by conducting polymerizations in toluene and analyzing the composition of the copolymer by 1H NMR spectroscopy. The results of parameter estimation for the second system are r12 = 0.84, r21 = 0.90, r13 = 1.75, r31 = 0.62, r23 = 1.10, r32 = 0.77. Self-estimated parameters better describe experimental data than parameters available in literature from separate binary experiments. For the second system, no significant pseudoazeotropic areas were found by simulation.