sistema di equazioni non lineare a tre variabili di stato

[raffa]

Ciao a tutti vi dico il mio problema. ho il seguente sistema di
equazioni differenziali non lineari del terzo ordine:

dx1/dt= x1* (r*(1-(x1/50))-(2/(8+x1)*x2));
dx2/dt=x2*(-1.5+(2*(2/8+x1)*x1)-((2/5+x1)*x3)
dx3/dt=x3*(-0.5+(2*(2/5+x2)*x2)

Dato questo sistema devo trovare in matlab i massimi locali della
prima variabile di stato. una volta trovati li dovei plottare in un
grafico in cui in ordinata metto questi massimi locali in ascissa il
parametro di controllo r che dovrà variare tra 0:3 con un passo pari a
0.1
Come posso risolvere questo?