Помощь в написании скрипта на shell

[mrceye]

В универе дали следующее задание:

Написать консольную программу для табулирования функции (например, sin(x)) на отрезке [a,b] с шагом h=(b-a)/1000000 с хранением значений таблицы в массиве. Проанализировать скорость работы программы при выводе на экран каждого элемента массива, каждого сотого, каждого десятитысячного.

Собственно у меня получился следующий код:

Код: [Выделить]

#!/bin/bash
a=$1
b=$2
h=$(echo "("$2"-"$1")/1000000" | bc -lq)

x=$a
i=0

while (($x >= $h)); do
x=$(echo "s("$x")" | bc -lq)
mas[$i]=$x
echo ${mas[$i]}
i=$((i+1))
done
При выполнении получаю следующую ошибку:

Код: [Выделить]

dim@home-pc:~$ bash prog.sh 1 2
prog.sh: строка 9: ((: 1 >= .00000100000000000000: ошибка синтаксиса: ожидается операнд (error token is ".00000100000000000000")
Как я понял проблема в программе bc, так как она возвращает значение вида .xxx, хотя может быть ошибаюсь. Какие могут быть способы для решения этой проблемы? Или может быть стоит стоит воспользоваться какой-нибудь другой программой для расчета?

Заранее спасибо за помощь! 

[wl]

У bash-а целочисленная арифметика, он не понимает floating point.
гугль по словам bash floating-point arithmetic выдает несколько вариантов решения проблемы.

Да, а bash - это Ваш выбор или условие задачи?
Просто я бы что-нибудь позаковыристей выбрал. R при возможности .

awk тоже неплохо справляется:

Код: [Выделить]

$ awk 'BEGIN {if(1>=0.1)print 10/4}'
2.5
Про perl, python, ruby вообще молчу (они в стандартной поставке)

[mrceye]

Да, bash — условие задачи.
Хм, пожалуй попробую awk. Спасибо за помощь

[alexander.pronin]

Просто я бы что-нибудь позаковыристей выбрал. R при возможности .

Какая жесть.

[wl]

Просто я бы что-нибудь позаковыристей выбрал. R при возможности .

Какая жесть.

Не надо грязи! R мне на прошлом месте работы 6 лет служил верой и правдой, будучи основной рабочей лошадкой.