def kestoNopeudenFunktiona(vMin, vMax, dt, dv): v = [] t = [] while vMin < vMax: v.append(vMin) t.append(matka(vMin, dt, 48*3600)) vMin += dv pylab.plot(v,t) pylab.show()
def haarukointi(vMin, vMax, vEpsilon, dt, max_t): vKeski = (vMin + vMax) / 2 while vMax - vMin >= vEpsilon: if matka(vKeski, dt, max_t) < 0: # Eli: jos kuuhun ei paasta... vMin = vKeski # ...nopeuden on oltava enemman kuin vKeski else: # Muuten... vMax = vKeski # ...etsitty nopeus on pienempi kuin vKeski vKeski = (vMin + vMax) /2 #print vKeski return vKeski
def kestoAikavalinFunktiona(vAlku, dtMin, dtMax, ddt): dt = [] t = [] while dtMin < dtMax: dt.append(dtMin) t.append(matka(vAlku, dtMin, 48*3600)) dtMin += ddt print dtMin pylab.plot(dt,t) pylab.show()