Skip to content

Conversation

ChrisRackauckas-Claude
Copy link

Summary

  • Fixed incorrect partial derivative calculation in the Van der Pol Jacobian
  • The derivative ∂f₂/∂x was missing the y term: changed from μ * (-2x - 1) to μ * (-2xy - 1)

Details

The Van der Pol equation is:

f₁ = y
f₂ = μ * ((1 - x²) * y - x)

The Jacobian element J[2,1] = ∂f₂/∂x should be:

∂f₂/∂x = μ * (∂/∂x[(1 - x²) * y - x])
       = μ * (-2xy - 1)

The previous implementation incorrectly calculated this as μ * (-2x - 1), missing the y term.

Test Plan

🤖 Generated with Claude Code

The partial derivative ∂f₂/∂x was incorrect. For the Van der Pol equation:
  f₂ = μ * ((1 - x²) * y - x)

The correct partial derivative is:
  ∂f₂/∂x = μ * (-2xy - 1)

The previous implementation was missing the y term in the derivative.

This fixes the failing test in OrdinaryDiffEq.jl CI.

🤖 Generated with [Claude Code](https://claude.ai/code)

Co-Authored-By: Claude <[email protected]>
@ChrisRackauckas ChrisRackauckas merged commit 50a3c6a into SciML:master Aug 9, 2025
14 of 19 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants