Main Page: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
mNo edit summary
Replaced content with "<br><br>Here is more info on [http://ladyleem.hubweb.net/xe/?document_srl=53975 Single In Phoenix Arizona] review our web-page."
Line 1: Line 1:
{{About|specific activity radioactivity|the use in biochemistry|Enzyme assay#Specific activity}}
{{lowercase|title=dc (Unix)}}
{{multiple issues|
'''dc''' is a [[cross-platform]] [[reverse Polish notation|reverse-polish]] '''d'''esk '''c'''alculator which supports [[arbitrary-precision arithmetic]].<ref>{{man|1|dc||an arbitrary precision calculator}}</ref> It is one of the oldest [[Unix]] utilities, predating even the invention of the [[C (programming language)|C programming language]]; like other utilities of that vintage, it has a powerful set of features but an extremely terse syntax.<ref>{{cite web
{{technical|date=January 2014}}
|url=http://plan9.bell-labs.com/7thEdMan/vol2/dc
{{unreferenced|date=January 2014}}
|title=The sources for the manual page for 7th Edition Unix dc
}}
}}</ref><ref>{{cite web
|author=Ritchie, Dennis M.
|date=Sep. 1979
|url=http://cm.bell-labs.com/cm/cs/who/dmr/hist.html
|title=The Evolution of the Unix Timesharing System
}}</ref>
Traditionally, the more user-friendly (with its [[infix notation]]) [[bc programming language|bc]] calculator program was implemented on top of dc, although more modern implementations are related in the opposite fashion: dc uses bc's library for arithmetic.<ref>{{cite web
|url= http://directory.fsf.org/project/bc/
|title=Free Software Directory: BC
|accessdate=5 Jan 2009
}}</ref>


In [[nuclear sciences]] and technologies, "activity" is the SI quantity related to the phenomenon of natural and artificial [[radioactivity]]. The SI unit of "activity" is [[becquerel]] (Bq) while that of "'''specific activity'''" is Bq/g. The old unit of "activity" was [[curie]] (Ci) making "specific activity", Ci/g.
This article provides some examples in an attempt to give a general flavour of the language; for a complete list of commands and syntax, one should consult the [[man page]] for one's specific implementation.


== Half-life ==
==Basic operations==
To multiply four and five in dc (note that most of the whitespace is optional):


Experimentally-measured specific activity can be used to calculate the [[half-life]] of a radioactive element.
4 5 *
p


The definition of half-life ('''T<sub>1/2</sub>''') is that
This translates into "push four and five onto the stack, then, with the multiplication operator, pop two elements  from the stack, multiply them and push the result back on the stack."  Then the 'p' command is used to examine (print out to the screen) the top element on the stack.
the half life T<sub>1/2</sub> of an isotope is the length of time at which half
of a given quantity has decayed into another isotope, usually of a different element:
Or more generally:
Starting with '''''N<sub>0</sub>''''', atoms of an element, the number of atoms, '''''N''''',
remaining after time, '''''t''''', is given by:


: <math>N=N_0\left(\frac{1}{2}\right)^{t \over T_{1/2} }</math>
The [[arithmetic precision]] is changed with the command 'k', which sets the number of fractional digits (the number of digits following the [[radix point|point]]) to be used for arithmetic operations.  Since the default precision is zero, this sequence of commands produces '0' as a result:
2 3 / p


The natural log of both sides
By adjusting the precision with 'k', arbitrary number of decimal places can be produced.  This command sequence outputs '.66666'.


: <math>\ln(N)=\ln(N_0)+\left(\frac{t}{T_{1/2} }\right)\ln\left(\frac{1}{2}\right)</math>
5 k
2 3 / p


The derivative with respect to time, ''t''
To evaluate <math>\sqrt{(12 + (-3)^4)\over11}-22</math>: ('v' computes the square root of the top of the stack and '_' is used to input a negative number):


: <math>\frac{1}{N}\frac{dN}{dt}=\frac{\ln\left(\frac{1}{2}\right)}{T_{1/2} }</math>
12 _3 4 ^ + 11 / v 22 -
p


Multiplying both sides by ''N''
To swap the top two elements of the stack, use the 'r' command.  To duplicate the top element, use the 'd' command.


: <math>\frac{dN}{dt}=\frac{N\ln\left(\frac{1}{2}\right)}{T_{1/2} }</math>
==Input/Output==
To read a line from [[stdin]], use the '?' command.  This will evaluate the line as if it were a ''dc'' command, and so it is necessary that it be syntactically correct and potentially be a security problem since the '!' ''dc'' command will allow arbitrary command execution.


Yields
As mentioned above, 'p' will print the top of the stack with a newline after it.  'n' will pop the top of the stack and output it without a trailing newline.  'f' will dump the entire stack with one entry per line.


: <math>\frac{dN}{dt}=\frac{-0.693\,N}{T_{1/2} }</math>
''dc'' also support arbitrary input and output [[radix|radices]].  The 'i' command will pop the top of the stack and use it for the input base.  Hex digits must be in upper case to avoid collisions with ''dc'' commands and are not limited to A-F if the input radix is larger than 16.  The 'o' command does the same for the output base, but keep in mind that the input base will affect the parsing of every numeric value afterwards so it is usually advisable to set the output base first. To read the values, the 'K', 'I' and 'O' will push the current precision, input radix and output radix on to the top of the stack.


''dN''/''dt'' represents the decay rate of atoms. The negative sign shows that the rate is negative, so the number of atoms is decreasing with time. Rearranging terms:
As an example, to convert from hex to binary:


: <math>T_{1/2}=\frac{-0.693\,N}{\frac{dN}{dt} }</math>
16i2o DEADBEEFp


=== Example: half-life of Rb-87 ===
outputs <tt>11011110101011011011111011101111</tt>.


One gram of [[rubidium]]-87 and a radioactivity count rate that, after taking solid angle effects into account, is consistent with a decay rate of 3200 decays per second corresponds to a specific activity of {{val|3.2|e=6|u=Bq/kg}}. Rubidium's atomic weight is 87, so one gram is one 87th of a mole, or ''N''={{val|6.9|e=21}} atoms. Plugging in the numbers:
==Language Features==
===Registers===
In addition to these basic arithmetic and stack operations, dc includes support for [[Macro (computer science)|macros]], conditionals and storing of results for later retrieval.


: <math>T_{1/2}=\frac{-0.693(6.9\times 10^{21})}{-3200\text{ s}^{-1} }=1.5\times 10^{18}\text{ s or 47 billion years}</math>
The mechanism underlying macros and conditionals is the '''register''', which in dc is a storage location with a single character name which can be stored to and retrieved from: 'sc' pops the top of the stack and stores it in register c, and 'lc' pushes the value of register c onto the stack. For example:


== Formulation ==
3 sc 4 lc * p


Radioactivity is expressed as the decay rate of a particular [[radionuclide]] with decay constant ''λ'' and the number of atoms ''N'':
Registers can also be treated as secondary stacks, so values can be pushed and popped between them and the main stack using the 'S' and 'L' commands.


: <math>-\frac{dN}{dt}= \lambda N</math>
===Strings===
String values are enclosed in '[' and ']' characters and may be pushed on the stack and stored in registers.  The 'a' command will convert a the low order byte of the numeric value into an [[ASCII#ASCII_printable_characters|ASCII]] character, or if the top of the stack is a string it will replace it with the first character of the string.  There are no ways to build up strings or perform string manipulation other than executing it with the 'x' command, or printing it with the 'P' command.


Mass of the radionuclide is given by
The '#' character begins a comment to the end of the line.


: <math>\frac{N}{N_A} [\text{mol}] \times {m} [\text {g } \text{mol}^{-1}]</math>
===Macros===
Macros are then implemented by allowing registers and stack entries to be strings as well as numbers. A string can be printed, but it can also be executed (i.e. processed as a sequence of dc commands). So for instance we can store a macro to add one and then multiply by 2 into register m:


where ''m'' is [[mass number]] of the radionuclide and ''N<sub>A</sub>'' is [[Avogadro's constant]].
[1 + 2 *] sm


Specific radioactivity ''S'' is defined as radioactivity per unit mass of the radionuclide:
and then (using the 'x' command which executes the top of the stack) we can use it like this:


: <math>S [\text {Bq/g}] = \frac{\lambda N}{m N/N_A} = \frac{\lambda N_A}{m}</math>
3 lm x p


In addition, decay constant ''λ'' is related to the half-life ''T<sub>1/2</sub>'' by the following equation:
===Conditionals===
Finally, we can use this macro mechanism to provide conditionals. The command '=r' will pop two values from the stack, and execute the macro stored in register 'r' only if they are equal. So this will print the string 'equal' only if the top of the stack is equal to 5:
<pre>
[[equal]p] sm 5 =m
</pre>


: <math>{\lambda} = \frac{ln2}{T_{1/2}}</math>
Other conditionals are '>', '!>', '<', '!<', '!=', which will execute the specified macro if the top two values on the stack are greater, less than or equal to ("not greater"), less than,  greater than or equal to ("not less than"), and not equals, respectively.


Thus, specific radioactivity can also be described by
===Loops===
Looping is then possible by defining a macro which (conditionally) reinvokes itself.  A simple factorial of the top of the stack might be implemented as:


: <math>S = \frac{ln2 \times {N_A}}{T_{1/2} \times {m}}</math>
# F(x): return x!
# if x-1 > 1
#    return x * F(x-1)
# otherwise
#    return x
[d1-d1<F*]dsFxp


This equation is simplified by
The '1Q' command will exit from a macro, allowing an early return. 'q' will quit from two levels of macros (and ''dc'' itself if there are less than two levels on the call stack).  'z' will push the current stack depth before the 'z' operation.


: <math>S [\text {Bq/g}] \simeq \frac{4.17\times 10^{23}}{T_{1/2} [s]\times m}</math>
==Examples==
As an example of a relatively simple program in dc, this command (in 1 line):


When the unit of half-life converts a year
dc -e <nowiki>'[[Enter a number (metres), or 0 to exit]psj]sh[q]sz[lhx?d0=z10k39.370079*.5+0k12~1/rn[ feet ]
Pn[ inches]P10Pdx]dx'</nowiki>


: <math> S [\text {Bq/g}] = \frac{ln2 \times {N_A}}{T_{1/2} [s] \times {m}} =  \frac{ln2 \times {N_A}}{T_{1/2}[year] \times365\times24\times60\times60 \times m} \simeq \frac{1.32\times 10^{16} }{T_{1/2}[year] \times m}</math>
will convert distances from metres to feet and inches; the bulk of it is concerned with prompting for input, printing output in a suitable format and looping round to convert another number.


For example, specific radioactivity of [[Isotopes of radium|radium 226]] with a half-life of 1600 years is obtained by
As an example, here is an implementation of the [[Euclidean algorithm]] to find the [[Greatest common divisor|GCD]]:


: <math> \frac{1.32\times 10^{16} }{1600[year] \times 226} \simeq {3.7} \times 10^{10} [\text {Bq/g}] </math>
dc -e <nowiki>'??[dSarLa%d0<a]dsax+p' # shortest
dc -e '[a=]P?[b=]P?[dSarLa%d0<a]dsax+[GCD:]Pp' # easier-to-read version</nowiki>


This value derived from radium 226 was defined as unit of radioactivity known as [[Curie]] (Ci).
Computing the [[factorial]] of an input value, <math>n! = \prod_{i=1}^n i</math>
<pre>
dc -e '?[q]sQ[d1=Qd1-lFx*]dsFxp'
</pre>


[[Category:Units of radioactivity]]
A more complex example performs [[Diffie-Hellman key exchange]].  This was popular as a [[signature block]] among [[cypherpunk]]s during the [[ITAR]] debates<ref>{{cite web
|url=http://www.cypherspace.org/adam/rsa/perl-dh.html
|title=Diffie-Hellman in 2 lines of Perl
|accessdate=5 Jan 2009
|author=Adam Back
}}</ref>:
<pre>
#!/usr/bin/perl -- -export-a-crypto-system-sig Diffie-Hellman-2-lines
($g,$e,$m)=@ARGV,$m||die"$0 gen exp mod\n";print`echo "16dio1[d2%Sa2/d0<X+d
*La1=z\U$m%0]SX$e"[$g*]\EszlXx+p|dc`
</pre>
 
A commented version is slightly easier to understand and shows how to use loops, conditionals, and the 'q' command to return from a macro.  With a modern version of dc, the '|' command can be used to do arbitrary precision modular exponentiation without needing to write the X function.
<source lang=perl>
#!/usr/bin/perl
 
my ($g,$e,$m) = map { "\U$_" } @ARGV;
die "$0 gen exp mod\n" unless $m;
 
print `echo $g $e $m | dc -e '
# Hex input and output
16dio
# Read m, e and g from stdin on one line
?SmSeSg
 
# Function z: return g * top of stack
[lg*]sz
 
# Function Q: remove the top of the stack and return 1
[sb1q]sQ
 
# Function X(e): recursively compute g^e % m
# It is the same as Sm^Lm%, but handles arbitrarily large exponents.
# Stack at entry: e
# Stack at exit: g^e % m
# Since e may be very large, this uses the property that g^e % m ==
# if( e == 0 )
# return 1
# x = (g^(e/2)) ^ 2
# if( e % 2 == 1 )
# x *= g
# return x %
[
d 0=Q # return 1 if e==0 (otherwise, stack: e)
d 2% Sa # Store e%2 in a (stack: e)
2/ # compute e/2
lXx # call X(e/2)
d* # compute X(e/2)^2
La1=z # multiply by g if e%2==1
lm % # compute (g^e) % m
] SX
 
le # Load e from the register
lXx # compute g^e % m
p # Print the result
'`;
</source>
 
== See also ==
* [[bc programming language]]
* [[Calculator input methods]]
* [[HP calculators]]
* [[Orpie]], an RPN calculator
* [[Stack machine]]
 
==References==
{{reflist}}
 
==External links==
*Package [http://packages.debian.org/search?keywords=dc&searchon=names&exact=1&suite=all&section=all dc] in [[Debian GNU/Linux]] repositories
*[http://gnuwin32.sourceforge.net/packages/bc.htm Native Windows port] of ''[[bc programming language|bc]]'', which includes ''dc''.
 
{{unix commands}}
 
[[Category:cross-platform software]]
[[Category:Unix software]]
[[Category:free mathematics software]]
[[Category:Numerical programming languages]]
[[Category:stack-oriented programming languages]]
 
[[cs:Dc (programovací jazyk)]]
[[de:Dc (Unix)]]
[[fr:Dc (logiciel)]]
[[pl:Dc (informatyka)]]
[[ru:Dc]]
[[es:DC (Unix)]]

Revision as of 08:59, 12 August 2014

Template:Lowercase dc is a cross-platform reverse-polish desk calculator which supports arbitrary-precision arithmetic.[1] It is one of the oldest Unix utilities, predating even the invention of the C programming language; like other utilities of that vintage, it has a powerful set of features but an extremely terse syntax.[2][3] Traditionally, the more user-friendly (with its infix notation) bc calculator program was implemented on top of dc, although more modern implementations are related in the opposite fashion: dc uses bc's library for arithmetic.[4]

This article provides some examples in an attempt to give a general flavour of the language; for a complete list of commands and syntax, one should consult the man page for one's specific implementation.

Basic operations

To multiply four and five in dc (note that most of the whitespace is optional):

4 5 *
p

This translates into "push four and five onto the stack, then, with the multiplication operator, pop two elements from the stack, multiply them and push the result back on the stack." Then the 'p' command is used to examine (print out to the screen) the top element on the stack.

The arithmetic precision is changed with the command 'k', which sets the number of fractional digits (the number of digits following the point) to be used for arithmetic operations. Since the default precision is zero, this sequence of commands produces '0' as a result:

2 3 / p

By adjusting the precision with 'k', arbitrary number of decimal places can be produced. This command sequence outputs '.66666'.

5 k
2 3 / p

To evaluate : ('v' computes the square root of the top of the stack and '_' is used to input a negative number):

12 _3 4 ^ + 11 / v 22 -
p

To swap the top two elements of the stack, use the 'r' command. To duplicate the top element, use the 'd' command.

Input/Output

To read a line from stdin, use the '?' command. This will evaluate the line as if it were a dc command, and so it is necessary that it be syntactically correct and potentially be a security problem since the '!' dc command will allow arbitrary command execution.

As mentioned above, 'p' will print the top of the stack with a newline after it. 'n' will pop the top of the stack and output it without a trailing newline. 'f' will dump the entire stack with one entry per line.

dc also support arbitrary input and output radices. The 'i' command will pop the top of the stack and use it for the input base. Hex digits must be in upper case to avoid collisions with dc commands and are not limited to A-F if the input radix is larger than 16. The 'o' command does the same for the output base, but keep in mind that the input base will affect the parsing of every numeric value afterwards so it is usually advisable to set the output base first. To read the values, the 'K', 'I' and 'O' will push the current precision, input radix and output radix on to the top of the stack.

As an example, to convert from hex to binary:

16i2o DEADBEEFp

outputs 11011110101011011011111011101111.

Language Features

Registers

In addition to these basic arithmetic and stack operations, dc includes support for macros, conditionals and storing of results for later retrieval.

The mechanism underlying macros and conditionals is the register, which in dc is a storage location with a single character name which can be stored to and retrieved from: 'sc' pops the top of the stack and stores it in register c, and 'lc' pushes the value of register c onto the stack. For example:

3 sc 4 lc * p

Registers can also be treated as secondary stacks, so values can be pushed and popped between them and the main stack using the 'S' and 'L' commands.

Strings

String values are enclosed in '[' and ']' characters and may be pushed on the stack and stored in registers. The 'a' command will convert a the low order byte of the numeric value into an ASCII character, or if the top of the stack is a string it will replace it with the first character of the string. There are no ways to build up strings or perform string manipulation other than executing it with the 'x' command, or printing it with the 'P' command.

The '#' character begins a comment to the end of the line.

Macros

Macros are then implemented by allowing registers and stack entries to be strings as well as numbers. A string can be printed, but it can also be executed (i.e. processed as a sequence of dc commands). So for instance we can store a macro to add one and then multiply by 2 into register m:

[1 + 2 *] sm

and then (using the 'x' command which executes the top of the stack) we can use it like this:

3 lm x p

Conditionals

Finally, we can use this macro mechanism to provide conditionals. The command '=r' will pop two values from the stack, and execute the macro stored in register 'r' only if they are equal. So this will print the string 'equal' only if the top of the stack is equal to 5:

[[equal]p] sm 5 =m

Other conditionals are '>', '!>', '<', '!<', '!=', which will execute the specified macro if the top two values on the stack are greater, less than or equal to ("not greater"), less than, greater than or equal to ("not less than"), and not equals, respectively.

Loops

Looping is then possible by defining a macro which (conditionally) reinvokes itself. A simple factorial of the top of the stack might be implemented as:

# F(x): return x!
# if x-1 > 1
#    return x * F(x-1)
# otherwise
#    return x
[d1-d1<F*]dsFxp

The '1Q' command will exit from a macro, allowing an early return. 'q' will quit from two levels of macros (and dc itself if there are less than two levels on the call stack). 'z' will push the current stack depth before the 'z' operation.

Examples

As an example of a relatively simple program in dc, this command (in 1 line):

dc -e '[[Enter a number (metres), or 0 to exit]psj]sh[q]sz[lhx?d0=z10k39.370079*.5+0k12~1/rn[ feet ]
Pn[ inches]P10Pdx]dx'

will convert distances from metres to feet and inches; the bulk of it is concerned with prompting for input, printing output in a suitable format and looping round to convert another number.

As an example, here is an implementation of the Euclidean algorithm to find the GCD:

dc -e '??[dSarLa%d0<a]dsax+p' # shortest
dc -e '[a=]P?[b=]P?[dSarLa%d0<a]dsax+[GCD:]Pp' # easier-to-read version

Computing the factorial of an input value,

dc -e '?[q]sQ[d1=Qd1-lFx*]dsFxp'

A more complex example performs Diffie-Hellman key exchange. This was popular as a signature block among cypherpunks during the ITAR debates[5]:

#!/usr/bin/perl -- -export-a-crypto-system-sig Diffie-Hellman-2-lines
($g,$e,$m)=@ARGV,$m||die"$0 gen exp mod\n";print`echo "16dio1[d2%Sa2/d0<X+d
*La1=z\U$m%0]SX$e"[$g*]\EszlXx+p|dc`

A commented version is slightly easier to understand and shows how to use loops, conditionals, and the 'q' command to return from a macro. With a modern version of dc, the '|' command can be used to do arbitrary precision modular exponentiation without needing to write the X function.

#!/usr/bin/perl

my ($g,$e,$m) = map { "\U$_" } @ARGV;
die "$0 gen exp mod\n" unless $m;

print `echo $g $e $m | dc -e '
# Hex input and output
16dio
# Read m, e and g from stdin on one line
?SmSeSg

# Function z: return g * top of stack
[lg*]sz

# Function Q: remove the top of the stack and return 1
[sb1q]sQ

# Function X(e): recursively compute g^e % m
# It is the same as Sm^Lm%, but handles arbitrarily large exponents.
# Stack at entry: e
# Stack at exit: g^e % m
# Since e may be very large, this uses the property that g^e % m == 
#	if( e == 0 )
#		return 1
#	x = (g^(e/2)) ^ 2
#	if( e % 2 == 1 )
#		x *= g
#	return x %
[
	d 0=Q		# return 1 if e==0 (otherwise, stack: e)
	d 2% Sa		# Store e%2 in a (stack: e)
	2/		# compute e/2
	lXx		# call X(e/2)
	d*		# compute X(e/2)^2
	La1=z		# multiply by g if e%2==1
	lm %		# compute (g^e) % m
] SX

le	# Load e from the register
lXx	# compute g^e % m
p	# Print the result
'`;

See also

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

External links

Template:Unix commands

cs:Dc (programovací jazyk) de:Dc (Unix) fr:Dc (logiciel) pl:Dc (informatyka) ru:Dc es:DC (Unix)