Bc programming language

bc is "an arbitrary precision calculator language" with syntax similar to the C programming language. It is generally used by typing the command  on a Unix command prompt and entering a mathematical expression, such as (1 + 3) * 2, whereupon 8 will be outputted.

There are currently two main dialects: the rigorously defined POSIX bc, and its direct descendant, the much expanded GNU bc (also, GNU bc is available for a wider range of platforms, such as Microsoft Windows).

Both forms of bc can be executed as either a mathematical scripting language or as an interactive mathematical shell.

POSIX bc
The POSIX standardised bc language is traditionally written as a program in the dc programming language to provide a higher level of access to the features of the dc language without the complexities of dc's terse syntax.

In this form, the bc language contains single letter variable, array and function names and most standard arithmetic operators as well as the familiar control flow constructs, (,  and  ) from C. Unlike C, an   clause may not be followed by an  .

Functions are defined using a   keyword and values are returned from them using a   followed by the return value in parentheses. The   keyword (optional in C) is used to declare a variable as local to a function.

All numbers and variable contents are fixed precision floating-point numbers whose precision (in decimal places) is determined by the global   variable.

The numeric base of input (in interactive mode), output and program constants may be specified by setting the reserved   (input base) and   (output base) variables.

Output is generated by deliberately not assigning the result of a calculation to a variable.

Comments may be added to bc code by use of the C   and   (start and end comment) symbols.

Exactly as C
The following POSIX bc operators behave exactly like their C counterparts: +    -     *     / +=    -=    *=    /= ++    --    <     >

==   !=    <=    >=

[ ]  { } 

Similar to C
The modulus operators:

%    %=

... behave exactly like their C counterparts only when the global   variable is set to 0, i.e. all calculations are integer-only. When   is greater than 0 the modulus is calculated relative to the smallest positive value greater than zero.

Only resembling C
The operators:

^    ^=

... resemble the C bitwise exclusive-or operators, but are in fact the bc integer exponentiation operators.

'Missing' operators relative to C
The bitwise, boolean and conditional operators:

&    |     ^     &&    ||    ^^

&=   |=    ^=    &&=   ||=   ^^=

<<   >>

<<=  >>=

?:

... are not available in POSIX bc.

Built-in functions
The   function for calculating square roots is POSIX bc's only built-in mathematical function. Other functions are available in an external standard library.

Standard library functions
bc's standard library contains functions for calculating sine, cosine, arctangent, natural logarithm, the exponential function and the two parameter Bessel function J. Most standard mathematical functions (including the other inverse trigonometric functions) can be constructed using these. See external links for implementations of many other functions.

GNU bc
GNU bc derives from the POSIX standard and includes many enhancements. It is entirely separate from dc-based implementations of the POSIX standard and is instead written in C. Nevertheless, it is fully backwards compatible as all POSIX bc programs will run unmodified as GNU bc programs.

GNU bc variables, arrays and function names may contain more than one character, some more operators have been included from C, and notably, an   clause may be followed by an  .

Output is achieved either by deliberately not assigning a result of a calculation to a variable (the POSIX way) or by using the added   statement.

Furthermore, a   statement allows the interactive input of a number into a running calculation.

In addition to C-style comments, a   character will cause everything after it until the next new-line to be ignored.

The value of the last calculation is always stored within the additional built-in   variable.

Extra operators
The following logical operators are additional to those in POSIX bc:

&&    ||      !

... and are available for use in conditional statements (such as within an   statement). Note, however, that there are still no equivalent bitwise or assignment operations.

Functions
All functions available in GNU bc are inherited from POSIX. No further functions are provided as standard with the GNU distribution.

Example code
Since the bc   operator only allows an integer power to its right, one of the first functions a bc user might write is a power function with a floating point exponent. Both of the below assume the standard library has been included:

A 'Power' function in POSIX bc
/* A function to return the integer part of x */ define i(x) { auto s   s = scale scale = 0 x /= 1  /* round x down */ scale = s   return (x) }

/* Use the fact that x^y == e^(y*log(x)) */ define p(x,y) { if (y == i(y)) { return (x ^ y)   } return ( e( y * l(x) ) ) }

An equivalent 'Power' function in GNU bc
define int(number) { auto oldscale oldscale = scale scale = 0 number /= 1 /* round number down */ scale = oldscale return number } define power(number,exponent) { if (exponent == int(exponent)) { return number ^ exponent } else { return e( exponent * l(number) ) } }
 * 1) A function to return the integer part of a number
 * 1) Use the fact that number^exponent == e^(exponent*log(number))

Calculating Pi to 10000 places (using the K.Takano equation (1982))
$ bc -l -q scale = 10000; (12*a(1/49)+32*a(1/57)-5*a(1/239)+12*a(1/110443))*4

A translated C function
Because the syntax of bc is very similar to that of C, published algorithms written in C can often be translated into BC quite easily, which immediately provides the arbitrary precision of BC. For example, in the Journal Of Statistical Software (July 2004, Volume 11, Issue 5), George Marsaglia published the following C code for the cumulative normal distribution:

With a few minutes of work, this can be translated to the following GNU bc code:

define normal(x) {    auto s,t,b,q,i,const; const=0.5*l(8*a(1)); s=x; t=0; b=x; q=x*x; i=1; while(s!=t) {s=(t=s)+(b*=q/(i+=2))}; return .5+s*e(-.5*q-const); }