scipy-numerical

Solve physics initial value problems and discretize PDEs with SciPy numerical methods.

33|6|Updated Mar 17, 2026
One-click install
npx skills add https://github.com/xjtulyc/awesome-rosetta-skills --skill scipy-numerical
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: scipy-numerical
Source: https://github.com/xjtulyc/awesome-rosetta-skills/tree/main/skills/01-physics/scipy-numerical
Command: npx skills add https://github.com/xjtulyc/awesome-rosetta-skills --skill scipy-numerical

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires scipy, numpy, matplotlib, and includes assets (resource) components.

What problem does it solve?

It helps you compute and analyze physics systems numerically when closed-form solutions are unavailable, covering time evolution, spectral content, and parameter optimization.

Core Features & Use Cases

  • ODE/IVP solving: Simulate nonlinear dynamics with scipy.integrate.solve_ivp, including method selection for non-stiff vs stiff regimes (e.g., RK45, Radau, BDF).
  • PDE discretization: Solve a 1D diffusion/heat equation using sparse finite-difference matrices and time-stepping schemes like Crank–Nicolson.
  • Spectral analysis with FFT: Perform FFT-based power spectral density estimation and compare methods such as Welch PSD for noisy signals.
  • Optimization & fitting: Fit physically motivated models (e.g., damped cosines) with uncertainty estimation via covariance matrices.
  • Sparse linear algebra: Compute eigenvalues/eigenstates for quantum-like operators using sparse eigensolvers (e.g., eigsh).

Quick Start

Use the scipy-numerical skill to simulate a Lorenz attractor with solve_ivp and visualize the 3D trajectory for method RK45.

Frequently Asked Questions about scipy-numerical

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I solve stiff ODE initial value problems using scipy.integrate?

Solve stiff ODE initial value problems by configuring scipy.integrate.solve_ivp with stiff integration methods like Radau or BDF, selecting appropriate numerical tolerances for stable nonlinear dynamics time evolution.

What is the best way to solve a 1D diffusion PDE with sparse matrices?

The best way to solve a 1D diffusion PDE is applying sparse finite-difference matrices for spatial discretization and time-stepping schemes like Crank–Nicolson to compute heat equation evolution robustly.

How do I analyze noisy signals using FFT and Welch PSD in Python?

Analyze noisy signals by computing FFT-based power spectral density estimation and comparing methods like Welch PSD using SciPy to extract frequency content from noisy measurement data.

Can I fit damped cosine models to data and estimate parameter uncertainty?

Fit physically motivated models like damped cosines using nonlinear parameter estimation with SciPy, calculating parameter uncertainty via covariance matrices extracted during the curve fitting process.

Does scipy support eigenvalue computation for quantum-inspired sparse operators?

SciPy supports eigenvalue computation for quantum-like operators using sparse eigensolvers such as eigsh, enabling efficient eigenstate calculation for sparse matrix assembly and factorization problems.