velocityold

velocitynew

time

Eulers method for solving dv/dt=rhs(v)

plot

burnvelocity = 1.4 (slider from 0 to 5, step 0.1)
burnrate = 0.2 (slider from 0 to 1, step 0.1)
drag = 0.1 (slider from 0 to 1, step 0.1)
Force = 1.1 (slider from 0 to 5, step 0.1)
mass = 0.7 (slider from 0 to 1, step 0.1)
friction = 0.5 (slider from 0 to 1, step 0.1)
deltat = 0.49 (slider from 0 to 1, step 0.01)
startv = 0 (slider from -5 to 5, step 1)
B3: startv
rhs(v) = Force / mass - ((drag * v) * abs(v)) - (friction * (1 - ℯ^((-v)))) + (burnvelocity * burnrate)
C3: B3 + (deltat * rhs(B3))
B4: C3
A3: 0
A4: A3 + deltat
D3 = Segment((A3, B3), (A3 + deltat, C3))
C4: B4 + (deltat * rhs(B4))
D4 = Segment((A4, B4), (A4 + deltat, C4))
B5: C4
A5: A4 + deltat
C5: B5 + (deltat * rhs(B5))
D5 = Segment((A5, B5), (A5 + deltat, C5))
B6: C5
A6: A5 + deltat
C6: B6 + (deltat * rhs(B6))
D6 = Segment((A6, B6), (A6 + deltat, C6))
B7: C6
A7: A6 + deltat
C7: B7 + (deltat * rhs(B7))
D7 = Segment((A7, B7), (A7 + deltat, C7))
B8: C7
A8: A7 + deltat
C8: B8 + (deltat * rhs(B8))
D8 = Segment((A8, B8), (A8 + deltat, C8))
B9: C8
A9: A8 + deltat
C9: B9 + (deltat * rhs(B9))
D9 = Segment((A9, B9), (A9 + deltat, C9))
B10: C9
A10: A9 + deltat
C10: B10 + (deltat * rhs(B10))
D10 = Segment((A10, B10), (A10 + deltat, C10))
B11: C10
A11: A10 + deltat
C11: B11 + (deltat * rhs(B11))
D11 = Segment((A11, B11), (A11 + deltat, C11))
B12: C11
A12: A11 + deltat
C12: B12 + (deltat * rhs(B12))
D12 = Segment((A12, B12), (A12 + deltat, C12))
B13: C12
A13: A12 + deltat
C13: B13 + (deltat * rhs(B13))
D13 = Segment((A13, B13), (A13 + deltat, C13))
B14: C13
A14: A13 + deltat
C14: B14 + (deltat * rhs(B14))
D14 = Segment((A14, B14), (A14 + deltat, C14))
B15: C14
A15: A14 + deltat
C15: B15 + (deltat * rhs(B15))
D15 = Segment((A15, B15), (A15 + deltat, C15))