Ik heb nog even verder gesimuleerd met LTspice. nu met de differentiaalvergelijking zelf gebruikt.
blok X7 beschrijft de dome: verband tussen r en versnelling a1. a1 wordt dan gecopieerd naar a en dan vervolgens geintegreerd om snelheid v te krijgen en nogmaals geintegreerd on afstand r te krijgen en dan is het loopje rond. als je dat gaat simuleren blijft alles op 0 staan dus dat is keurig de ene oplossing van de differentiaalvergelijking r=0. in het modelletje zit echter ook nog de optie om de dome te verplaatsen tov het balletje via een x coordinaat rdome. die staat normaal op 0 dus als r=0 dan zit het balletje keurig precies midden op de dome in rust.
maar nu kan ik rdome een kleine tijdelijke verplaatsing geven (een stapje van 1 nm positief en negatief gedurende 1ms op t=10ms) en dan zie je dat het systeem direct begint met het uitvoeren van de andere oplossing waarbij r steeds sneller toeneemt.
dus het klopt allemaal precies met de theorie en je ziet dus dat er ergens een verstoring nodig is om de bal te laten rollen.