Bitcoin Forum

package main

import (
	"os"
	"fmt"
	"strconv"
	"math/rand"
	"math"
)



var block = 0//block count ( for genesis block = 0 )

var maxBlock = []int { 
	262980,
	525960,
	788940,
	1051920,
	1314900,
	1577980,
	1840850,
	2103840,
	2366820,
	2629800,
	2892780,
	4733640,
	4996620,
	5259600,
}


var percent_year = 0.0314159
var percent_day_multi = math.Pow(1+percent_year, float64(1.0/365.0))
var percent = (percent_day_multi - 1) / 720.0

var PoS_percent = float64(0.01);

var Coin = int64(1000000)
var Cent = int64(10000)

var mplus, wMplus, sMplus float64

var wDiff = float64(1.0)
var sDiff = 0.0001

var supply = 0.0//moneysupply
var fees = 0.0

var wHour = 29//PoW blocks' rate per 1 hour
var shour = 1// PoS blocks' rate per 1 hour

// 1st year - hyper-maining.
var nblock_like_premine = 1000000
var bonus_premine = float64(20);


var cepochka_ostatok = float64(0);
var freq_pos = float64(0.2);
var max_pos2pow_dobavka = float64(0.25);

var ncount_pow = 0;
var ncount_pos = 0;

func main() {
	rnd := rand.New(rand.NewSource(755875))
	for {
		x := 0.0;
		if block < nblock_like_premine {x=1.0 - float64(block)/float64(nblock_like_premine)}
		wReward := bonus_premine * x + supply * percent;
		sReward := float64(0);
		if rnd.Float64() < freq_pos{ // PoS-block
			age_block := rnd.Int63n(60*720) + 30*720
			coin := (0.3 + rnd.Float64()*0.4)*supply/(60*720*freq_pos) // 60 - srednij vozrast v dnyah, 720*freq_pos - blockov v den'. 0.3+rnd*0.4, v pos zadeistvovany 30%-70% vseh sredstv
			coin_age := coin * float64(age_block) / float64(maxBlock[0]) //pseudorandom  coin_age in days
			sReward = float64(coin_age) * PoS_percent
			cepochka_ostatok += wReward - sReward // nedoemmissia, ili izlishki
//			fmt.Println("Age:", age_block, "Coin:", coin, "PoS: ", sReward)
			ncount_pos++
		}else {// PoW-block
			if cepochka_ostatok > 0{ // razmazywaem nevvedenyu emissiu iz pos po sledujshim pow
				if cepochka_ostatok > wReward * max_pos2pow_dobavka{
					cepochka_ostatok -= wReward * max_pos2pow_dobavka
					wReward += wReward * max_pos2pow_dobavka
				}else{
					wReward += cepochka_ostatok;
					cepochka_ostatok = 0
				}
			}
			ncount_pow++
		}

		if sReward == 0{
			supply += wReward
		}else{
			supply += sReward
		}

//		if block%(maxBlock[0]/10) == 0{
//		fmt.Println("moneysupply: ", strconv.FormatFloat(supply, 'f' , -1, 64), " at block #", block, "year #", block / maxBlock[0])
//		fmt.Println("PoS_blocks =", ncount_pos)
//		fmt.Println("PoW_blocks =", ncount_pow)
//		fmt.Println("==========================================================")
//		}
		block++
		if block > maxBlock[10] { break }
	}
		fmt.Println("moneysupply: ", strconv.FormatFloat(supply, 'f' , -1, 64), " at block #", block, "year #", block / maxBlock[0])
		fmt.Println("PoS_blocks =", ncount_pos," PoW_blocks =", ncount_pow)
		fmt.Println("==========================================================")

	os.Exit(0)
}

x	ln(Г(х))	x*ln(x)
1E+00	0	0
1E+01	12,80182748	23,02585093
1E+02	359,1342054	460,5170186
1E+03	5905,220423	6907,755279
1E+04	82099,7175	92103,40372
1E+05	1051287,709	1151292,546
1E+06	12815504,57	13815510,56
1E+07	151180949,4	161180956,5
1E+08	1742068066	1842068074
1E+09	19723265828	20723265837
1E+10	2,20259E+11	2,30259E+11
1E+11	2,43284E+12	2,53284E+12
1E+12	2,6631E+13	2,7631E+13
1E+13	2,89336E+14	2,99336E+14
1E+14	3,12362E+15	3,22362E+15
1E+15	3,35388E+16	3,45388E+16
1E+16	3,58414E+17	3,68414E+17
1E+17	3,81439E+18	3,91439E+18

x := wDiff + 1000.0
w1 := math.Log(x) * x + 1.0
y := float64(block)
 w2 := math.Log(y) * y + 1.0
ycur := float64(block) / float64(year);
guile := 20.0
ygui := 3.0;
w3 := guile * ycur / ygui;
if (ycur > ygui) {w3 = guile * (2.0 - ycur/ygui)}
if (ycur > ygui*2) {w3 = 0}
 wReward = math.Abs((w1 * w2)) / 100000000000.0 + w3;

// The message start string is designed to be unlikely to occur in normal data.
// The characters are rarely used upper ASCII, not valid as UTF-8, and produce
// a large 4-byte int at any alignment.
unsigned char pchMessageStart[4] = { 0xf9, 0xbe, 0xb4, 0xd9 };