Numerisk lösning Diffekvation

Hej! 

Jag skulle behöva hjälp med hur jag kan ta hjälp av exempelvis Geogebra för att lösa en andra ordningens icke linjär diffekvation numeriskt. Jag har förstått att man bör göra någon slags variabelsubstitution av andraderivatan för att sedan få ett system av två första ordningens diffar. Min fråga är dock hur jag kan ta hjälp av Geogebra eller liknande för att plotta detta för att se hur olika begynnelsevillkor påverkar lösningen.

Frågan handlar om en matematisk pendel:

$y\text{'}\text{'}+\frac{g}{l}\sin \left(y\right)=0$