In your case you may want to try with arbitrary precision floats (BigFloat)
So just define x = BigFloat[1, 0.045126, 6.9315] w1 = w2 = w3 = BigFloat[1 1 1] w4 = BigFloat[1 -10 -10] and call your function as compute(f1[1], f2[1], f3[1], d[1], BigFloat(0), BigFloat(153))
