Sets C culturas /pastnat30, cfeno30, pastsem30, pastnat60, cfeno60, pastsem60, pastnat150, pastsem150, pastnat180, pastnat410/ A actividades de produção animal PRODI /BovM, BovML/ U unidades de utilização /U1/ F alimentos /pastn, feno, pasts/ FP(F) alimentos produzidos na empresa /pastn, feno, pasts/ Tr transferência de alimentos /T1-2, T2-3, T3-4, T4-1/ ZF Períodos de alimentação /R1*R4/ ZT trabalho alocado a actividade pecuária /vaqueiroM, vaqueiroML/ UC(U,C) crossing possibilities UUtilização com culturas /U1.pastnat30, U1.cfeno30, U1.pastsem30, U1.pastnat60, U1.cfeno60, U1.pastsem60, U1.pastnat150, U1.pastsem150, U1.pastnat180, U1.pastnat410/ CF(C,F) crossing possibilities culturas alimentos da empresa /pastnat30.pastn, cfeno30.feno, pastsem30.pasts, pastnat60.pastn, cfeno60.feno, pastsem60.pasts, pastnat150.pastn, pastsem150.pasts, pastnat180.pastn, pastnat410.pastn/ Parameter R(A) Rendimento das actividades pecuárias /BovM 284.6, BovML 310/; Parameter K(U,C) Ajuda por hectare em PRODI /U1.pastnat30 97.5, U1.cfeno30 36.8, U1.pastsem30 119.6, U1.pastnat60 78, U1.cfeno60 29.44, U1.pastsem60 95.7, U1.pastnat150 48.8, U1.pastsem150 59.8, U1.pastnat180 19.5, U1.pastnat410 0/ ; Scalar CV custo hora contrat vaqueiro /10.6/; Parameter ChaC(U,C) custo Total por ha cultura /U1.pastnat30 40, U1.cfeno30 439.8, U1.pastsem30 122, U1.pastnat60 40, U1.cfeno60 439.8, U1.pastsem60 122, U1.pastnat150 40, U1.pastsem150 122, U1.pastnat180 40, U1.pastnat410 40/; Table AF(U,C,F,ZF) produção alimento por cultura e por hectare e período R1 R2 R3 R4 U1.Pastnat30.pastn 250 180 700 U1.cfeno30.feno 4400 U1.Pastsem30.pasts 300 250 1700 U1.Pastnat60.pastn 250 180 700 U1.cfeno60.feno 4400 U1.Pastsem60.pasts 300 250 1700 U1.Pastnat150.pastn 250 180 700 U1.Pastsem150.pasts 300 250 1700 U1.Pastnat180.pastn 250 180 700 U1.Pastnat410.pastn 250 180 700 ; Table COF(U,A,F,ZF) Alim disponível por classe anim e período de alimentação R1 R2 R3 R4 U1.BovM.pastn -1 -1 -1 -1 U1.BovM.feno -1 -1 -1 -1 U1.BovM.pasts -1 -1 -1 -1 U1.BovML.pastn -1 -1 -1 -1 U1.BovML.feno -1 -1 -1 -1 U1.BovML.pasts -1 -1 -1 -1 ; Table TRF(F,TR,ZF) transferências de alimentos entre períodos R1 R2 R3 R4 pastn.T1-2 -1 1 pastn.T2-3 -1 1 pastn.T3-4 -1 0.85 feno.T3-4 -1 1 feno.T4-1 1 -1 feno.T1-2 -1 1 pasts.T1-2 -1 1 pasts.T2-3 -1 1 pasts.T3-4 -1 0.85 ; Table UFL(A,F,ZF) UFL por kgMS (past e feno) e kg alimento (ACC e TLuz) R1 R2 R3 R4 BovM.pastn 0.8 0.75 0.5 0.35 BovM.pasts 0.85 0.85 0.6 0.40 BovM.feno 0.6 0.6 0.6 0.6 BovML.pastn 0.8 0.75 0.5 0.35 BovML.pasts 0.85 0.85 0.6 0.40 BovML.feno 0.6 0.6 0.6 0.6 ; Table DR(A,F,ZF) Coeficiente matéria seca R1 R2 R3 R4 BovM.pastn -1 -1 -1 -1 BovM.feno -1 -1 -1 -1 BovM.pasts -1 -1 -1 -1 BovML.pastn -1 -1 -1 -1 BovML.feno -1 -1 -1 -1 BovML.pasts -1 -1 -1 -1 ; Table NUFL(A,ZF) necessidades UFL por unidade pecuária R1 R2 R3 R4 bovM -356 -550 -500 -600 bovML -373 -580 -526 -632 ; Table CMI(A,ZF) capacidade máx ing por unid pecuária R1 R2 R3 R4 bovM 695 1040 1260 1388 bovML 704 1053 1275 1418 ; Parameter HT(A,ZT) horas trabalho por unid pecuária /bovM.vaqueiroM 13, bovML.vaqueiroML 13/ ; Variables Z resultado de exploração X(U,C) ha por cultura Q(A) nº unidades pecuárias Transf(F,TR) transferência de alimentos Consm(A,F,ZF) consumo de alimentos Contvaq(ZT) contrat vaqueiro ; Positive variables X,Q,Transf,Consm,Contvaq ; Equations income função objectivo land1 restrição área feno land2 restrição área pastsem land3 restrição área total ajPRODIn30 ajuda haPRODI pastnat 30 ajPRODIf30 ajuda haPRODI cfeno 30 ajPRODIs30 ajuda haPRODI pastsem 30 ajPRODIn60 ajuda haPRODI pastnat 60 ajPRODIf60 ajuda haPRODI cfeno 60 ajPRODIs60 ajuda haPRODI pastsem 60 ajPRODIn150 ajuda haPRODI pastnat 150 ajPRODIs150 ajuda haPRODI pastnat 150 ajPRODIn180 ajuda haPRODI pastnat 180 ajPRODIn410 ajuda haPRODI pastnat 410 CNland Blabor eq trab para bovinos Feedbal(F,ZF) restrição de alimentos UFLBal(A,ZF) UFL balanço de consumo alimentos Maxing(A,ZF) cap máx ing ; INCOME .. Sum(A,Q(A)*R(A))+ sum((U,C),K(U,C)*X(U,C))-(Contvaq('vaqueiroM')*CV)- (Contvaq('vaqueiroML')*CV)- Sum((U,C),ChaC(U,C)*X(U,C)) =E= Z ; Land1 .. Sum((U)$UC(U,'cfeno30'),X(U,'cfeno30'))+ Sum((U)$UC(U,'cfeno60'),X(U,'cfeno60')) =L= 50 ; Land2 .. Sum((U)$UC(U,'pastsem60'),X(U,'pastsem60')) +Sum((U)$UC(U,'pastsem60'),X(U,'pastsem60')) +Sum((U)$UC(U,'pastsem150'),X(U,'pastsem150')) =L= 72 ; Land3 .. Sum(U,X(U,'pastnat30'))+Sum(U,X(U,'cfeno30'))+Sum(U,X(U,'pastsem30'))+ Sum(U,X(U,'pastnat60'))+Sum(U,X(U,'cfeno60'))+Sum(U,X(U,'pastsem60')) +Sum(U,X(U,'pastnat150'))+Sum(U,X(U,'pastsem150'))+Sum(U,X(U,'pastnat180')) +Sum(U,X(U,'pastnat410')) =L= 410 ; ajPRODIn30 .. sum(U,X(U,'pastnat30')) =L= 30 ; ajPRODIf30 .. sum(U,X(U,'cfeno30')) =L= 30 ; ajPRODIs30 .. sum(U,X(U,'pastsem30')) =L= 30 ; ajPRODIn60 .. sum(U,X(U,'pastnat60')) =L= 30 ; ajPRODIf60 .. sum(U,X(U,'cfeno60')) =L= 20 ; ajPRODIs60 .. sum(U,X(U,'pastsem60')) =L= 30 ; ajPRODIn150 .. sum(U,X(U,'pastnat150')) =L= 90 ; ajPRODIs150 .. sum(U,X(U,'pastsem150')) =L= 12 ; ajPRODIn180 .. sum(U,X(U,'pastnat180')) =L= 30 ; ajPRODIn410 .. sum(U,X(U,'pastnat410')) =L= 230 ; CNland .. Sum(A,Q(A)*1.2) - Sum(U,X(U,'pastnat30'))-Sum(U,X(U,'cfeno30'))-Sum(U,X(U,'pastsem30'))- Sum(U,X(U,'pastnat60'))-Sum(U,X(U,'cfeno60'))-Sum(U,X(U,'pastsem60')) -Sum(U,X(U,'pastnat150'))-Sum(U,X(U,'pastsem150'))-Sum(U,X(U,'pastnat180')) =E= 0 ; BLabor .. Sum((A),HT(A,'vaqueiroM')*Q('bovM'))+ Sum((A),HT(A,'vaqueiroML')*Q('bovML')) =L= (contvaq('vaqueiroM')+ contvaq('vaqueiroML')) ; Feedbal(F,ZF) .. Sum((U,C),AF(U,C,F,ZF)*X(U,C)$UC(U,C))+ Sum(TR,TRF(F,TR,ZF)*Transf(F,TR))+ Sum((U,A),COF(U,A,F,ZF)*Consm(A,F,ZF)$A(A)) =G= 0 ; UFLBAL(A,ZF) .. NUFL(A,ZF)*Q(A)+Sum(FP,UFL(A,FP,ZF)*consm(A,FP,ZF)$A(A)) =G= 0 ; Maxing(A,ZF) .. (CMI(A,ZF)*Q(A))+Sum(FP,DR(A,FP,ZF)*Consm(A,FP,ZF)$A(A)) =G= 0 ; model PRODI/all/ ; Solve PRODI using LP maximizing Z;