calc.sh

   1 #!/bin/sh
   2 # vim: set tabstop=4 softtabstop=4 shiftwidth=4 expandtab:
   3 #
   4 # Trabajo Práctico I de Análisis Numérico I
   5 # Realiza los cálculos del punto c (y a y b).
   6 # Copyright (C) 2002 Leandro Lucarella <leandro@lucarella.com.ar>
   7 #
   8 # Este programa es Software Libre; usted puede redistribuirlo
   9 # y/o modificarlo bajo los términos de la "GNU General Public
  10 # License" como lo publica la "FSF Free Software Foundation",
  11 # o (a su elección) de cualquier versión posterior.
  12 #
  13 # Este programa es distribuido con la esperanza de que le será
  14 # útil, pero SIN NINGUNA GARANTIA; incluso sin la garantía
  15 # implícita por el MERCADEO o EJERCICIO DE ALGUN PROPOSITO en
  16 # particular. Vea la "GNU General Public License" para más
  17 # detalles.
  18 #
  19 # Usted debe haber recibido una copia de la "GNU General Public
  20 # License" junto con este programa, si no, escriba a la "FSF
  21 # Free Software Foundation, Inc.", 59 Temple Place - Suite 330,
  22 # Boston, MA  02111-1307, USA.
  23 #
  24 # $URL$
  25 # $Date$
  26 # $Rev$
  27 # $Author$
  28 #
  29
  30 # Calculo soluciones del punto a, b y c.
  31 for (( n=4; n <= 32; n*=2 )); do
  32     ./77891 $n | tail -n1 | ./resultado.awk > ca$n.txt
  33 done
  34
  35 # Calculo iteraciones y S en función de w.
  36 for (( n=4; n <= 32; n*=2 )); do
  37     echo "Procesando N = $n ..."
  38     for (( i=0; i < 100; i+=1 )); do
  39         w=`printf "1.%02d\n" $i`
  40         echo -n "$w "
  41         # Iteraciones en función de w.
  42         ./77891 $n $w | wc -l | awk '{printf "%-4d", $1-1}'
  43         echo -n " "
  44         # S en función de w.
  45         ./77891 $n $w | tail -n1 | awk '{print $3}'
  46     done > cb$n.txt
  47 done
  48
  49 # Calculo los puntos e-g
  50 ./77891 16 1.35 0.001 | awk '{print $1 " " log($2) " " $3 " " $116}' > eg.txt
  51 # Calculo los puntos f-g
  52 ./77891 16 1.85 0.001 | awk '{print $1 " " log($2) " " $3 " " $116}' > fg.txt
  53
  54 # Calculo la Parte II
  55 ej2="./77891 8 1.446467 0.001 71 79 78 77"
  56 $ej2 | tail -n1 | ./isolineas.awk > isolineas.txt
  57 $ej2 | tail -n1 | awk '{ printf "Radio espectral: %.2f\nX25 = %.3f ± 0.001%%\n", $3, $28 }' > parte2.txt