PESOPROPIO.m

function[R,D]=PESOPROPIO(Lt,Lc,Dnt,Dnc);

phor=2.5/(100*100*100); %[ton/m^3]
pac=7.8/(100*100*100);%[ton/m^3]
l1=3232.08;
l2=2740;
l3=2262.74;
l4=1812.07;
l5=1341.64;
d=2.5;
a=2*pi*(d/2)^2;

acs1=18150;
acs2=19800;
acs3=690*110;
acs4=110^2*2;


Pcb1=a*pac*l1;
Pcb2=a*pac*l2;
Pcb3=a*pac*l3;
Pcb4=a*pac*l4;
Pcb5=a*pac*l5;


At=60*45*2+25*(450-90); 
Vtm=At*100;
Vte=At*50;
Ptm=Vtm*phor;
Pte=Vte*phor;


c=zeros(1,4);
cc=zeros(1,4);
c2=zeros(1,11); 
D=[Pcb1/2,c,Pcb2/2,c,Pcb3/2,c,Pcb4/2,c,Pcb5/2,c2,Pcb5/2,c,Pcb4/2,c,Pcb3/2,c,Pcb2/2,c,Pcb1/2,cc, Pcb1/2,c,Pcb2/2,c,Pcb3/2,c,Pcb4/2,c,Pcb5/2,c2,Pcb5/2,c,Pcb4/2,c,Pcb3/2,c,Pcb2/2,c,Pcb1/2];
%Columns

Pcs1=acs1*80*phor;
Pcs2=acs2*80*phor;
Pcs3=acs3*80*phor;
Pcs4=acs4*80*phor;
R(1)=-Pte-D(1);
R(110)=-Pte-D(110);
for i=2:109 
    if (D(i)~=0)
        R(i)=-(D(i)+Ptm);
    else
        R(i)=-Ptm;
    end
end
R(111:242)=0;
R=R';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
d2=2.5/10^6;
f=82.571;
h=81.7*110*110*2;
V34=0;
V33=165*110*80+Pcb1;
V32=2*V33;
V31=3*V33+Pcb2;
V30=4*V33;
V29=5*V33+Pcb3;
V28=6*V33;
V27=7*V33+Pcb4;
V26=8*V33;
V25=9*V33;
V24=((V33/2)+(2*f/2)*90*110+V25)+Pcb5;
V23=(V24+2*f*90*110);
V22=(V23+2*f*90*110);
V21=(V22+2*f*90*110);
V20=(V21+2*f*90*110);
V19=(V20+2*f*90*110);
V18=(V19+2*f*90*110);
V17=(V18+2*f*90*110);
V16=(V17+2*f*90*110);
V15=(V16+2*f*90*110);
V14=(V15+2*f*90*110);
V13=(V14+2*f*90*110);
V12=(V13+2*f*90*110);
V11=(V12+2*f*90*110);
V10=(V11+2*f*90*110);
V9=(V10+2*f*90*110);
V8=(V9+f*90*110+40*(690-17.91)*110);
V7=(V8+30*110*(690-8.96)+50*(690-10.38)*110);
V6=(V7+20*110*(669.24-4.15)+3/4*81.7*110*110*2);
V5=(V6+h);
V4=(V5+h);
V3=(V4+h);
V2=(V3+h);
V1=(V2+h)
V=-d2*[V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 V20 V21 V22 V23 V24 V25 V26 V27 V28 V29 V30 V31 V32 V33];
R(210:242)=V;
R(177:209)=V;
return