{VERSION 6 0 "Windows 2000" "6.0" } {USTYLETAB {PSTYLE "Heading 4" -1 20 1 {CSTYLE "" -1 -1 "Times" 1 10 0 0 0 1 1 1 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Ord ered List 5" -1 200 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 144 2 0 2 2 -1 1 }{PSTYLE "Ordered List 1" -1 201 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Bullet Item" -1 15 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Lef t Justified Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Help" -1 10 1 {CSTYLE "" -1 -1 "Courier" 1 9 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Author" -1 19 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 8 8 2 0 2 0 2 2 -1 1 }{PSTYLE "Diagnostic" -1 9 1 {CSTYLE "" -1 -1 "Courier" 1 10 64 128 64 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Headi ng 2" -1 4 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }1 1 0 0 8 2 2 0 2 0 2 2 -1 1 }{PSTYLE "Ordered List 3" -1 202 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 72 2 0 2 2 -1 1 }{PSTYLE "Maple Plot" -1 13 1 {CSTYLE "" -1 -1 "Tim es" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Dash Item" -1 16 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 3" -1 5 1 {CSTYLE "" -1 -1 "Times " 1 12 0 0 0 1 1 1 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 } {PSTYLE "Ordered List 4" -1 203 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 108 2 0 2 2 -1 1 }{PSTYLE "Maple \+ Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Line Printed Output" -1 6 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "List Item" -1 14 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Warning" -1 7 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Error" -1 8 1 {CSTYLE "" -1 -1 "Courier" 1 10 255 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Fixed Width" -1 17 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }1 1 0 0 8 4 2 0 2 0 2 2 -1 1 }{PSTYLE "Title " -1 18 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 1 2 2 2 1 0 0 1 }3 1 0 0 12 12 2 0 2 0 2 2 -1 1 }{PSTYLE "Ordered List 2" -1 204 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 36 2 0 2 2 -1 1 }{CSTYLE "Help Variable" -1 25 "Courier" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Text" -1 200 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Bold" -1 39 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Page Number" -1 33 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic Small" -1 201 "Times" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Nonterminal" -1 24 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Default" -1 38 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Comment" -1 21 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Input" -1 0 "Courier" 1 12 255 0 0 1 2 1 2 2 1 2 0 0 0 1 }{CSTYLE "2D Math Small" -1 7 "Times " 1 1 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Fixed" -1 23 "Courier " 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Popup" -1 31 "Times" 1 12 0 128 128 1 1 2 1 2 2 2 0 0 0 1 }{CSTYLE "Plot Title" -1 27 "Times" 1 10 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Input" -1 19 "Times" 1 12 255 0 0 1 2 2 2 2 1 2 0 0 0 1 }{CSTYLE "Copyright" -1 34 "Times" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Input Placeholder" -1 202 "Courier" 1 12 200 0 200 1 2 1 2 2 1 2 0 0 0 1 }{CSTYLE "2D Math Bold Small" -1 10 "Times" 1 1 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Mat h" -1 2 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Notes" -1 37 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underli ned Bold" -1 41 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 0 0 0 1 }{CSTYLE "Hyp erlink" -1 17 "Times" 1 12 0 128 128 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "2 D Math Symbol 2" -1 16 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 } {CSTYLE "Plot Text" -1 28 "Times" 1 8 0 0 0 1 2 2 2 2 2 2 0 0 0 1 } {CSTYLE "Help Italic" -1 42 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 } {CSTYLE "Output Labels" -1 29 "Times" 1 8 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Heading" -1 26 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Normal" -1 30 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Comment" -1 18 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Maple Name" -1 35 "Times" 1 12 104 64 92 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Output" -1 20 "Times" 1 12 0 0 255 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Dictionary Hyperlink" -1 45 "Times" 1 12 147 0 15 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "Help Emphasized" -1 22 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Italic Bold" -1 40 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "LaTeX" -1 32 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Menus" -1 36 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Prompt" -1 1 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined" -1 44 "Times" 1 12 0 0 0 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined Italic" -1 43 "Times" 1 12 0 0 0 1 1 2 1 2 2 2 0 0 0 1 }{CSTYLE "2D Math Bold" -1 5 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic" -1 3 "Times " 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{PSTYLE "" -1 205 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 2 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 99 "restart: with(Statis tics): with(plots): with(plottools): with(Maplets): with(Maplets[Eleme nts]): \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 22 "PlotOptions:= proc()\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 23 "local \+ r1,jadro,miesz:\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 48 " if Maplets: -Tools:-Get(Roz1_1) = 'true' then \n" }{MPLTEXT 1 0 120 " r1:=Normal (Maplets:-Tools:-Get(m::realcons),Maplets:-Tools:-Get(s::positive)); \+ miesz:=0; \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 107 " elif Maplets:-Tools:-Get(Roz1_2) = 'true' then r1:=Ex ponential(Maplets:-Tools:-Get(l::positive));\n" }{MPLTEXT 1 0 12 "mies z:=0; \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 137 " elif Maplets:-Tools :-Get(Roz1_3) = 'true' then r1:=Uniform(Maplets:-Tools:-Get(a: :realcons),Maplets:-Tools:-Get(b::realcons));\n" }{MPLTEXT 1 0 11 "mie sz:=0;\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 104 "elif Maplets:-Tools: -Get(Roz1_5) = 'true' then r1:=ChiSquare(Maplets:-Tools:-Get(n 1::posint));\n" }{MPLTEXT 1 0 12 " miesz:=0;\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 138 "elif Maplets:-Tools:-Get(Roz1_6) = 'true' then \+ r1:=Weibull(Maplets:-Tools:-Get(t1::realcons),Maplets:-Tools:-Get(t 2::realcons));\n" }{MPLTEXT 1 0 11 "miesz:=0;\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 9 "end if;\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 56 "if \+ Maplets:-Tools:-Get(Jadro) = \" Gaussowskie\" then \n" }{MPLTEXT 1 0 13 " jadro:=1;\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 85 " elif Maplet s:-Tools:-Get(Jadro) = \" Epanechnikowa (optymalne)\" then jadro:=2 ;\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 67 " elif Maplets:-Tools:-Get( Jadro) = \" Biwagowe\" then jadro:=3;\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 67 "elif Maplets:-Tools:-Get(Jadro) = \" Trojkatne\" then jadro:=4;\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 69 "elif Maplets:-T ools:-Get(Jadro) = \" Prostokatne\" then jadro:=5;\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 10 "end if; \n" }{MPLTEXT 1 0 97 " est_jadrowa( r1,jadro,Maplets:-Tools:-Get(N::integer),Maplets:-Tools:-Get(h::positi ve),miesz,\n" }{MPLTEXT 1 0 192 " Maplets:-Tools:-Get(e::positive) ,Maplets:-Tools:-Get(m11::realcons), Maplets:-Tools:-Get(s11::posi tive),Maplets:-Tools:-Get(m22::realcons), Maplets:-Tools:-Get(s22: :positive) );\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 11 "end proc:\n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 31 "Estymacja_jadrowa := Maplet( \n" }{MPLTEXT 1 0 138 " Window[W1](title=\"Estymacja jadrowa \25 6 by P.Koscielniak, J.Ombach and J.Szczepanski\", height=650, width=73 0, menubar=MB1, \n" }{MPLTEXT 1 0 24 " 'layout' = BL1),\n" } {MPLTEXT 1 0 175 " MenuBar[MB1](Menu(\"Zamknij\", MenuItem(\"Zamknij \", Shutdown() ) ),Menu(\"Parametry\", MenuItem(\"Przywroc\", 'onclick '='A3' ), MenuItem(\"Wyczysc\", 'onclick'='A2' ) )),\n" }{MPLTEXT 1 0 4 " \n" }{MPLTEXT 1 0 22 "BoxLayout[BL1]( \n" }{MPLTEXT 1 0 7 " \+ \n" }{MPLTEXT 1 0 66 " BoxColumn(inset=0, spacing=0, background= COLOR(RGB,1,1,4/5),\n" }{MPLTEXT 1 0 83 " BoxRow(Label(\"Licznosc probki N=\", 'font' = Font(\"default\", bold, 12)),\n" }{MPLTEXT 1 0 102 " TextField[N](\"100\",7),Label(\"Szerokosc pasma h=\", 'fo nt' = Font(\"default\", bold, 12)),\n" }{MPLTEXT 1 0 69 " TextF ield[h](\"0.9\",4),Button(\"Wykonaj\", 'onclick'='A1',\n" }{MPLTEXT 1 0 69 " background=COLOR(RGB,209/255,209/255,83/85 ), \n" }{MPLTEXT 1 0 63 " 'font' = Font(\"default \", bold, 12))),\n" }{MPLTEXT 1 0 3 " \n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 14 "BoxRow( \n" }{MPLTEXT 1 0 128 "BoxColumn(Label(\" Rozklad, z ktorego losujemy probke\",'font' = Font(\"default\", bold, \+ 12)), background=COLOR(RGB,1,1,4/5), \n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 235 "BoxRow(RadioButton['Roz1_1'](\"Normalny N(m,s)\", 'v alue'='true','group'='BG1'), Label(\" m=\", 'font' =Font(\"default\", \+ bold, 12)),TextField[m](\"0\",4), Label(\" s=\", 'font' = Font(\"defau lt\", bold, 12)), TextField[s](\"1\",4)),\n" }{MPLTEXT 1 0 3 " \n" } {MPLTEXT 1 0 161 "BoxRow(RadioButton['Roz1_2'](\"Wykladniczy o paramet rze l\", 'group'='BG1'), Label(\" l=\", 'font' =Font(\"default \", bold, 12)),TextField[l](\"2\",4)),\n" }{MPLTEXT 1 0 4 " \n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 274 "BoxRow( RadioButton['Roz1_3']( \"Jednastajny na odcinku [a,b]\", 'group'='BG1'), Label(\" a=\", 'font ' =Font(\"default\", bold, 12)),TextField[a](\"0\",4),Label(\" b=\", ' font' = Font(\"default\", bol d, 12)),TextField[b](\"1\",4)), \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 145 " BoxRow(RadioButton['Roz1_5'](\"Chi-kwadrat[n1]\",'group'='BG1'), Label (\" n1=\", 'font' =Font(\"default\", bold, 12)),TextField[n1](\"3\",4) ), \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 224 "BoxRow(RadioButton['Roz 1_6'](\"Weibulla[t1,t2]\",'group'='BG1'), Label(\" t1=\", 'font' =Font (\"default\", bold, 12)),TextField[t1](\"1\",4), Label(\" t2=\", 'font ' = Font(\"default\", bold, 12)), TextField[t2](\"2\",4)), \n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 162 "BoxRow(Rad ioButton['Roz1_7'](\"Mieszany eN(m1,s1)+(1-e)N(m2,s2)\",'group'='BG1') ,Label(\" e=\", 'font' =Font(\"default\", bold, 12)),TextField[e](\"0. 5\",4) ), \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 327 "BoxRow( Label(\" m1=\", 'font' =Font(\"default\", bold, 12)),TextField[m11](\"0\",4), \+ Label(\" s1=\", 'font' = Font(\"default\", bold, 12)), TextField[s11]( \"1\",4), Label(\" m2=\", 'font' =Font(\"default\", bold, 12)),TextFie ld[m22](\"5\",4), Label(\" s2=\", 'font' = Font(\"default\", bold, 12) ), TextField[s22](\"1\",4) ) \n" }{MPLTEXT 1 0 4 " \n" }{MPLTEXT 1 0 5 "), \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 " \n" }{MPLTEXT 1 0 70 " BoxColumn( Label(\"Jadro\", 'font' = Font(\"def ault\", bold, 12)), \n" }{MPLTEXT 1 0 63 " DropDow nBox['Jadro']([\" Gaussowskie\", \n" }{MPLTEXT 1 0 76 " \+ \" Epanechnikowa (optymalne)\", \n" } {MPLTEXT 1 0 59 " \" Biwagowe \", \n" }{MPLTEXT 1 0 59 " \" Trojkatne\",\n" }{MPLTEXT 1 0 95 " \+ \" Prostokatne\"]), background=COLOR(RGB,1,1,2/5), \n" } {MPLTEXT 1 0 107 "BoxRow(Label(\"Bledy dla wybranego i optymalnego jad ra wynosza\",'font' = Font(\"default\", bold, 12))), \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 16 "Plotter[P1]() \n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 5 " ))\n" } {MPLTEXT 1 0 17 " )),\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 23 "ButtonGroup['BG1'](),\n" }{MPLTEXT 1 0 2 " \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 47 "Action['A1']( Evaluate('P1' ='PlotOptions')), \n" }{MPLTEXT 1 0 57 "Action['A2']( SetOption('N'=\" \"),SetOption('h'=\"\")),\n" }{MPLTEXT 1 0 62 "Action['A3']( SetOption ('N'=\"100\"),SetOption('h'=\"0.9\"))\n" }{MPLTEXT 1 0 4 "):\n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 85 "########################### END \+ OF MAPLET DEFINITION ##############################\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 20 "jadro_opt:=proc(x)\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 40 "if x < -evalf(sqrt(5)) then return(0);\n" }{MPLTEXT 1 0 41 "elif x > evalf(sqrt(5)) then return(0);\n" }{MPLTEXT 1 0 51 "e lse return( evalf(3/(4*sqrt(5))*(1-(x^2)/5)) );\n" }{MPLTEXT 1 0 10 " end if;\n" }{MPLTEXT 1 0 11 "end proc:\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 20 "jadro_Biw:= proc(x)\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 27 "if x < -1 then retur n(0);\n" }{MPLTEXT 1 0 28 "elif x > 1 then return(0);\n" }{MPLTEXT 1 0 34 "else return( 15/16*(1-x^2)^2 );\n" }{MPLTEXT 1 0 10 " end if;\n " }{MPLTEXT 1 0 11 "end proc:\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 21 "jadro_troj:=proc(x)\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 27 "if x < -1 then return(0);\n" }{MPLTEXT 1 0 28 "elif x > 1 then return(0);\n" }{MPLTEXT 1 0 27 "else return( 1-abs(x) );\n" }{MPLTEXT 1 0 10 " end if;\n" }{MPLTEXT 1 0 11 "end proc:\n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 22 "jadro_prost :=proc(x)\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 27 "if x < -1 then ret urn(0);\n" }{MPLTEXT 1 0 28 "elif x > 1 then return(0);\n" }{MPLTEXT 1 0 22 "else return( 1/2 );\n" }{MPLTEXT 1 0 10 " end if;\n" } {MPLTEXT 1 0 11 "end proc:\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 64 "e st_jadrowa := proc(rozklad,jadro,N,h,miesz,e,m11,s11,m22,s22)\n" } {MPLTEXT 1 0 3 " \n" }{MPLTEXT 1 0 82 "local dane,dane_sort,X,a,b,K,ff ,ff_opt,punkty,f,bb,blad,gestosc,bb_opt,blad_opt:\n" }{MPLTEXT 1 0 2 " \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 69 "gestosc:= t->e*PDF(Normal(m 11,s11),t)+(1-e)*PDF(Normal(m22,s22),t):\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 18 "if miesz=0 then \n" }{MPLTEXT 1 0 29 "X:=RandomVariab le(rozklad):\n" }{MPLTEXT 1 0 23 "f:=x->PDF(rozklad,x):\n" }{MPLTEXT 1 0 7 "else \n" }{MPLTEXT 1 0 33 "X:=RandomVariable(PDF=gestosc):\n" } {MPLTEXT 1 0 14 "f:=gestosc: \n" }{MPLTEXT 1 0 10 "end if: \n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 20 "dane:=Sample(X,N):\n" }{MPLTEXT 1 0 24 "dane_sort:=sort(dane):\n" }{MPLTEXT 1 0 20 "a:=dane_sort[1]-2: \n" }{MPLTEXT 1 0 20 "b:=dane_sort[N]+2:\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 49 "if jadro=1 then K:=x->exp(-(x^2)/2)/sqrt(2*Pi);\n" } {MPLTEXT 1 0 39 "elif jadro=2 then K:=x->jadro_opt(x);\n" }{MPLTEXT 1 0 39 "elif jadro=3 then K:=x->jadro_Biw(x);\n" }{MPLTEXT 1 0 41 "elif \+ jadro=4 then K:=x->jadro_troj(x); \n" }{MPLTEXT 1 0 41 "elif jadro=5 t hen K:=x->jadro_prost(x);\n" }{MPLTEXT 1 0 9 "end if;\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 50 "ff := x ->1/(N*h)*add( K((x-dane[i])/h),i =1..N):\n" }{MPLTEXT 1 0 62 "ff_opt := x ->1/(N*h)* add(jadro_opt((x-dane[i])/h),i =1..N):\n" }{MPLTEXT 1 0 2 "\n" } {MPLTEXT 1 0 28 "bb := x->(ff(x) - f(x))^2:\n" }{MPLTEXT 1 0 37 "bb_op t := x->(ff_opt(x) - f(x))^2: \n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 38 "blad:=evalf(Int(bb,a..b,5,_NCrule));\n" }{MPLTEXT 1 0 46 "blad_opt :=evalf(Int(bb_opt,a..b,5,_NCrule));\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 83 "plots[display](\{plot(ff,a..b,title = ca t(\"\",convert([blad,blad_opt],string))),\n" }{MPLTEXT 1 0 46 " \+ plot(f,a..b,color=green)\});\n" }{MPLTEXT 1 0 11 "end proc:\n " }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }{MPLTEXT 1 0 2 "\n" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 7 "" 1 "" {TEXT 203 51 "Wa rning, the name changecoords has been redefined\n" }}{PARA 7 "" 1 "" {TEXT 203 59 "Warning, the assigned name arrow now has a global bindin g\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "Display(Estymacja_ja drowa);" }}{PARA 205 "" 1 "" {TEXT 200 38 "Initializing Java runtime e nvironment." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {MARK "0 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }