1. p> p>
Cambiar a G=1/3*((5*Pp-2*P1)/(P1 2*Pp)-P1*c1^2/(Pp*cp^2));
El error en la expresión provoca una oscilación de alta frecuencia y el tamaño del paso de integración es muy pequeño (aproximadamente 10 ^ -8), por lo que no se puede calcular durante mucho tiempo.
2. A juzgar por los resultados, el tiempo de simulación es de 50 segundos (ha entrado en estado estable).
3. Cuando ocurre ocupado, puede considerar permitir que la función ode* no devuelva parámetros. En este caso, la curva de respuesta se dibujará automáticamente y la operación se puede detener en cualquier momento a través de la GUI. Facilitar la observación del proceso de solución.