ROSMIR

//Solucion tarea 2

exec('load_Elev.sci', -1) // data de las funciones [x, y, z]
exec('interpolate_grid.sci',-1); //Interpolacion de los datos x,y,z

//Coordenadas
[x,y,z]=load_Elev('surf-excacion.csv');

//Matriz con la data [x, y, z]
M= [x; y; z];

//Realizacion de la grilla
 grilla=interpolate_grid(M,1);



// para el calculo del volumen total se calcula de forma individual cada volumen y luego se suman.
// Volumen removido para la Excavacion

     V1 =   (x(1,1)* y(1,1) * z(1,1));
     V2 =   (x(1,2)* y(1,2) * z(1,2));
     V3 =   (x(1,3)* y(1,3) * z(1,3));
     V4 =   (x(1,4)* y(1,4) * z(1,4));
     V5 =   (x(1,5)* y(1,5) * z(1,5));
     V6 =   (x(1,6)* y(1,6) * z(1,6));
     V7 =   (x(1,7)* y(1,7) * z(1,7));
     V8 =   (x(1,8)* y(1,8) * z(1,8));
     V9 =   (x(1,9)* y(1,9) * z(1,9));
     V10 =  (x(1,10)* y(1,10) * z(1,10));
     V11 =  (x(1,11)* y(1,11) * z(1,11));
     V12 =  (x(1,12)* y(1,12) * z(1,12));
     V13 =  (x(1,13)* y(1,13) * z(1,13));
     V14 =  (x(1,14)* y(1,14) * z(1,14));
     V15 =  (x(1,15)* y(1,15) * z(1,15));
     V16 =  (x(1,16)* y(1,16) * z(1,16));
     V17 =  (x(1,17)* y(1,17) * z(1,17));
     V18 =  (x(1,18)* y(1,18) * z(1,18));
     V19 =  (x(1,19)* y(1,19) * z(1,19));
     V20 =  (x(1,20)* y(1,20) * z(1,20));
     V21 =  (x(1,21)* y(1,21) * z(1,21));
     V22 =  (x(1,22)* y(1,22) * z(1,22));
     V23 =  (x(1,23)* y(1,23) * z(1,23));
     V24 =  (x(1,24)* y(1,24) * z(1,24));
     V25 =  (x(1,25)* y(1,25) * z(1,25));
     V26 =  (x(1,26)* y(1,26) * z(1,26));
     V27 =  (x(1,27)* y(1,27) * z(1,27));
     V28 =  (x(1,28)* y(1,28) * z(1,28));
     V29 =  (x(1,29)* y(1,29) * z(1,29));
     V30 =  (x(1,30)* y(1,30) * z(1,30));
     V31 =  (x(1,31)* y(1,31) * z(1,31));
     V32 =  (x(1,32)* y(1,32) * z(1,32));
     V33 =  (x(1,33)* y(1,33) * z(1,33));
     V34 =  (x(1,34)* y(1,34) * z(1,34));
     V35 =  (x(1,35)* y(1,35) * z(1,35));
     V36 =  (x(1,36)* y(1,36) * z(1,36));
     V37 =  (x(1,37)* y(1,37) * z(1,37));
     V38 =  (x(1,38)* y(1,38) * z(1,38));
     V39 =  (x(1,39)* y(1,39) * z(1,39));
     V40 =  (x(1,40)* y(1,40) * z(1,40));
     V41 =  (x(1,41)* y(1,41) * z(1,41));
     V42 =  (x(1,42)* y(1,42) * z(1,42));
     V43 =  (x(1,43)* y(1,43) * z(1,43));
     V44 =  (x(1,44)* y(1,44) * z(1,44));
     V45 =  (x(1,45)* y(1,45) * z(1,45));
     V46 =  (x(1,46)* y(1,46) * z(1,46));
     V47 =  (x(1,47)* y(1,47) * z(1,47));
     V48 =  (x(1,48)* y(1,48) * z(1,48));
     V49 =  (x(1,49)* y(1,49) * z(1,49));
     V50 =  (x(1,50)* y(1,50) * z(1,50));
     V51 =  (x(1,51)* y(1,51) * z(1,51));
     V52 =  (x(1,52)* y(1,52) * z(1,52));
     V53 =  (x(1,53)* y(1,53) * z(1,53));
     V54 =  (x(1,54)* y(1,54) * z(1,54));
     V55 =  (x(1,55)* y(1,55) * z(1,55));
     V56 =  (x(1,56)* y(1,56) * z(1,56));
     V57 =  (x(1,57)* y(1,57) * z(1,57));
     V58 =  (x(1,58)* y(1,58) * z(1,68));
     V59 =  (x(1,59)* y(1,59) * z(1,59));
     V60 =  (x(1,60)* y(1,60) * z(1,60));
     V61 =  (x(1,61)* y(1,61) * z(1,61));
     V62 =  (x(1,62)* y(1,62) * z(1,62));
     V63 =  (x(1,63)* y(1,63) * z(1,63));
     V64 =  (x(1,64)* y(1,64) * z(1,64));
     V65 =  (x(1,65)* y(1,65) * z(1,65));
     V66 =  (x(1,66)* y(1,66) * z(1,66));
     V67 =  (x(1,67)* y(1,67) * z(1,67));
     V68 =  (x(1,68)* y(1,68) * z(1,68));
     V69 =  (x(1,69)* y(1,69) * z(1,69));
     V70 =  (x(1,70)* y(1,70) * z(1,70));
     V71 =  (x(1,71)* y(1,71) * z(1,71));
     V72 =  (x(1,72)* y(1,72) * z(1,72));
     V73 =  (x(1,73)* y(1,73) * z(1,73));
     V74 =  (x(1,74)* y(1,74) * z(1,74));
     V75 =  (x(1,75)* y(1,75) * z(1,75));
     V76 =  (x(1,76)* y(1,76) * z(1,76));
     V77 =  (x(1,77)* y(1,77) * z(1,77));
     V78 =  (x(1,78)* y(1,78) * z(1,78));

 VT=(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+V34+V35+V36+V37+V38+V39+V40+V41+V4+V43+V44+V45+V46+V47+V48+V49+V50+V51+V52+V53+V54+V55+V56+V57+V58+V59+V60+V61+V62+V63+V64+V65+V66+V67+V68+V69+V70+V71+V72+V73+V74+V75+V76+V77+V78);

// Comentarios para guiarnos a realizar la grilla: head_comments('interpolate_grid');

// Graficar los datos de x, y, z.
plot3d(grilla(1), grilla(2), grilla(3));