MPAライブラリ 使用例1



[ 簡単な説明 ]

MPAライブラリ の機能全般をテストします。
出力結果と突き合わせて、各種関数の使用方法を見てください。


プログラム・ソース("mpa.c")           top (先頭に戻る)
/*		mpa.c		Multiple-Precision Arithmetic	*/
#include	"mpa.h"

int main(void)
{
	MPA v, w, x, y, z;
	int k;
	ULONG i = 0;
	time_t now;

	printf("m_set_a()   2.0 --> x\n");
	x = m_set_a("2.0");	m_print("x = ", x, 1);

	printf("m_set_a()   0.11 --> y\n");
	y = m_set_a("0.11");	m_print("y = ", y, 1);
	printf("m_prt_b() exec.\n");
							m_prt_b("y = ", y, 1);
	printf("m_prt_m() exec.\n");
							m_prt_m("y = ", y, 1);

	printf("m_sqrt()    sqrt(x) --> z\n");
	z = m_sqrt(x);		m_print("z = sqrt(x) = ", z, 1);

	printf("m_sqr()     sqr(z) --> x\n");
	x = m_sqr(z);		m_print("x = z * z = ", x, 1);

	printf("m_pwr_s()   x^3 --> z\n");
	z = m_pwr_s(x, 3L);	m_print("z = x ^ 3 = ", z, 1);

	printf("m_div1()    x/y --> x\n");
	m_div1(&x, y);		m_print("x = x / y = ", x, 1);

	printf("m_div1()    y/x --> y\n");
	m_div1(&y, x);		m_print("y = y / x = ", y, 1);

	printf("m_div1_s()  x/2 --> x\n");
	m_div1_s(&x, 2);	m_print("x = x / 2 = ", x, 1);

	printf("m_div1_s()  y/-3 --> y\n");
	m_div1_s(&y, -3);	m_print("y = y / -3 = ", y, 1);

	printf("m_add()     x+y --> z\n");
	z = m_add(x, y);	m_print("z = x + y = ", z, 1);

	printf("m_add1()    z+x --> z\n");
	m_add1(&z, x);		m_print("z = z + x = ", z, 1);

	printf("m_sub1()    z-x --> z\n");
	m_sub1(&z, x);		m_print("z = z - x = ", z, 1);

	printf("m_mul()     x*y --> z\n");
	z = m_mul(x, y);	m_print("z = x * y = ", z, 1);

	printf("m_mul1()    z*x --> x\n");
	m_mul1(&z, x);		m_print("z = z * x = ", z, 1);

	printf("m_set_a()   1.2 --> z\n");
	z = m_set_a("1.2");	m_print("z = ", z, 1);

	printf("m_inv()     1/z --> z\n");
	z = m_inv(z);		m_print("z = 1 / z = ", z, 1);

	printf("m_kaijo()   10! --> x\n");
	x = m_kaijo(10);	m_print("x = 10! = ", x, 1);

	printf("m_iset_s()  x --> (long)\n");
						printf("x = 10! = %ld\n", m_iset_s(x));

	printf("m_kaijo()   15! --> y\n");
	y = m_kaijo(15);	m_print("y = 15! = ", y, 1);

	printf("m_gcd()     gcd(x,y) --> z\n");
	z = m_gcd(x, y);	m_print("z = gcd(x, y) = ", z, 1);

	printf("m_lcm()     lcm(x,y) --> z\n");
	z = m_lcm(x, y);	m_print("z = lcm(x, y) = ", z, 1);

	printf("m_cmp() exec.\n");
	printf("(x : y) = %d    (y : z) = %d    (z : x) = %d\n",
		m_cmp(x, y), m_cmp(y, z), m_cmp(z, x)); 

	printf("m_iset(x)   x,y,z --> (double)\n");
	printf("x = %e, y = %e, z = %e\n", m_iset(x), m_iset(y), m_iset(z));

	m_print("pi()=", pi(), 0);
	m_print("_PI = ", _PI, 0);

	m_print("e()=", e(), 0);
	m_print("_E =", _E, 0);
	m_print("_2_PI", _2_PI, 0);
	m_print("_PI_2", _PI_2, 0);
	m_print("_PI_4", _PI_4, 0);

	x = _PI_4;
	printf("m_sin()     sin(π/4) --> y\n");
	y = m_sin(x);		m_print("sin(pi/4) = ", y, 0);

	printf("m_cos()     cos(π/4) --> z\n");
	z = m_cos(x);		m_print("cos(pi/4) = ", z, 0);

	printf("m_set_a()   0.9999 --> x\n");
	x = m_set_a("0.9999");	m_print("x = ", x, 1);

	printf("m_asin()    arc sin(x) --> z\n");
	z = m_asin(x);		m_print("asin(x) = ", z, 1);

	printf("m_acos()    arc cos(x) --> y\n");
	y = m_acos(x);		m_print("acos(x) = ", y, 1);

	printf("m_set_a()   0.0111 --> x\n");
	x = m_set_a("0.0111");	m_print("x = ", x, 1);

	printf("m_asin()    arc sin(x) --> z\n");
	z = m_asin(x);		m_print("asin(x) = ", z, 1);

	printf("m_acos()    arc cos(x) --> y\n");
	y = m_acos(x);		m_print("acos(x) = ", y, 1);

	printf("m_set_a()   -0.9999 --> x\n");
	x = m_set_a("-0.9999");	m_print("x = ", x, 1);

	printf("m_asin()    arc sin(x) --> z\n");
	z = m_asin(x);		m_print("asin(x) = ", z, 1);

	printf("m_acos()    arc cos(x) --> y\n");
	y = m_acos(x);		m_print("acos(x) = ", y, 1);

	printf("m_set_a()   -0.0111 --> x\n");
	x = m_set_a("-0.0111");	m_print("x = ", x, 1);

	printf("m_asin()    arc sin(x) --> z\n");
	z = m_asin(x);		m_print("asin(x) = ", z, 1);

	printf("m_acos()    arc cos(x) --> y\n");
	y = m_acos(x);		m_print("acos(x) = ", y, 1);

	printf("m_set_a()   123.456 --> x\n");
	x = m_set_a("123.456");	m_print("x = ", x, 1);

	printf("m_log()     ln(x) --> z\n");
	z = m_log(x);		m_print("log(x) = ", z, 1);

	printf("m_log10()   log10(x) --> y\n");
	y = m_log10(x);		m_print("log10(x) = ", y, 1);

	printf("m_set_a()   2 --> x\n");
	x = m_set_a("2");	m_print("x = ", x, 0);

	printf("m_x1n()     x^(1/3) --> z\n");
	z = m_x1n(x, 3);	m_print("z = 2 ^ (1/3) = ", z, 0);

	printf("m_sqr()     z^2 --> y\n");
	y = m_sqr(z);		m_print("y = z ^ 2 = ", y, 0);

	printf("m_mul1()	y*z --> y\n");
	m_mul1(&y, z);		m_print("z ^ 3 = ", y, 0);

	printf("m_set_a()   3 --> x\n");
	x = m_set_a("3");	m_print("x = ", x, 0);

	printf("m_x_y()     z^x --> y\n");
	y = m_x_y(z, x);	m_print("z ^ 3 = ", y, 0);

	printf("m_set_a()   3.2879 --> y\n");
	y = m_set_a("3.2879");	m_print("y = ", y, 0);

	printf("m_x_y()     x^y --> z\n");
	z = m_x_y(x, y);	m_print("x ^ y = ", z, 0);

	printf("m_exp()     e^x --> z\n");
	z = m_exp(x);		m_print("exp(x) = ", z, 1);

	printf("m_10()      10^x --> z\n");
	z = m_10(x);		m_print("10^(x) = ", z, 1);

	printf("m_set_a()   1.5 --> x\n");
	x = m_set_a("1.5");	m_print("x = ", x, 0);

	printf("m_hcos()    hcos(x) --> y\n");
	y = m_hcos(x);		m_print("hcos(x) = ", y, 1);

	printf("m_hsin()    hsin(x) --> y\n");
	y = m_hsin(x);		m_print("hsin(x) = ", y, 1);

	printf("m_htan()    htan(x) --> y\n");
	y = m_htan(x);		m_print("htan(x) = ", y, 1);

	printf("m_bino()    5 C 3 --> x\n");
	x = m_bino( 5,  3);	m_print(" 5 C 3  = ", x, 0);

	printf("m_bino()    10 C 5 --> x\n");
	x = m_bino(10,  5);	m_print("10 C 5  = ", x, 0);

	printf("m_bino()    50 C 30 --> x\n");
	x = m_bino(50, 30);	m_print("50 C 30 = ", x, 0);

	printf("m_set_a()   13783.1721389 --> x\n");
	x = m_set_a("13783.1721389");	m_print("x = ", x, 0);

	printf("m_round()   round(x,3) --> y\n");
	y = m_round(x, 3);	m_print("round(x, 3) = ", y, 0);

	printf("m_round()   round(x,0) --> y\n");
	y = m_round(x, 0);	m_print("round(x, 0) = ", y, 0);

	printf("m_round()   round(x,-3) --> y\n");
	y = m_round(x, -3);	m_print("round(x, -3) = ", y, 0);

	printf("m_srand() exec.\n");
	m_srand(time(&now));
	y = m_set_a("0.0001");
	for(k = 0; k < 10; k++)
	{
		i = 0;
		do
		{
			x = m_rand();
			i++;
		} while(m_cmp(x, y) >= 0);
		printf("count = %ld (target=10000) : ", i);
		m_print("r = ", x, 1);
	}

	printf("m_set_a()   3 --> x\n");
	x = m_set_a("3");	m_print("x = ", x, 0);

	printf("m_log()     ln(x) --> y\n");
	y = m_log(x);		m_print("ln(3)=", y, 0);

	printf("ln2()       ln(2) --> x\n");
	x = ln2();			m_print("ln(2)=", x, 0);

	printf("ln625()     ln(0.625) --> y\n");
	y = ln625();		m_print("ln(.625)=", y, 0);

	printf("m_mul_s()   x*4 --> z\n");
	z = m_mul_s(x, 4);	m_print("4*ln(2) = ", z, 0);

	printf("m_add1()    y+z --> y\n");
	m_add1(&y, z);		m_print("ln(10)=", y, 0);

	printf("m_div()     x/y --> z\n");
	z = m_div(x, y);	m_print("log10(2)=", z, 0);

	printf("ln75()      ln(0.75) --> z\n");
	z = ln75();			m_print("ln(.75)=", z, 0);

	printf("m_mul1_s()  x*2 --> x\n");
	m_mul1_s(&x, 2);	m_print("x = ", x, 0);

	printf("m_add1()    z+x --> z\n");
	m_add1(&z, x);		m_print("ln(3)=", z, 0);

	printf("m_div()     z/y --> x\n");
	x = m_div(z, y);	m_print("log10(3)=", x, 0);

	printf("π --> x\n");
	x = _PI;			m_print("x = ", x, 0);

	printf("m_div1_s()  x/3 --> x\n");
	m_div1_s(&x, 3);	m_print("x = ", x, 0);

	printf("m_tan()     tan(x) --> y\n");
	y = m_tan(x);		m_print("tan(x) =", y, 0);

	m_print("sqrt(3)=", _SQRT3, 0);

	printf("m_set_a()   0.5 --> x\n");
	x = m_set_a("0.5");	m_print("x = ", x, 0);

	printf("m_acos()    arc cos(x) --> y\n");
	y = m_acos(x);		m_print("acos(x)=", y, 0);

	printf("sqrt(3) --> x\n");
	x = _SQRT3;			m_print("x = ", x, 0);

	printf("m_div_s()   x/2 --> v\n");
	v = m_div_s(x, 2);	m_print("v = ", v, 0);

	printf("m_asin()    arc sin(v) --> y\n");
	y = m_asin(v);		m_print("asin(x)=", y, 0);

	printf("m_atan()    arc tan(x) --> y\n");
	y = m_atan(x);		m_print("atan(x)=", y, 0);

	printf("π --> z\n");
	z = _PI;			m_print("z = ", z, 0);

	printf("m_div1_s()  z/3 --> z\n");
	m_div1_s(&z, 3);	m_print("pi / 3 =", z, 0);

	return 1;
}

出力結果           top (先頭に戻る)
m_set_a()   2.0 --> x
x =  +2.
m_set_a()   0.11 --> y
y =  +0.1100 
m_prt_b() exec.
y = 
 zero = 0   sign = 1   exp = -1
000111000010100 011110101110000 101000111101011 100001010001111 010111000010100 011110101110000 101000111101011 100001010001111 010111000010100 011110101110000 
m_prt_m() exec.
y =  : <   -1> + 3604. 15728 20971 17039 11796 15728 20971 17039 11796 15728 20971
m_sqrt()    sqrt(x) --> z
z = sqrt(x) =  +1.4142 1356 2373 0950 4880 1688 7242 0969 8078 5696 7188 
m_sqr()     sqr(z) --> x
x = z * z =  +2.
m_pwr_s()   x^3 --> z
z = x ^ 3 =  +8.
m_div1()    x/y --> x
x = x / y =  +18.1818 1818 1818 1818 1818 1818 1818 1818 1818 1818 1818 
m_div1()    y/x --> y
y = y / x =  +0.0060 5000 
m_div1_s()  x/2 --> x
x = x / 2 =  +9.0909 0909 0909 0909 0909 0909 0909 0909 0909 0909 0909 
m_div1_s()  y/-3 --> y
y = y / -3 =  -0.0020 1666 6666 6666 6666 6666 6666 6666 6666 6666 6666 6667 
m_add()     x+y --> z
z = x + y =  +9.0888 9242 4242 4242 4242 4242 4242 4242 4242 4242 4242 
m_add1()    z+x --> z
z = z + x =  +18.1798 0151 5151 5151 5151 5151 5151 5151 5151 5151 5152 
m_sub1()    z-x --> z
z = z - x =  +9.0888 9242 4242 4242 4242 4242 4242 4242 4242 4242 4242 
m_mul()     x*y --> z
z = x * y =  -0.0183 3333 3333 3333 3333 3333 3333 3333 3333 3333 3333 3333 
m_mul1()    z*x --> x
z = z * x =  -0.1666 6666 6666 6666 6666 6666 6666 6666 6666 6666 6666 6667 
m_set_a()   1.2 --> z
z =  +1.2000 
m_inv()     1/z --> z
z = 1 / z =  +0.8333 3333 3333 3333 3333 3333 3333 3333 3333 3333 3333 3333 
m_kaijo()   10! --> x
x = 10! =  +362 8800.
m_iset_s()  x --> (long)
x = 10! = 3628800
m_kaijo()   15! --> y
y = 15! =  +1 3076 7436 8000.
m_gcd()     gcd(x,y) --> z
z = gcd(x, y) =  +362 8800.
m_lcm()     lcm(x,y) --> z
z = lcm(x, y) =  +1 3076 7436 8000.
m_cmp() exec.
(x : y) = -1    (y : z) = 0    (z : x) = 1
m_iset(x)   x,y,z --> (double)
x = 3.628800e+06, y = 1.307674e+12, z = 1.307674e+12
pi()= +3.1415 9265 3589 7932 3846 2643 3832 7950 2884 1971 6939 9375 1058 
2097 4944 5923 0781 6406 2862 0899 8628 0348 2534 2117 0679 8214 8086 5132 
8230 6647 0938 4460 9550 5822 3172 5359 4081 2848 1117 4502 8410 2701 9385 
2110 5559 6446 2294 8954 9303 8196 4428 8109 7566 5933 4461 2847 5648 2337 
8678 3165 2712 0190 9145 6485 6692 3460 3486 1045 4326 6482 1339 3607 2602 
4914 1273 7245 8700 6606 3155 8817 4881 5209 2096 2829 2540 9171 5364 3678 
9259 0360 0113 3053 0548 8204 6652 1384 1469 5194 1511 6094 3305 7270 3657 
5959 1953 0921 8611 7381 9326 1179 3105 1185 
_PI =  +3.1415 9265 3589 7932 3846 2643 3832 7950 2884 1971 6939 9375 
1058 2097 4944 5923 0781 6406 2862 0899 8628 0348 2534 2117 0679 8214 8086 
5132 8230 6647 0938 4460 9550 5822 3172 5359 4081 2848 1117 4502 8410 2701 
9385 2110 5559 6446 2294 8954 9303 8196 4428 8109 7566 5933 4461 2847 5648 
2337 8678 3165 2712 0190 9145 6485 6692 3460 3486 1045 4326 6482 1339 3607 
2602 4914 1273 7245 8700 6606 3155 8817 4881 5209 2096 2829 2540 9171 5364 
3678 9259 0360 0113 3053 0548 8204 6652 1384 1469 5194 1511 6094 3305 7270 
3657 5959 1953 0921 8611 7381 9326 1179 3105 1185 
e()= +2.7182 8182 8459 0452 3536 0287 4713 5266 2497 7572 4709 3699 9595 
7496 6967 6277 2407 6630 3535 4759 4571 3821 7852 5166 4274 2746 6391 9320 
0305 9921 8174 1359 6629 0435 7290 0334 2952 6059 5630 7381 3232 8627 9434 
9076 3233 8298 8075 3195 2510 1901 1573 8341 8793 0702 1540 8914 9934 8841 
6750 9244 7614 6066 8082 2648 0016 8477 4118 5374 2345 4424 3710 7539 0777 
4499 2069 5517 0276 1838 6062 6133 1384 5830 0075 2044 9338 2656 0297 6067 
3711 3200 7093 2870 9127 4437 4704 7230 6969 7720 9310 1416 9283 6819 0255 
1510 8657 4637 7211 1252 3897 8442 5056 9537 
_E = +2.7182 8182 8459 0452 3536 0287 4713 5266 2497 7572 4709 3699 9595 
7496 6967 6277 2407 6630 3535 4759 4571 3821 7852 5166 4274 2746 6391 9320 
0305 9921 8174 1359 6629 0435 7290 0334 2952 6059 5630 7381 3232 8627 9434 
9076 3233 8298 8075 3195 2510 1901 1573 8341 8793 0702 1540 8914 9934 8841 
6750 9244 7614 6066 8082 2648 0016 8477 4118 5374 2345 4424 3710 7539 0777 
4499 2069 5517 0276 1838 6062 6133 1384 5830 0075 2044 9338 2656 0297 6067 
3711 3200 7093 2870 9127 4437 4704 7230 6969 7720 9310 1416 9283 6819 0255 
1510 8657 4637 7211 1252 3897 8442 5056 9537 
_2_PI +6.2831 8530 7179 5864 7692 5286 7665 5900 5768 3943 3879 8750 2116 
4194 9889 1846 1563 2812 5724 1799 7256 0696 5068 4234 1359 6429 6173 0265 
6461 3294 1876 8921 9101 1644 6345 0718 8162 5696 2234 9005 6820 5403 8770 
4221 1119 2892 4589 7909 8607 6392 8857 6219 5133 1866 8922 5695 1296 4675 
7356 6330 5424 0381 8291 2971 3384 6920 6972 2090 8653 2964 2678 7214 5204 
9828 2547 4491 7401 3212 6311 7634 9763 0418 4192 5658 5081 8343 0728 7357 
8518 0720 0226 6106 1097 6409 3304 2768 2939 0388 3023 2188 6611 4540 7315 
1918 3906 1843 7223 4763 8652 2358 6210 2371 
_PI_2 +1.5707 9632 6794 8966 1923 1321 6916 3975 1442 0985 8469 9687 5529 
1048 7472 2961 5390 8203 1431 0449 9314 0174 1267 1058 5339 9107 4043 2566 
4115 3323 5469 2230 4775 2911 1586 2679 7040 6424 0558 7251 4205 1350 9692 
6055 2779 8223 1147 4477 4651 9098 2214 4054 8783 2966 7230 6423 7824 1168 
9339 1582 6356 0095 4572 8242 8346 1730 1743 0522 7163 3241 0669 6803 6301 
2457 0636 8622 9350 3303 1577 9408 7440 7604 6048 1414 6270 4585 7682 1839 
4629 5180 0056 6526 5274 4102 3326 0692 0734 7597 0755 8047 1652 8635 1828 
7979 5976 5460 9305 8690 9663 0589 6552 5593 
_PI_4 +0.7853 9816 3397 4483 0961 5660 8458 1987 5721 0492 9234 9843 7764 5524 
3736 1480 7695 4101 5715 5224 9657 0087 0633 5529 2669 9553 7021 6283 2057 
6661 7734 6115 2387 6455 5793 1339 8520 3212 0279 3625 7102 5675 4846 3027 
6389 9111 5573 7238 7325 9549 1107 2027 4391 6483 3615 3211 8912 0584 4669 
5791 3178 0047 7286 4121 4173 0865 0871 5261 3581 6620 5334 8401 8150 6228 
5318 4311 4675 1651 5788 9704 3720 3802 3024 0707 3135 2292 8841 0919 7314 
7590 0028 3263 2637 2051 1663 0346 0367 3798 5377 9023 5826 4317 5914 3989 
7988 2730 4652 9345 4831 5294 8276 2796 3699 
m_sin()     sin(π/4) --> y
sin(pi/4) =  +0.7071 0678 1186 5475 2440 0844 3621 0484 9039 2848 3593 7688 
4740 3658 8339 8689 9536 6239 2310 5351 9425 1937 6716 3820 7863 6750 6923 
1154 5614 8512 4624 1802 7925 3686 0632 2060 7485 4996 7915 7066 1133 2963 
7527 9637 7899 9752 5057 6391 0302 8573 5054 7799 8580 2985 1372 6729 8431 
0073 6425 8709 3204 4459 9304 7761 6461 5242 1543 5716 0725 4198 8130 1813 
9976 2570 3994 8436 2669 8273 1659 0441 4820 3103 0762 9176 1975 2737 2875 
1438 7998 0864 9177 8761 0168 7659 2850 5677 1873 0170 4249 4235 8019 3449 
9853 4950 2407 5152 7201 3895 1582 2712 3911 5342 4645 
m_cos()     cos(π/4) --> z
cos(pi/4) =  +0.7071 0678 1186 5475 2440 0844 3621 0484 9039 2848 3593 7688 
4740 3658 8339 8689 9536 6239 2310 5351 9425 1937 6716 3820 7863 6750 6923 
1154 5614 8512 4624 1802 7925 3686 0632 2060 7485 4996 7915 7066 1133 2963 
7527 9637 7899 9752 5057 6391 0302 8573 5054 7799 8580 2985 1372 6729 8431 
0073 6425 8709 3204 4459 9304 7761 6461 5242 1543 5716 0725 4198 8130 1813 
9976 2570 3994 8436 2669 8273 1659 0441 4820 3103 0762 9176 1975 2737 2875 
1438 7998 0864 9177 8761 0168 7659 2850 5677 1873 0170 4249 4235 8019 3449 
9853 4950 2407 5152 7201 3895 1582 2712 3911 5342 4648 
m_set_a()   0.9999 --> x
x =  +0.9999 
m_asin()    arc sin(x) --> z
asin(x) =  +1.5566 5407 3274 9542 1988 8699 6104 0463 8489 0353 9905 
m_acos()    arc cos(x) --> y
acos(x) =  +0.0141 4225 3519 9423 9934 2622 0812 3511 2953 0631 8564 5306 
m_set_a()   0.0111 --> x
x =  +0.0111 
m_asin()    arc sin(x) --> z
asin(x) =  +0.0111 0022 7963 7783 4419 4453 4654 8042 6505 6997 8567 9800 
m_acos()    arc cos(x) --> y
acos(x) =  +1.5596 9609 8831 1182 7503 6868 2261 5932 4936 3987 9902 
m_set_a()   -0.9999 --> x
x =  -0.9999 
m_asin()    arc sin(x) --> z
asin(x) =  -1.5566 5407 3274 9542 1988 8699 6104 0463 8489 0353 9905 
m_acos()    arc cos(x) --> y
acos(x) =  +3.1274 5040 0069 8508 3912 0021 3020 4438 9931 1339 8375 
m_set_a()   -0.0111 --> x
x =  -0.0111 
m_asin()    arc sin(x) --> z
asin(x) =  -0.0111 0022 7963 7783 4419 4453 4654 8042 6505 6997 8567 9800 
m_acos()    arc cos(x) --> y
acos(x) =  +1.5818 9655 4758 6749 6342 5775 1571 2017 7947 7983 7038 
m_set_a()   123.456 --> x
x =  +123.4560 
m_log()     ln(x) --> z
log(x) =  +4.8158 8481 7283 2638 8310 9232 1051 6652 5577 1721 5814 
m_log10()   log10(x) --> y
log10(x) =  +2.0915 1220 1627 7716 8106 9399 7770 6790 5794 6535 8419 
m_set_a()   2 --> x
x =  +2.
m_x1n()     x^(1/3) --> z
z = 2 ^ (1/3) =  +1.2599 2104 9894 8731 6476 7210 6072 7822 8350 5702 
5146 4701 5079 8008 1975 1121 5529 9676 5139 5948 3729 3965 6243 6255 0941 
5431 0256 0356 1566 5259 3990 2404 0613 7372 2845 9110 3042 6935 5246 9606 
4261 6625 0009 7747 4526 5654 8030 6867 1854 0551 8689 2458 7251 6764 1993 
7370 9695 0983 8278 3161 3991 5512 9313 6953 6618 3947 4634 4857 6570 3031 
1909 5895 9847 4110 5981 1629 0705 3590 8164 7801 1473 5213 2548 4771 2978 
8024 2208 5820 5325 7972 5266 6220 2669 0056 6560 8199 4715 6281 7640 5060 
6648 2677 3572 6704 1948 6207 6214 4296 5694 2050 7931 9172 
m_sqr()     z^2 --> y
y = z ^ 2 =  +1.5874 0105 1968 1994 7475 1705 6392 7230 8260 3914 9332 
7899 8530 0980 8285 7618 2521 6505 6242 1917 3273 5442 1326 2220 9570 2293 
4761 6813 2201 7903 4976 5989 8152 7522 7814 0011 0445 4144 6619 3755 1827 
8562 4368 7329 0512 4850 7229 2337 4497 7427 5364 6071 6753 7106 6638 1847 
5711 7523 5157 5061 7045 6027 2792 4265 0876 7815 7144 6940 6950 9962 4484 
5050 9638 7787 0022 3172 9232 2609 4249 6708 4796 6044 9914 1550 8010 9940 
1078 5308 1210 3380 1245 0308 5826 8253 6745 9016 7176 3948 8477 2523 6020 
8912 9549 1732 2106 5479 8947 8852 5868 4912 0268 8445 
m_mul1()	y*z --> y
z ^ 3 =  +2.
m_set_a()   3 --> x
x =  +3.
m_x_y()     z^x --> y
z ^ 3 =  +2.
m_set_a()   3.2879 --> y
y =  +3.2879 
m_x_y()     x^y --> z
x ^ y =  +37.0447 7602 8209 5410 5268 7108 8990 9712 7649 4117 1616 6763 
2474 6487 8342 4696 2965 0143 4008 2125 8230 7551 4705 5520 6619 4958 8197 
6887 4286 8790 9418 3266 6068 7344 5328 0097 6532 9193 6671 7616 2739 1047 
7574 9121 5843 0907 3849 4669 6115 2716 7750 7468 0170 7063 1476 3989 2591 
1116 6451 4360 2355 9631 4874 1226 9238 1453 4722 7057 5134 7358 2629 9562 
0503 9357 9182 5829 4120 9652 3815 7440 0577 6344 7036 5321 4818 1578 1700 
0945 7217 3351 9573 6302 3420 3604 2542 4408 5942 9282 7561 8067 2267 3999 
6546 2497 0576 0600 8739 8339 2282 4162 0509 8347 
m_exp()     e^x --> z
exp(x) =  +20.0855 3692 3187 6677 4092 8529 6545 8171 7896 9879 0784 
m_10()      10^x --> z
10^(x) =  +1000.
m_set_a()   1.5 --> x
x =  +1.5000 
m_hcos()    hcos(x) --> y
hcos(x) =  +2.3524 0961 5243 2473 2576 7667 9654 4164 4170 1739 6075 
m_hsin()    hsin(x) --> y
hsin(x) =  +2.1292 7945 5094 8174 9683 4387 4946 7763 1648 8317 8912 
m_htan()    htan(x) --> y
htan(x) =  +0.9051 4825 3644 8664 3824 2303 6964 5649 5597 2276 4113 5159 
m_bino()    5 C 3 --> x
 5 C 3  =  +10.
m_bino()    10 C 5 --> x
10 C 5  =  +252.
m_bino()    50 C 30 --> x
50 C 30 =  +47 1292 1224 3960.
m_set_a()   13783.1721389 --> x
x =  +1 3783.1721 3890 
m_round()   round(x,3) --> y
round(x, 3) =  +1 4000.
m_round()   round(x,0) --> y
round(x, 0) =  +1 3783.
m_round()   round(x,-3) --> y
round(x, -3) =  +1 3783.1720 
m_srand() exec.
count = 51668 (target=10000) : r =  +0.0000 3677 1498 1166 4982 9913 0615 9367 6224 1784 9841 5734 
count = 5624 (target=10000) : r =  +0.0000 0000 0000 0000 0000 0007 4913 7033 9371 9935 0333 3250 5688 8686 
1611 6429 0377 
count = 31634 (target=10000) : r =  +0.0000 6634 2050 6157 3049 5223 8277 8207 3807 7181 4870 3952 
count = 795 (target=10000) : r =  +0.0000 6778 2663 8479 9261 3968 4218 7070 5177 4123 9270 7145 
count = 7019 (target=10000) : r =  +0.0000 0328 1717 2552 7718 9258 0230 0496 9818 5614 1889 3428 
count = 987 (target=10000) : r =  +0.0000 9441 3725 9466 5110 7853 0270 4252 2777 1022 3402 9711 
count = 28103 (target=10000) : r =  +0.0000 1678 0019 7630 1925 9726 3947 9228 1933 0986 1103 0495 
count = 9336 (target=10000) : r =  +0.0000 8762 7651 3963 6098 5643 5846 7482 1359 7096 0955 3438 
count = 37395 (target=10000) : r =  +0.0000 8892 5395 3175 7347 5619 1378 8813 1862 7116 9164 9683 
count = 5346 (target=10000) : r =  +0.0000 8691 4705 1787 4205 1174 7410 8097 1749 0386 4625 5414 
m_set_a()   3 --> x
x =  +3.
m_log()     ln(x) --> y
ln(3)= +1.0986 1228 8668 1096 9139 5245 2369 2252 5704 6474 9055 7822 
7494 5173 4694 3336 3749 4293 2186 0896 6873 6157 5481 3732 0887 8797 0029 
0659 5786 5742 3680 0422 5930 5198 2105 2801 8707 6727 7410 6031 6276 9183 
3813 6717 9373 6988 4436 0959 9037 4257 0316 7959 1152 1145 5919 1775 0671 
3470 5494 0166 7755 8022 2203 1702 5294 6897 5606 9010 6521 5056 4286 8138 
0363 1737 3298 5777 8236 6991 6547 9213 1818 1490 2003 0103 8236 3012 2248 
6527 4819 8225 9910 9745 2490 8964 5805 3467 0088 4596 5085 7484 4411 9018 
8570 8764 7494 8670 7961 3085 8294 1160 2166 1212 
ln2()       ln(2) --> x
ln(2)= +0.6931 4718 0559 9453 0941 7232 1214 5817 6568 0755 0013 4360 2552 
5412 0680 0094 9339 3621 9696 9471 5605 8633 2699 6418 6875 4200 1481 0205 
7068 5733 6855 2023 5758 1305 5703 2670 7516 3507 5961 9307 2757 0828 3714 
3519 0307 0386 2389 1673 4711 2335 0115 3644 9795 5239 1204 7517 2681 5749 
3206 5155 5247 3413 9525 8829 5045 3007 0953 2636 6642 6541 0423 9157 8149 
5204 3740 4303 8550 0801 9441 7064 1671 5186 4471 2839 9681 7178 4546 9570 
2627 1631 0645 4615 0257 2074 0248 1637 7733 8963 8550 6952 6066 8341 1372 
7387 3722 9289 5649 3547 0257 6265 2098 8596 9320 
ln625()     ln(0.625) --> y
ln(.625)= -0.4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 
3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 
3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 
3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 
3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 
3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 
3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 1538 4615 
3846 1538 4615 3846 1538 4615 3846 1538 4615 3846 
m_mul_s()   x*4 --> z
4*ln(2) =  +2.7725 8872 2239 7812 3766 8928 4858 3270 6272 3020 0053 7441 
0210 1648 2720 0379 7357 4487 8787 7886 2423 4533 0798 5674 7501 6800 5924 
0822 8274 2934 7420 8094 3032 5222 2813 0683 0065 4030 3847 7229 1028 3313 
4857 4076 1228 1544 9556 6693 8844 9340 0461 4579 9182 0956 4819 0069 0726 
2997 2826 0622 0989 3655 8103 5318 0181 2028 3813 0546 6570 6164 1695 6631 
2598 0817 4961 7215 4200 3207 7766 8256 6686 0745 7885 1359 8726 8713 8187 
8281 0508 6524 2581 8460 1028 8296 0992 6551 0935 5855 4202 7810 4267 3364 
5490 9549 4891 7158 2597 4188 1030 5060 8395 4388 
m_add1()    y+z --> y
ln(10)= +2.3110 5026 0701 3196 9920 7390 0242 9424 4733 8404 6207 5902 
5594 7802 1181 5764 3511 2949 4172 4040 0884 9917 6952 4136 2886 2954 4385 
6207 4428 1396 2805 4248 1494 0606 8966 9144 5450 0184 2309 2613 7182 1775 
0242 0229 9689 6929 5710 5155 4229 5493 8922 9964 5335 9418 0203 6222 9187 
8381 8979 9083 6373 9809 6565 0702 6335 0489 9197 6700 5032 1548 7849 5092 
7982 6971 3423 2600 0354 1669 3151 4410 5147 6130 4038 9821 4111 4867 6649 
3665 6662 4985 7966 4613 9490 3680 7146 5012 6320 2009 2664 3195 0421 1826 
0875 5703 3353 2542 8751 2649 6415 1214 6856 9772 
m_div()     x/y --> z
log10(2)= +0.2999 2735 0065 3336 5072 0924 9334 0384 9196 2942 6325 6509 5274 
7016 9720 3688 1936 8182 2929 9933 9794 4117 7688 6174 2882 1709 8213 5257 
5742 1864 2326 0078 2456 3854 2121 0531 5696 2405 2699 2328 9410 4407 8590 
8838 2146 4124 9969 1398 5722 7038 6319 4037 5045 9058 1565 1579 5400 7840 
1984 3570 4380 9508 5033 8972 9111 3014 7214 9605 7475 1826 9672 5250 7992 
2700 2300 1234 5126 3746 9848 7221 2453 8871 7421 6106 3193 0744 5741 3399 
3525 1062 2653 3099 3796 2305 2866 9132 6304 3633 7923 1285 2440 4451 1681 
3224 6763 0431 6454 6974 8028 7156 5985 9823 6758 
ln75()      ln(0.75) --> z
ln(.75)= -0.2876 8207 2451 7809 2743 9219 0059 9382 7431 5035 0971 0897 7610 
5650 6665 6853 4929 2950 7207 8046 4338 1108 9917 9105 2862 9603 2932 9751 
8350 5725 0030 3624 5585 7412 9301 2539 6325 0287 4513 2582 9237 2473 3615 
0320 1240 3784 0342 2387 0385 0412 9913 9330 8438 9332 6490 3259 4691 8028 
0919 0144 2738 8805 6848 5956 4795 9116 6299 6262 6763 8025 6561 0177 5935 
8671 4182 2829 8863 4612 2335 4915 1524 8882 6939 5576 1127 1344 6845 2613 
0434 5036 1379 9484 8023 5183 4690 9808 5379 3331 2015 6420 7721 7663 4174 
6009 9950 9908 3337 4008 2221 1370 2031 5982 0240 
m_mul1_s()  x*2 --> x
x =  +1.3862 9436 1119 8906 1883 4464 2429 1635 3136 1510 0026 8720 5105 
0824 1360 0189 8678 7243 9393 8943 1211 7266 5399 2837 3750 8400 2962 0411 
4137 1467 3710 4047 1516 2611 1406 5341 5032 7015 1923 8614 5514 1656 7428 
7038 0614 0772 4778 3346 9422 4670 0230 7289 9591 0478 2409 5034 5363 1498 
6413 0311 0494 6827 9051 7659 0090 6014 1906 5273 3285 3082 0847 8315 6299 
0408 7480 8607 7100 1603 8883 4128 3343 0372 8942 5679 9363 4356 9093 9140 
5254 3262 1290 9230 0514 4148 0496 3275 5467 7927 7101 3905 2133 6682 2745 
4774 7445 8579 1298 7094 0515 2530 4197 7194 
m_add1()    z+x --> z
ln(3)= +1.0986 1228 8668 1096 9139 5245 2369 2252 5704 6474 9055 7822 
7494 5173 4694 3336 3749 4293 2186 0896 6873 6157 5481 3732 0887 8797 0029 
0659 5786 5742 3680 0422 5930 5198 2105 2801 8707 6727 7410 6031 6276 9183 
3813 6717 9373 6988 4436 0959 9037 4257 0316 7959 1152 1145 5919 1775 0671 
3470 5494 0166 7755 8022 2203 1702 5294 6897 5606 9010 6521 5056 4286 8138 
0363 1737 3298 5777 8236 6991 6547 9213 1818 1490 2003 0103 8236 3012 2248 
6527 4819 8225 9910 9745 2490 8964 5805 3467 0088 4596 5085 7484 4411 9018 
8570 8764 7494 8670 7961 3085 8294 1160 2166 1212 
m_div()     z/y --> x
log10(3)= +0.4753 7360 2794 2208 4886 7418 8198 8450 0018 9836 9634 5581 8457 
6309 1427 7689 6206 9882 5677 5646 5772 0379 2546 2624 2377 9766 9832 6261 
3765 9374 3585 7648 6755 7432 8311 7124 0336 3364 0419 9213 3789 4337 4763 
0278 5866 9358 4906 6112 2528 5683 3410 5218 6398 5836 5471 9092 9927 5563 
3950 1028 7890 0338 0198 3874 2740 3750 9027 1809 4307 5399 5793 8054 1654 
1007 9172 4161 5901 6843 5958 8397 9924 0819 0408 7363 8695 8165 9106 9288 
4546 0534 6150 8596 5044 1762 9089 2719 2950 7115 5209 5566 8716 7981 0092 
1719 1841 8637 5605 8843 9736 3633 4525 3391 7773 
π --> x
x =  +3.1415 9265 3589 7932 3846 2643 3832 7950 2884 1971 6939 9375 1058 
2097 4944 5923 0781 6406 2862 0899 8628 0348 2534 2117 0679 8214 8086 5132 
8230 6647 0938 4460 9550 5822 3172 5359 4081 2848 1117 4502 8410 2701 9385 
2110 5559 6446 2294 8954 9303 8196 4428 8109 7566 5933 4461 2847 5648 2337 
8678 3165 2712 0190 9145 6485 6692 3460 3486 1045 4326 6482 1339 3607 2602 
4914 1273 7245 8700 6606 3155 8817 4881 5209 2096 2829 2540 9171 5364 3678 
9259 0360 0113 3053 0548 8204 6652 1384 1469 5194 1511 6094 3305 7270 3657 
5959 1953 0921 8611 7381 9326 1179 3105 1185 
m_div1_s()  x/3 --> x
x =  +1.0471 9755 1196 5977 4615 4214 4610 9316 7628 0657 2313 3125 0352 
7365 8314 8641 0260 5468 7620 6966 6209 3449 4178 0705 6893 2738 2695 5044 
2743 5549 0312 8153 6516 8607 4390 8453 1360 4282 7039 1500 9470 0900 6461 
7370 1853 2148 7431 6318 3101 2732 1476 2703 2522 1977 8153 7615 8549 4112 
6226 1055 0904 0063 6381 8828 5564 1153 4495 3681 8108 8827 3779 7869 0867 
4971 3757 9081 9566 8868 7718 6272 4960 5069 7365 4276 4180 3057 1788 1226 
3086 3453 3371 1017 6849 6068 2217 3794 7156 5064 7170 5364 7768 5756 7885 
8653 0651 0307 2870 5793 9775 3726 4368 3728 
m_tan()     tan(x) --> y
tan(x) = +1.7320 5080 7568 8772 9352 7446 3415 0587 2366 9428 0525 3810 
3806 2805 5806 9794 5193 3016 9088 0003 7081 1461 8675 7248 5756 7562 6141 
4154 0670 3029 9699 4509 4998 9524 7881 1655 5120 9437 3648 5280 9323 1902 
3055 8206 7974 8201 0108 4674 9232 6501 5312 3432 6690 3322 8866 5067 2254 
6689 2183 7971 2270 4713 1660 3678 6158 8019 0499 8653 7379 8593 8946 7650 
3475 0657 6050 7566 1834 8129 6061 0094 7602 1871 9032 5083 1458 2952 3959 
8329 9778 9824 5082 8871 4463 8329 1734 7224 1639 8458 7855 3976 6795 8063 
8183 5366 6110 8431 7378 0894 3783 1610 2088 3055 
sqrt(3)= +1.7320 5080 7568 8772 9352 7446 3415 0587 2366 9428 0525 3810 
3806 2805 5806 9794 5193 3016 9088 0003 7081 1461 8675 7248 5756 7562 6141 
4154 0670 3029 9699 4509 4998 9524 7881 1655 5120 9437 3648 5280 9323 1902 
3055 8206 7974 8201 0108 4674 9232 6501 5312 3432 6690 3322 8866 5067 2254 
6689 2183 7971 2270 4713 1660 3678 6158 8019 0499 8653 7379 8593 8946 7650 
3475 0657 6050 7566 1834 8129 6061 0094 7602 1871 9032 5083 1458 2952 3959 
8329 9778 9824 5082 8871 4463 8329 1734 7224 1639 8458 7855 3976 6795 8063 
8183 5366 6110 8431 7378 0894 3783 1610 2088 3055 
m_set_a()   0.5 --> x
x =  +0.5000 
m_acos()    arc cos(x) --> y
acos(x)= +1.0441 0455 6239 7321 5969 5881 2062 1568 2595 3269 7156 9589 
2437 2287 9444 3838 5373 7147 6046 9064 8524 9293 8265 2102 7667 7711 2247 
4344 7916 4012 4296 7714 0241 2066 1806 3012 7414 2268 1776 4365 5237 5285 
9282 6917 3716 2180 6796 9114 3210 7444 4376 0928 4740 9154 2934 0463 5825 
7205 4074 8951 9284 9274 6581 8029 4250 4611 7711 0577 7332 7691 0523 1209 
9922 0387 5883 8552 5829 4964 5613 8924 3503 8335 0780 0204 8315 9409 7144 
8465 2418 3213 0490 6652 4789 3933 6193 4677 1991 1333 3673 9667 9171 4032 
1801 5717 3469 4267 2178 9107 1924 2271 9894 3376 
sqrt(3) --> x
x =  +1.7320 5080 7568 8772 9352 7446 3415 0587 2366 9428 0525 3810 3806 
2805 5806 9794 5193 3016 9088 0003 7081 1461 8675 7248 5756 7562 6141 4154 
0670 3029 9699 4509 4998 9524 7881 1655 5120 9437 3648 5280 9323 1902 3055 
8206 7974 8201 0108 4674 9232 6501 5312 3432 6690 3322 8866 5067 2254 6689 
2183 7971 2270 4713 1660 3678 6158 8019 0499 8653 7379 8593 8946 7650 3475 
0657 6050 7566 1834 8129 6061 0094 7602 1871 9032 5083 1458 2952 3959 8329 
9778 9824 5082 8871 4463 8329 1734 7224 1639 8458 7855 3976 6795 8063 8183 
5366 6110 8431 7378 0894 3783 1610 2088 3055 
m_div_s()   x/2 --> v
v =  +0.8660 2540 3784 4386 4676 3723 1707 5293 6183 4714 0262 6905 1903 1402 
7903 4897 2596 6508 4544 0001 8540 5730 9337 8624 2878 3781 3070 7077 0335 
1514 9849 7254 7499 4762 3940 5827 7560 4718 6824 2640 4661 5951 1527 9103 
3987 4100 5054 2337 4616 3250 7656 1716 3345 1661 4433 2533 6127 3344 6091 
8985 6135 2356 5830 1839 3079 4009 5249 9326 8689 9296 9473 3825 1737 5328 
8025 3783 0917 4064 8030 5047 3801 0935 9516 2541 5729 1476 1979 9164 9889 
4912 2541 4435 7231 9164 5867 3612 0819 9229 3927 6988 3397 9031 9091 7683 
3055 4215 8689 0447 1891 5805 1044 1527 6246 
m_asin()    arc sin(v) --> y
asin(x)= +1.0441 0455 6239 7321 5969 5881 2062 1568 2595 3269 7156 9589 
2437 2287 9444 3838 5373 7147 6046 9064 8524 9293 8265 2102 7667 7711 2247 
4344 7916 4012 4296 7714 0241 2066 1806 3012 7414 2268 1776 4365 5237 5285 
9282 6917 3716 2180 6796 9114 3210 7444 4376 0928 4740 9154 2934 0463 5825 
7205 4074 8951 9284 9274 6581 8029 4250 4611 7711 0577 7332 7691 0523 1209 
9922 0387 5883 8552 5829 4964 5613 8924 3503 8335 0780 0204 8315 9409 7144 
8465 2418 3213 0490 6652 4789 3933 6193 4677 1991 1333 3673 9667 9171 4032 
1801 5717 3469 4267 2178 9107 1924 2271 9894 3376 
m_atan()    arc tan(x) --> y
atan(x)= +1.0441 0455 6239 7321 5969 5881 2062 1568 2595 3269 7156 9589 
2437 2287 9444 3838 5373 7147 6046 9064 8524 9293 8265 2102 7667 7711 2247 
4344 7916 4012 4296 7714 0241 2066 1806 3012 7414 2268 1776 4365 5237 5285 
9282 6917 3716 2180 6796 9114 3210 7444 4376 0928 4740 9154 2934 0463 5825 
7205 4074 8951 9284 9274 6581 8029 4250 4611 7711 0577 7332 7691 0523 1209 
9922 0387 5883 8552 5829 4964 5613 8924 3503 8335 0780 0204 8315 9409 7144 
8465 2418 3213 0490 6652 4789 3933 6193 4677 1991 1333 3673 9667 9171 4032 
1801 5717 3469 4267 2178 9107 1924 2271 9894 3376 
π --> z
z =  +3.1415 9265 3589 7932 3846 2643 3832 7950 2884 1971 6939 9375 1058 
2097 4944 5923 0781 6406 2862 0899 8628 0348 2534 2117 0679 8214 8086 5132 
8230 6647 0938 4460 9550 5822 3172 5359 4081 2848 1117 4502 8410 2701 9385 
2110 5559 6446 2294 8954 9303 8196 4428 8109 7566 5933 4461 2847 5648 2337 
8678 3165 2712 0190 9145 6485 6692 3460 3486 1045 4326 6482 1339 3607 2602 
4914 1273 7245 8700 6606 3155 8817 4881 5209 2096 2829 2540 9171 5364 3678 
9259 0360 0113 3053 0548 8204 6652 1384 1469 5194 1511 6094 3305 7270 3657 
5959 1953 0921 8611 7381 9326 1179 3105 1185 
m_div1_s()  z/3 --> z
pi / 3 = +1.0471 9755 1196 5977 4615 4214 4610 9316 7628 0657 2313 3125 
0352 7365 8314 8641 0260 5468 7620 6966 6209 3449 4178 0705 6893 2738 2695 
5044 2743 5549 0312 8153 6516 8607 4390 8453 1360 4282 7039 1500 9470 0900 
6461 7370 1853 2148 7431 6318 3101 2732 1476 2703 2522 1977 8153 7615 8549 
4112 6226 1055 0904 0063 6381 8828 5564 1153 4495 3681 8108 8827 3779 7869 
0867 4971 3757 9081 9566 8868 7718 6272 4960 5069 7365 4276 4180 3057 1788 
1226 3086 3453 3371 1017 6849 6068 2217 3794 7156 5064 7170 5364 7768 5756 
7885 8653 0651 0307 2870 5793 9775 3726 4368 3728