oscilující dipól

Needs["Graphics`PlotField`"];
Needs["Graphics`Graphics`"];

\(
\)

k = 2 Pi;

    ω = 2 Pi;

    r[x_, y_] := Sqrt[x^2 + y^2];

    ϑ[x_, y_] := ArcTan[x/y];

    contourDiOsc = 
        Table[pomDiOsc, {pomDiOsc, -0.9, 0.9, 
            0.2}];


    DiOscMovie = 
        Table[
ContourPlot[
((Cos[
                        k r[x, y] - ω t])/(k r[x, y]) - 
                  Sin[k r[x, y] - ω t]) (Sin[ϑ[x, 
                      y]])^2, 
{x, -1, 
              1}, 
{y, -1, 
              1}, 
PlotPoints -> 
              110, 
Contours -> 
              contourDiOsc, 
ContourShading -> 
              False, 
ContourLines -> 
              True, 
FrameTicks -> 
              None
], 
{t, 0, 1 - 1/20, 
            1/20}
];

Export["dipole-oscilujici.avi", DiOscMovie,  "AVI", 
        ImageSize -> {400, 300}];

Export["dipole-oscilujici.gif", DiOscMovie,  "GIF", 
        ConversionOptions -> {"Loop" -> True}];

DiOscPicture = 
        ContourPlot[
((Cos[
                      k r[x, y] - ω 0.6])/(k r[x, y]) - 
                Sin[k r[x, y] - ω 0.6]) (Sin[ϑ[x, 
                    y]])^2, 
{x, -1, 
            1}, 
{y, -1, 
            1}, 
PlotPoints -> 
            110, 
Contours -> 
            contourDiOsc, 
ContourShading -> 
            False, 
ContourLines -> 
            True, 
FrameTicks -> 
            None
];

Export["dipole-oscilujici.jpg", DiOscPicture, "JPEG", 
        ImageSize -> {640, 480}];