Skip to content

Instantly share code, notes, and snippets.

using DiffEqOperators
using OrdinaryDiffEq
A = -10*speye(3)
M = speye(3)
num_steps = 1
function update_A(A,u,p,t)
using OrdinaryDiffEq, DiffEqOperators,Sundials
n = 600
M = rand(n,n)
K = rand(n,n)
ϵ = 1e-2
p = (ϵ,M,K)
u0 = ones(n)