x(0)=1
x(1600)=1/2

C=1, k=ln[1/2] /(-1600)= (ln1 - ln2) / (-1600) = ln2 / 1600
x(100)=x(0)+int(0, 100) x'(t)dt= x(0)+int (0, 100) (-k*C*e^(-k*t))dt =
=1-ln2 / 1600 * int (0, 100) (e^(-t*ln2 / 1600))dt ~ 0.9576 , т.е. останется ~95.76% радия и распадется, соответственно, ~4.24%