Algoritmia - Isoclinas

import numpy as np

import matplotlib.pyplot as plt

# Define la ecuación diferencial (ejemplo: dy/dx = x^2 - y^2)

def differential_equation(x, y):

    return x**2 - y**2

# Define el rango de valores para x y y

x_values = np.linspace(-2, 2, 100)

y_values = np.linspace(-2, 2, 100)

# Crea una cuadrícula para x y y

X, Y = np.meshgrid(x_values, y_values)

# Calcula las derivadas para cada punto en la cuadrícula

U = 1

V = differential_equation(X, Y)

# Normaliza las flechas para que tengan la misma longitud

N = np.sqrt(U**2 + V**2)

U /= N

V /= N

# Grafica las isoclinas usando flechas

plt.quiver(X, Y, U, V, scale=20, color='blue')

# Configura el aspecto del gráfico

plt.xlabel('x')

plt.ylabel('y')

plt.title('Isoclinas de la ecuación diferencial')

plt.grid(True)

plt.show()