• Bat Andrea
    Bat Andrea
    2018-04-08

    BesselZeros = N@Table[BesselJZero[m, n], {m, 0, 5}, {n, 1, 5}];

    λ[m_, n_] := BesselZeros[[m + 1, n]];

    u[m_, n_, r_, θ, t, c_, A_, B_, C_, D_]
    := (A Cos[c λ[m, n] t] + B Sin[c λ[m, n] t])
    BesselJ[m, λ[m, n] r] (C Cos[m θ] + D Sin[m θ]);

    diskChladni = With[
    {m = 1,
    n = 2,
    c = 1,
    cols = RGBColor /@ {"#D7F2F7", "#354D62"}},
    ParallelTable[
    Show[
    Graphics[{Thickness[.005], cols[[1]], Circle[]}],
    ContourPlot[
    Cos[t] u[2, 3, Norm[{-y, x}], ArcTan[-y, x], 0, c, 1/2, 0, 1/2, 0]
    + Sin[t] u[3, 2, Norm[{-y, x}], ArcTan[-y, x], 0, c, 1/2, 0, 1/2, 0] == 0,
    {x, -1, 1}, {y, -1, 1}, PlotPoints -> 100,
    ContourStyle -> Directive[Thickness[.005], cols[[1]]],
    RegionFunction -> Function[{x, y, z}, x^2 + y^2 < .99]],
    Background -> cols[[-1]], PlotRange -> Sqrt[2],
    ImageSize -> 540],
    {t, 0., π - #, #}] &[π/200]
    ];

    Export[NotebookDirectory[] <> "diskChladni.gif", diskChladni, "DisplayDurations" -> {3/100}]

    0