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 371, BovML 395.5/; 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 VARIÁVEL por ha cultura /U1.pastnat30 0.0, U1.cfeno30 305.7, U1.pastsem30 39.5, U1.pastnat60 0.0, U1.cfeno60 305.7, U1.pastsem60 39.5, U1.pastnat150 0.0, U1.pastsem150 39.5, U1.pastnat180 0.0, U1.pastnat410 0.0/; 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 -344 -535 -486 -584 bovML -360 -562 -510 -612 ; Table CMI(A,ZF) capacidade máx ing por unid pecuária R1 R2 R3 R4 bovM 707 1061 1282 1370 bovML 718 1078 1300 1390 ; 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.18) - 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;