/* Polynom x^4 + x^3 + x^2 + x [program = x + x * x + x * (x * x + (x + 0)) * x,res = 8.842604099730679,count = 1] [program = (x + x / x) * (x + x * x * x),res = 8.842604099730679,count = 140] [program = (x + x * x * x) * (x + x / x),res = 8.842604099730679,count = 118] [program = (x + x * x * x) * (1 + x),res = 8.842604099730679,count = 76] [program = (1 + x) * (x + x * x * x),res = 8.842604099730679,count = 70] [program = (x + x * x * x) * (x - -1),res = 8.842604099730679,count = 38] [program = (x - -1) * (x + x * x * x),res = 8.842604099730679,count = 21] resultMap = [8.842604099730679 = 6] Using pow2,pow3,pow4 instead of pow_mod: [program = pow2(x * (x - -4 + -3)) - (pow3(x) - x),res = 8.8426,count = 1] [program = pow3(x) + pow4(x) + (x + x * x),res = 8.8426,count = 92] [program = x + x * x + (pow3(x) + pow4(x)),res = 8.8426,count = 90] [program = (1 + x) * (x + pow3(x)),res = 8.8426,count = 109] [program = (x + pow3(x)) * (1 + x),res = 8.8426,count = 84] [program = (x + 1) * (x + pow3(x)),res = 8.8426,count = 73] [program = (x + pow3(x)) * (x + 1),res = 8.8426,count = 54] Cf p4_2.conf */ import util. data(p4_2,Data,Vars,Unknown,Ops,Constants,MaxSize,Params) :- Data = [ [[X],Y] : [X,Y] in chunks_of([ 1.34409048,8.84260401, 3.90205296,310.37239860, -3.15213438,74.18748889, 3.35929573,179.90161287, -3.19920446,79.04550313, 4.96816180,761.51040708, 4.30692881,446.83751274, -4.24947609,263.16439510, 4.57876301,561.07211193, 0.97705025,3.77570926, -1.07333129,0.16938441, 3.19436503,150.11449372, 0.95060075,3.52981359, 4.49502293,523.77459626, 1.70531921,18.02981955, 0.08083968,0.08794574, 3.95182523,325.17236887, 1.52566301,12.82245996, 0.34589428,0.52123531, 1.60021173,14.81558533, 4.84594306,693.58618055, -2.35157391,20.75414994, 0.86667165,2.83294645, 4.40422388,485.48187407, 2.70009240,82.82698775, -2.38631390,22.14659192, -1.45913850,2.09632019, 4.64952478,594.12282109, -1.57368729,3.13859193, 1.91581155,26.08916435 ],2)], Unknown = [1.34409048], Vars = ['x'], % Ops = [+,-,*,/,log,pow_mod2], Ops = [+,-,*,/,log,pow2,pow3,pow4,pow5], Constants = 1..10, MaxSize = 51, Params = new_map([ approx=1, init_size=1000, show_best=1, num_gens=180 % show_only_improvements=true % stop_criteria=generations, % remove_dups=true ]).