//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-nat.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 nivelacion.
 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));

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);


........................................................................
//Realizacion del grafico por grillas
plot3d(grilla(1),grilla(2),grilla(3));