Dag Freek,
Algemeen:
"het is in een ruwe fase die nu niet werkt"
Bedoel je dan dat je dit in één keer in elkaar hebt gestoken zonder voorafgaande vereenvoudigde versies (bijv nog zonder die luchtweerstand) te bouwen en op correcte werking te testen? Want dat is namelijk een absolute aanrader. In een volledig model vallen kleine foutjes in de wirwar van verbanden nauwelijks meer op en wordt debuggen gauw een enorme puzzel.
Dus mijn advies is om eerst even opnieuw te beginnen zonder die luchtweerstand, en controleren dat dat een mooie parabool oplevert. Of nog een stap verder terug, gewoon verticaal afvuren en zorgen dat dat voorspelde resultaten oplevert.
Verder, om te debuggen helpt het ook als je v-t en s-t grafieken hebt. Daarin vallen wellicht gekke fouten te herkennen. Anderen hier hebben deze software overigens wel en kunnen die wellicht draaien
natuurkundig:
Eén fout (waarschijnlijk, dwz als ik dit systeem van realties goed begrijp) valt me al wel op. Die zal niet veroorzaken dat je model NIET werkt, maar wel dat het natuurkundig niet correct werkt.
Je gebruikt als ik het goed begrijp vy om een Fw,y te berekenen, en vx voor een Fw,x. Maar vanwege het kwadratische verband gaat dat mis. De bal heeft alleen maar last van een diagonale weerstandskracht als gevolg van zijn diagonale snelheid. Dié diagonale Fw moet je dus eerst berekenen, op basis van (onder andere) V², en dán pas kun je gaan ontbinden.
Groet, Jan