(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     23780,        507]
NotebookOptionsPosition[     23085,        489]
NotebookOutlinePosition[     23533,        506]
CellTagsIndexPosition[     23490,        503]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "V\[YAcute]po\[CHacek]et", " ", "EGM", " ", "dle", " ", "Eriksson"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Wolfram", " ", 
    RowBox[{"Mathematica", " ", ".", "nb"}], " ", "file"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Jakub", " ", "Benda"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"05", "/", "2025"}], "*)"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"ClearAll", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Definice", " ", "z\[AAcute]kladn\[IAcute]ch", " ", "parametr\[URing]", 
     " ", "simulace"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"cc", "[", 
      RowBox[{"i_", "?", "NumericQ"}], "]"}], ":=", 
     RowBox[{"PDF", "[", 
      RowBox[{
       RowBox[{"LogNormalDistribution", "[", 
        RowBox[{"2.43822", ",", "0.655125"}], "]"}], ",", "i"}], "]"}]}], 
    ";"}], 
   RowBox[{"(*", 
    RowBox[{
    "Distribu\[CHacek]n\[IAcute]", " ", "funkce", " ", 
     "pravd\[EHacek]podobnosti", " ", "velikosti", " ", "proudu"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ccmin", "=", "2"}], ";"}], " ", 
   RowBox[{"(*", 
    RowBox[{
    "Spodn\[IAcute]", " ", "hranice", " ", "velikosti", " ", "proudu"}], 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ccmax", "=", "80"}], ";"}], " ", 
   RowBox[{"(*", 
    RowBox[{
    "Horn\[IAcute]", " ", "hranice", " ", "velikosti", " ", "proudu"}], 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Ng", "=", "2"}], ";", 
    RowBox[{"(*", 
     RowBox[{"Hustota", " ", "blesk\[URing]", " ", "\[CapitalCHacek]R", " ", 
      RowBox[{"(", 
       RowBox[{"blesky", "/", "km2"}], ")"}]}], "*)"}], "\[IndentingNewLine]", 
    RowBox[{"p1", "=", 
     RowBox[{"{", 
      RowBox[{"0.", ",", "30.32"}], "}"}]}], ";", " ", 
    RowBox[{"(*", 
     RowBox[{"Pozice", " ", "zemn\[IAcute]c\[IAcute]ho", " ", "lana"}], 
     "*)"}], "\[IndentingNewLine]", 
    RowBox[{"p2", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"12", "-", "6"}], ",", 
       RowBox[{"23.28", "-", "4.33"}]}], "}"}]}], ";"}], " ", 
   RowBox[{"(*", 
    RowBox[{
    "pozice", " ", "f\[AAcute]zov\[EAcute]ho", " ", "vodi\[CHacek]e"}], 
    "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "N\[AAcute]sleduje", " ", "definice", " ", "funkce", " ", 
     "\[SHacek]\[IAcute]\[RHacek]ky", " ", "selh\[AAcute]n\[IAcute]", " ", 
     "ochrany", " ", "zna\[CHacek]eno", " ", "Dc1"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"Dc1", "[", 
      RowBox[{
       RowBox[{"cc_", "?", "NumericQ"}], ",", 
       RowBox[{"p1y_", "?", "NumericQ"}], ",", 
       RowBox[{"p2x_", "?", "NumericQ"}]}], "]"}], ":=", 
     RowBox[{"Module", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
        "rad1", ",", "rad2", ",", "prusecik", ",", "prusecikCalc", ",", 
         "prusecikcalc2", ",", "prusecik2", ",", "x1", ",", "x2", ",", 
         "prumy", ",", "rg", ",", "w"}], "}"}], ",", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"(*", 
        RowBox[{
        "Z\[AAcute]kladn\[IAcute]", " ", "rovnice", " ", "pro", " ", 
         "polom\[EHacek]ry", " ", "kru\[ZHacek]nic", " ", 
         "doskokov\[YAcute]ch", " ", "vzd\[AAcute]lenost\[IAcute]"}], "*)"}], 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"rad1", "=", 
         RowBox[{"0.67", "*", 
          SuperscriptBox[
           RowBox[{"p1", "[", 
            RowBox[{"[", "2", "]"}], "]"}], "0.45"], "*", 
          SuperscriptBox["cc", "0.74"]}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"rad2", "=", 
         RowBox[{"0.67", "*", 
          SuperscriptBox[
           RowBox[{"p2", "[", 
            RowBox[{"[", "2", "]"}], "]"}], "0.45"], "*", 
          SuperscriptBox["cc", "0.74"]}]}], ";", "\[IndentingNewLine]", 
        "\[IndentingNewLine]", 
        RowBox[{"(*", 
         RowBox[{
         "V\[YAcute]po\[CHacek]et", " ", "pr\[URing]se\[CHacek]\[IAcute]ku", 
          " ", "doskokov\[YAcute]ch", " ", "vzd\[AAcute]lenost\[IAcute]", " ",
           "obou", " ", "vodi\[CHacek]\[URing]"}], "*)"}], 
        "\[IndentingNewLine]", 
        RowBox[{"prusecikCalc", "=", 
         RowBox[{"NSolve", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{"x", "-", 
                 RowBox[{"p1", "[", 
                  RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "2"], "+", 
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{"y", "-", "p1y"}], ")"}], "2"], "-", 
              SuperscriptBox["rad1", "2"]}], "==", "0"}], "&&", 
            RowBox[{
             RowBox[{
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{"x", "-", "p2x"}], ")"}], "2"], "+", 
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{"y", "-", 
                 RowBox[{"p2", "[", 
                  RowBox[{"[", "2", "]"}], "]"}]}], ")"}], "2"], "-", 
              SuperscriptBox["rad2", "2"]}], "==", "0"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"x", ",", "y"}], "}"}], ",", "Reals"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"prusecik", "=", 
         RowBox[{
          RowBox[{"Select", "[", 
           RowBox[{"prusecikCalc", ",", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"x", ">", 
                 RowBox[{"p2", "[", 
                  RowBox[{"[", "1", "]"}], "]"}]}], "&&", 
                RowBox[{"y", ">", 
                 RowBox[{"p1", "[", 
                  RowBox[{"[", "2", "]"}], "]"}]}]}], ")"}], "/.", "#"}], 
             "&"}]}], "]"}], "//", "Flatten"}]}], ";", "\[IndentingNewLine]", 
        "\[IndentingNewLine]", 
        RowBox[{"(*", 
         RowBox[{
         "Stanoven\[IAcute]", " ", "hodnot", " ", "krajn\[IAcute]ch", " ", 
          "bod\[URing]", " ", "\[UAcute]se\[CHacek]ky", " ", "Dc1"}], "*)"}], 
        "\[IndentingNewLine]", 
        RowBox[{"x1", "=", 
         RowBox[{"x", "/.", "prusecik"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"x2", "=", 
         RowBox[{
          RowBox[{"p2", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "+", "rad2"}]}], ";", 
        "\[IndentingNewLine]", "\[IndentingNewLine]", 
        RowBox[{"(*", 
         RowBox[{
         "O\[SHacek]et\[RHacek]en\[IAcute]", " ", "krajn\[IAcute]ch", " ", 
          "podm\[IAcute]nek", " ", "a", " ", "return", " ", 
          "v\[YAcute]sledku"}], "*)"}], "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"rad1", "+", "rad2"}], "<", 
            RowBox[{"Sqrt", "[", 
             RowBox[{
              SuperscriptBox[
               RowBox[{"p2", "[", 
                RowBox[{"[", "1", "]"}], "]"}], "2"], "+", 
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"p1", "[", 
                  RowBox[{"[", "2", "]"}], "]"}], "-", 
                 RowBox[{"p2", "[", 
                  RowBox[{"[", "2", "]"}], "]"}]}], ")"}], "2"]}], "]"}]}], "||", 
           RowBox[{"prusecik", "==", 
            RowBox[{"{", "}"}]}]}], ",", 
          RowBox[{"x1", "=", "rad1"}]}], "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"x2", "-", "x1"}], ">=", "0"}], ",", 
          RowBox[{"w", "=", 
           RowBox[{"x2", "-", "x1"}]}], ",", 
          RowBox[{"w", "=", "0."}]}], "]"}]}]}], "\[IndentingNewLine]", 
      "\[IndentingNewLine]", "]"}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.9529352798718643`*^9, 3.9529353369978857`*^9}, {
   3.952935378287838*^9, 3.9529354797869053`*^9}, {3.9529355107741776`*^9, 
   3.9529355856657996`*^9}, {3.952935676384483*^9, 3.952935832631583*^9}, {
   3.95293590980945*^9, 3.9529359200581665`*^9}, {3.952935970788432*^9, 
   3.952936395174972*^9}, {3.9529366168481197`*^9, 3.9529368276931295`*^9}, {
   3.9529368669672813`*^9, 3.9529369350006475`*^9}, {3.9529370057514486`*^9, 
   3.95293707321174*^9}, {3.952937394441228*^9, 3.952937441067309*^9}, {
   3.952937616489936*^9, 3.952937619905262*^9}, {3.95293774624911*^9, 
   3.9529377585378733`*^9}, {3.9529377893624754`*^9, 3.952937816007844*^9}, 
   3.9529381512734313`*^9, {3.9529382191782064`*^9, 3.952938445275852*^9}, {
   3.9529384984945936`*^9, 3.9529385099788246`*^9}, 3.9529385718210335`*^9, {
   3.9529387477882605`*^9, 3.9529388103361716`*^9}, {3.9529388970392036`*^9, 
   3.952938905364169*^9}, {3.95293893587106*^9, 3.9529389812262774`*^9}, {
   3.9529390793229337`*^9, 3.952939080634*^9}, {3.952939162972596*^9, 
   3.9529393784927917`*^9}, {3.9529394657230654`*^9, 3.9529395738259*^9}, {
   3.95293961604317*^9, 3.9529397092879868`*^9}, {3.9530007515728397`*^9, 
   3.953000777820141*^9}, {3.953000835200285*^9, 3.9530008419500885`*^9}, 
   3.9530009343997097`*^9, {3.953001153945405*^9, 3.9530011860803127`*^9}, {
   3.9530012235176363`*^9, 3.9530012593151255`*^9}, {3.953001291571084*^9, 
   3.953001334365924*^9}, {3.9530013793305893`*^9, 3.953001423274746*^9}, 
   3.953001455475423*^9, 3.95300156944923*^9, {3.9530016266546288`*^9, 
   3.9530016619309654`*^9}, {3.953001740116957*^9, 3.9530017572437954`*^9}, {
   3.953001792623351*^9, 3.9530017952872906`*^9}, {3.9530018407000275`*^9, 
   3.953001843317258*^9}, {3.9530018983861294`*^9, 3.95300202053675*^9}, {
   3.953002896097517*^9, 3.9530029062275343`*^9}, {3.953002999267458*^9, 
   3.953003010392926*^9}, {3.9530030755624647`*^9, 3.953003097152747*^9}, {
   3.953003200346046*^9, 3.9530032991217394`*^9}, {3.9530033643900037`*^9, 
   3.9530033816785936`*^9}, {3.953003418698546*^9, 3.9530034626615677`*^9}, {
   3.953003697292552*^9, 3.953003711385458*^9}, {3.953003795597426*^9, 
   3.953003819307809*^9}, {3.953003855873612*^9, 3.953003876839399*^9}, {
   3.95300393738096*^9, 3.9530041364021354`*^9}, {3.953004203541115*^9, 
   3.9530042070261183`*^9}, {3.9530042728162217`*^9, 3.953004274613542*^9}, {
   3.9530043977096343`*^9, 3.95300439924273*^9}, {3.9530049171081123`*^9, 
   3.9530050089313774`*^9}, {3.9530051203565636`*^9, 
   3.9530051395915747`*^9}, {3.953005169680179*^9, 3.9530053265717278`*^9}, {
   3.9530053724465585`*^9, 3.953005410704521*^9}, {3.9530054438690863`*^9, 
   3.95300546044202*^9}, {3.9530055208155966`*^9, 3.9530055363060865`*^9}, {
   3.9530055799622974`*^9, 3.953005944609658*^9}, {3.953006060567442*^9, 
   3.953006176151255*^9}, {3.953045027456458*^9, 3.9530452291560903`*^9}, 
   3.9530452665350275`*^9, {3.953045324189503*^9, 3.9530453882548323`*^9}, {
   3.9530454215585723`*^9, 3.953045427716935*^9}, {3.9530454780181847`*^9, 
   3.9530454841174617`*^9}, 3.953045522181984*^9, {3.953045571846213*^9, 
   3.9530456187577667`*^9}, {3.9530456650035915`*^9, 3.953045683400284*^9}, {
   3.953045721246173*^9, 3.9530457712621984`*^9}, {3.953045849177559*^9, 
   3.9530460602867436`*^9}, {3.9530461070763483`*^9, 
   3.9530461211943235`*^9}, {3.9530461872791505`*^9, 
   3.9530461951331873`*^9}, {3.953046253197362*^9, 3.9530462597107*^9}, {
   3.9530463487768188`*^9, 3.9530464823410473`*^9}, {3.9530465505896263`*^9, 
   3.953046551453376*^9}, {3.9530474629579287`*^9, 3.953047463431675*^9}, {
   3.9530481270620236`*^9, 3.953048133681076*^9}, 3.953048192800478*^9, {
   3.9530493343736663`*^9, 3.953049338651176*^9}, {3.9530494204848633`*^9, 
   3.953049481811053*^9}, {3.9530506089530225`*^9, 3.953050616709542*^9}, {
   3.9530506775172873`*^9, 3.9530506816812325`*^9}, {3.953050725778658*^9, 
   3.9530507440892878`*^9}, {3.953050800753995*^9, 3.9530508631453605`*^9}, {
   3.9530508954729586`*^9, 3.953050927630884*^9}, 3.953051583068654*^9, {
   3.953052577364326*^9, 3.9530526002725525`*^9}, {3.953052655177138*^9, 
   3.953052681076871*^9}, {3.9530535831721373`*^9, 3.953053583923006*^9}, {
   3.953054027909335*^9, 3.9530540683847246`*^9}, {3.953095538460096*^9, 
   3.9530955807269135`*^9}, {3.9530956942536354`*^9, 3.953095732234138*^9}, {
   3.953095766563799*^9, 3.9530958322921762`*^9}, {3.9530958691928186`*^9, 
   3.9530959110007463`*^9}, {3.9530959689323874`*^9, 3.953096015008075*^9}, {
   3.9530962757872734`*^9, 3.9530965278272285`*^9}, {3.9530968267618093`*^9, 
   3.9530968336723003`*^9}, {3.953097172432436*^9, 3.9530971761164474`*^9}, {
   3.9530974963405046`*^9, 3.953097634512038*^9}, {3.953097807381246*^9, 
   3.9530978337563963`*^9}, {3.9531012845627317`*^9, 3.953101315038169*^9}, {
   3.9531013557277837`*^9, 3.9531014232054048`*^9}, {3.953101457827037*^9, 
   3.953101462293116*^9}, {3.9531014932611294`*^9, 3.953101550234132*^9}, {
   3.953102446221075*^9, 3.9531024552388897`*^9}, 3.95310341861605*^9, {
   3.9531036446735425`*^9, 3.9531036474532204`*^9}, {3.953118323012701*^9, 
   3.9531183275867777`*^9}, {3.9531224117110343`*^9, 3.953122412120874*^9}, {
   3.9531276675761404`*^9, 3.953127668348832*^9}, {3.9531278002919717`*^9, 
   3.9531278009564943`*^9}, {3.9531278468387547`*^9, 
   3.9531278475836515`*^9}, {3.953196024987048*^9, 3.9531960253049283`*^9}, 
   3.9531986247249455`*^9, {3.9531996731383142`*^9, 3.9531996741907997`*^9}, {
   3.9532024504191275`*^9, 3.9532026764946675`*^9}, {3.9532027695160275`*^9, 
   3.953202827208037*^9}, {3.953202860617907*^9, 3.953202996276241*^9}, {
   3.9532055627276897`*^9, 3.9532055797707663`*^9}, {3.9532116967885823`*^9, 
   3.9532117185521483`*^9}, 3.9534614122641*^9, 3.953461692396603*^9, {
   3.953463714764241*^9, 3.9534637152048473`*^9}, {3.953463973236395*^9, 
   3.953464009159013*^9}, {3.953464312090768*^9, 3.953464313513056*^9}, {
   3.9534704617835054`*^9, 3.9534704874781446`*^9}, {3.953470558311071*^9, 
   3.953470578985241*^9}, 3.953471360939825*^9, {3.9534716245472183`*^9, 
   3.9534716278531065`*^9}, {3.953550870750945*^9, 3.953550895818293*^9}, {
   3.9535510534310055`*^9, 3.9535510543181877`*^9}, 3.953551096466219*^9, {
   3.9535517351519175`*^9, 3.953551748218913*^9}, {3.9545701557099867`*^9, 
   3.954570210366741*^9}, {3.954570243864868*^9, 3.9545703540595927`*^9}, {
   3.954570465342715*^9, 3.9545704743436146`*^9}, {3.954570758739609*^9, 
   3.9545707868844748`*^9}, {3.9545708172138042`*^9, 
   3.9545708230702696`*^9}, {3.9545708726617403`*^9, 3.95457089909772*^9}, {
   3.9545709742493067`*^9, 3.9545710648421106`*^9}, {3.9545711642646737`*^9, 
   3.954571229858988*^9}, {3.9545712609148865`*^9, 3.954571276057654*^9}, 
   3.9545713881657677`*^9, {3.95457141929935*^9, 3.9545714384604454`*^9}, {
   3.9545715827323713`*^9, 3.954571768137377*^9}, {3.954572062939503*^9, 
   3.954572073339546*^9}, {3.954572224299217*^9, 3.9545722259046946`*^9}, {
   3.9545724948498793`*^9, 3.954572495385999*^9}, {3.9545782400581203`*^9, 
   3.9545783077818265`*^9}, {3.954578384280878*^9, 3.95457839147386*^9}, {
   3.9545798680991535`*^9, 3.9545798715100136`*^9}, {3.954579998500524*^9, 
   3.954580053891636*^9}, {3.954580330258462*^9, 3.9545803314447656`*^9}, {
   3.9545804291004386`*^9, 3.954580429588647*^9}, {3.9545838707438602`*^9, 
   3.954583999296364*^9}, {3.9545840517484655`*^9, 3.954584102828826*^9}, {
   3.95458416033027*^9, 3.9545843054593563`*^9}, {3.9545843406503143`*^9, 
   3.954584355060046*^9}, {3.954584702723419*^9, 3.9545847180843334`*^9}, {
   3.9545862812132096`*^9, 3.9545863115105734`*^9}, {3.954586399409381*^9, 
   3.954586400018473*^9}, {3.95458663600809*^9, 3.9545866507832303`*^9}, {
   3.9545869599685698`*^9, 3.9545869606548047`*^9}, {3.9545882086999426`*^9, 
   3.9545882137633085`*^9}, {3.9545883032390594`*^9, 3.954588321461116*^9}, {
   3.9545886532110567`*^9, 3.9545886537271023`*^9}, {3.9545949372707376`*^9, 
   3.9545949427055273`*^9}, {3.95464451619975*^9, 3.9546445170146933`*^9}, {
   3.9546447940260005`*^9, 3.9546448109730153`*^9}, 3.9546461220866885`*^9, 
   3.9546463515935497`*^9, {3.9546537145054197`*^9, 3.954653739123309*^9}, {
   3.954653783415163*^9, 3.9546539238804502`*^9}, {3.954653975999436*^9, 
   3.9546539859925547`*^9}, {3.9546540380293503`*^9, 
   3.9546540802958107`*^9}, {3.9546541475462246`*^9, 
   3.9546541586460457`*^9}, {3.954654522859562*^9, 3.954654531825238*^9}, {
   3.954656827081189*^9, 3.9546568599862175`*^9}, {3.9546568942017827`*^9, 
   3.9546568947388983`*^9}, {3.954657430595415*^9, 3.9546574440226655`*^9}, 
   3.9546754411086364`*^9, {3.954676367307604*^9, 3.9546763769329667`*^9}, {
   3.9563064594364605`*^9, 3.9563064717662344`*^9}, {3.9563067379308577`*^9, 
   3.95630686122511*^9}, {3.9563069067387404`*^9, 
   3.956306988418821*^9}},ExpressionUUID->"9c8230d8-0396-42f1-8292-\
ebb93a21922d"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"V\[YAcute]po\[CHacek]et", " ", "SFR"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"SFR1", "=", 
     RowBox[{"0.2", "*", "Ng", "*", 
      RowBox[{"NIntegrate", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Dc1", "[", 
          RowBox[{"i", ",", 
           RowBox[{"p1", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"p2", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "]"}], "*", 
         RowBox[{"cc", "[", "i", "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "ccmin", ",", "ccmax"}], "}"}], ",", 
        RowBox[{"MaxRecursion", "\[Rule]", "20"}], ",", 
        RowBox[{"Method", "\[Rule]", "\"\<LocalAdaptive\>\""}]}], "]"}]}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{"SFR", "=", "SFR1"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Ur\[CHacek]en\[IAcute]", " ", "po\[CHacek]tu", " ", "blesk\[URing]", " ",
      "do", " ", "zemn\[IAcute]c\[IAcute]ho", " ", "lana"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"pocetbleskuDoZemLana", "=", 
    RowBox[{"0.2", "*", "Ng", "*", 
     RowBox[{"NIntegrate", "[", 
      RowBox[{
       RowBox[{"0.67", "*", 
        SuperscriptBox[
         RowBox[{"p1", "[", 
          RowBox[{"[", "2", "]"}], "]"}], "0.45"], "*", 
        SuperscriptBox["i", "0.74"], "*", 
        RowBox[{"cc", "[", "i", "]"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "ccmin", ",", "ccmax"}], "}"}], ",", 
       RowBox[{"MaxRecursion", "\[Rule]", "20"}], ",", 
       RowBox[{"Method", "\[Rule]", "\"\<LocalAdaptive\>\""}]}], 
      "]"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.954580063620654*^9, 3.9545801177087812`*^9}, {
   3.9545844790067816`*^9, 3.954584512264228*^9}, {3.954584557339294*^9, 
   3.9545845643247395`*^9}, {3.954584797325539*^9, 3.954584818716674*^9}, {
   3.95458841152674*^9, 3.954588499614276*^9}, {3.954654968762662*^9, 
   3.954654970656413*^9}, {3.9546763818356857`*^9, 3.954676384355595*^9}, 
   3.9547377396998787`*^9, {3.956306998531057*^9, 
   3.956307016739554*^9}},ExpressionUUID->"c4b95d54-3fb0-4497-b093-\
40e5be6b2b96"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   FractionBox[
    RowBox[{
     RowBox[{"ArcTan", "[", 
      FractionBox[
       RowBox[{"p2", "[", 
        RowBox[{"[", "1", "]"}], "]"}], 
       RowBox[{
        RowBox[{"p1", "[", 
         RowBox[{"[", "2", "]"}], "]"}], "-", 
        RowBox[{"p2", "[", 
         RowBox[{"[", "2", "]"}], "]"}]}]], "]"}], "*", "180"}], "Pi"], "//", 
   "N"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "Ur\[CHacek]en\[IAcute]", " ", "ochrann\[EAcute]ho", " ", "\[UAcute]hlu"}],
    "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Sqrt", "[", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"p2", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "2"], "+", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"p1", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "-", 
       RowBox[{"p2", "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], ")"}], "2"]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{
   "Absolutn\[IAcute]", " ", "vzd\[AAcute]lenost", " ", "mezi", " ", 
    "ochrann\[YAcute]m", " ", "a", " ", "chr\[AAcute]n\[EHacek]n\[YAcute]m", 
    " ", "vodi\[CHacek]em"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.9545725545299516`*^9, 3.9545725906333413`*^9}, {
   3.9545726719959326`*^9, 3.95457270173346*^9}, 3.9545727433823643`*^9, {
   3.9545727735320196`*^9, 3.954572778578215*^9}, {3.9545779856689916`*^9, 
   3.954578066893379*^9}, {3.9545781883439827`*^9, 3.9545782104072633`*^9}, {
   3.9545802417158413`*^9, 3.9545802724187994`*^9}, {3.954582098050541*^9, 
   3.954582132342205*^9}, {3.9563070294066124`*^9, 
   3.956307037993352*^9}},ExpressionUUID->"acb6071b-46ce-4b4d-826e-\
3ee13100f3cc"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Export", " ", "dat", " ", "a", " ", "plot"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"dc1proud", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"Dc1", "[", 
          RowBox[{"i", ",", 
           RowBox[{"p1", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"p2", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "2", ",", "40", ",", "0.1"}], "}"}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{"Export", "[", 
    RowBox[{"\"\<data.xlsx\>\"", ",", "dc1proud"}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"ListPlot", "[", "dc1proud", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.9545818888102818`*^9, 3.954582061052946*^9}, {
   3.954582177483175*^9, 3.954582239775855*^9}, {3.954582280186812*^9, 
   3.9545823078023453`*^9}, {3.954586430242647*^9, 3.954586432649931*^9}, {
   3.9545865507373233`*^9, 3.954586581399579*^9}, 3.95458667734653*^9, 
   3.9545868851248784`*^9, 3.954588229038267*^9, {3.954588290240892*^9, 
   3.954588297243288*^9}, {3.9546448566250772`*^9, 3.9546449067336135`*^9}, {
   3.9546456352594943`*^9, 3.9546457054323716`*^9}, {3.9546457425696735`*^9, 
   3.954645749175644*^9}, {3.9546459610329547`*^9, 3.9546459938080287`*^9}, {
   3.954646198099763*^9, 3.954646199213978*^9}, {3.9546462296792555`*^9, 
   3.9546462920251603`*^9}, {3.954646943763897*^9, 3.9546469540909367`*^9}, {
   3.954647949020691*^9, 3.9546479499836326`*^9}, {3.9546763908162355`*^9, 
   3.954676393237277*^9}, {3.95630705306143*^9, 
   3.956307064350825*^9}},ExpressionUUID->"db1e32d3-3893-41b5-b7da-\
2f27ce49ba5f"]
},
WindowSize->{574.8, 573.6},
WindowMargins->{{Automatic, -4.7999999999999545`}, {Automatic, 0}},
Magnification:>1.1 Inherited,
FrontEndVersion->"12.3 for Microsoft Windows (64-bit) (July 9, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"230a5cbd-8d56-4519-b98b-ca85356d3744"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 16901, 327, 1035, "Input",ExpressionUUID->"9c8230d8-0396-42f1-8292-ebb93a21922d"],
Cell[17462, 349, 2194, 52, 247, "Input",ExpressionUUID->"c4b95d54-3fb0-4497-b093-40e5be6b2b96"],
Cell[19659, 403, 1640, 44, 131, "Input",ExpressionUUID->"acb6071b-46ce-4b4d-826e-3ee13100f3cc"],
Cell[21302, 449, 1779, 38, 115, "Input",ExpressionUUID->"db1e32d3-3893-41b5-b7da-2f27ce49ba5f"]
}
]
*)