| Current File : //proc/thread-self/root/kunden/usr/share/emacs/27.2/lisp/calc/calcalg3.elc |
;ELC
;;; Compiled
;;; in Emacs version 27.2
;;; with all optimizations.
;;; This file uses dynamic docstrings, first added in Emacs 19.29.
;;; This file does not contain utf-8 non-ASCII characters,
;;; and so can be loaded in Emacs versions earlier than 23.
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
(byte-code "\300\301!\210\300\302!\207" [require calc-ext calc-macs] 2)
#@59 Apply BINOP to the elements of the lists ARGS1 and ARGS2.
(defalias 'math-map-binop #[(binop args1 args2) "\205 @\n@\"\303 A\nA#B\207" [args1 binop args2 math-map-binop] 5 (#$ . 480)])
(defalias 'calc-find-root #[(var) "\300\301`\"\207" [calc-do #[nil "\302 \203 \303\202\n \304 \305\232\204 \306\232\203) \307\310\311\312\313!\312\314!\312\310!F#\202` \315\316 \"\203@ \315\317 \"\204@ \320\321 \322Q!\202C \320 !\211\242\323=\203R \323\324 A@\"\210\307\314\311\312\310! \312\314!F#))\207" [func var calc-is-hyperbolic calcFunc-wroot calcFunc-root "" "$" calc-enter-result 2 "root" calc-top-n 3 1 string-match ",\\|[^ ] +[^ ]" "\\[" math-read-expr "[" "]" error "Bad format in expression: %s"] 9]] 3 nil "sVariable(s) to solve for: "])
(defalias 'calc-find-minimum #[(var) "\300\301`\"\207" [calc-do #[nil "\303 \203 \304 \203 \305\202 \306\202 \304 \203 \307\202 \310\303 \203% \311\202& \312\n\313\232\2044 \n\314\232\203F \315\316 \317\320!\317\321!\317\316!F#\202} \322\323\n\"\203] \322\324\n\"\204] \325\326\n\327Q!\202` \325\n!\211\242\330=\203o \330\331\nA@\"\210\315\321 \317\316!\n\317\321!F#)*\207" [tag func var calc-is-inverse calc-is-hyperbolic calcFunc-wmaximize calcFunc-maximize calcFunc-wminimize calcFunc-minimize "max" "min" "" "$" calc-enter-result 2 calc-top-n 3 1 string-match ",\\|[^ ] +[^ ]" "\\[" math-read-expr "[" "]" error "Bad format in expression: %s"] 9]] 3 nil "sVariable(s) to minimize over: "])
(defalias 'calc-find-maximum #[(var) "\301 \210\302!\207" [var calc-invert-func calc-find-minimum] 2 nil "sVariable to maximize over: "])
(defalias 'calc-poly-interp #[(arg) "\300\301`\"\207" [calc-do #[nil "\302\303! :\204 \304=\204 \303=\203 \305\306\303\211\"B\202' \203' \307\310!\210\311 \2039 \312\313\314\315\302\313!E#\202C \312\313\316\317\302\313!E#)\207" [data arg calc-top 2 0 vec calc-top-list error "Bad prefix argument" calc-is-hyperbolic calc-enter-result 1 "rati" calcFunc-ratint "poli" calcFunc-polint] 7]] 3 nil "P"])
#@29 History for calc-curve-fit.
(defvar calc-curve-fit-history nil (#$ . 2488))
(defalias 'calc-curve-fit #[(arg &optional calc-curve-model calc-curve-coefnames calc-curve-varnames) "\300\301`\"\207" [calc-do #[nil "\306\307 \203 \310\202 \311 \203 \312\202 \313\306\314\306\211\211\211\211\211\211\315 @ A B C D EF\204r\316\317C 8A\203D \320\202E \321\n\203M \322\202N \323$\210\324 \211D\325U\203` \326 \210\2020 D\327U\203r CT G\246C\2020 D\330>\203 \n?\2020 D\331U\203\263 A\203\221 \306A\2020 \332\333!\311 \204\246 \307 \204\246 G\334U\204\254 \335A\202\257 A)\2020 D\336=\203\300 \333@\202\303 \314@G\204\325 @T@\332@!\202/G:\204\342 G\314=\203 @\337\\@\332@!\340 !\203\377 \341 \332@S!C\"\202\342 \332@S!E\202/\343G!\211G\314V\203+\342\344G@T\"B@G\\@\202/\345\346!\210\340 !\204? AA\203?\345\347!\210 G\337Z\306\211HID\350U\203m\351
\211T\n\205Z\314#\210\352H\342\333IABB\"\211F\2022 D\353Y\203\231D\354X\203\231\351\333\355D\356\357#\n\205\210\314#\210\360HAIA@\"\211F\2022 D\361U\203\275\351\333 A@GS\n\205\254\314#\210\360HAIA@\"\211F\2022 D\362U\203\351\351
\211T\n\205\315\333#\210\352HA@\363\364\365\366I\342HAAB#\"\"\211F\2022 D\367U\203\351
\211T\n\205\371\333#\210\352HA@\363\370\365\371\342HAABI#\"\"\211F\2022 D\372U\203.\373\211BF\374\375!\210\376E!\210\2020 D\377U\203I\373\211BF\374\375!\210\201X E!\210\2020 D\201Y U\203\214\373\211BF\374\375!\210A\203\202A;\204\202\342AA@\342\201Z \201[ \337A8AAA@A#BEA\201\\ E!\210\2020 D\201] \267\202;\351
\211T\n\205\236\333#\210\352HA@\363\201^ \365D\201_ =\203\271\201` \202\274\201a \365\201b \342HAABI#\"\"\"\211F\2022 \351
\211T\n\205\333\314#\210\352H\342\333IABB\"FD\201c =\203\374\201d FD\202\201e \201f FE\211F\2022 \351
\211T\n\205\314#\210\352H\342\333\365D\201g =\203+\201h \202.\201i I\"ABB\"\211F\2022 D\201j U\203\221\351
\337
_T\n\205O\314#\210HI J\211 KA@FJAJKAA\211K\203\215\201k F\201l K@\201e \355J@KA@E\337EE\"F\202^*\2020 D\201m U\203\275\201n \201o !F\201p I\201q H\351\333HGS\n\205\270\333#\210\2020 D\201r >\203j\306\211 L MD\201s =\2032N O\314\211 P Q\201t \323\201u \306\201v $\211FG\333U\204\372\345\201w !\210F@F\373MP\314V\203 \342NGPZ\201x N!\233BH\202.Q\314V\203.\342\201y Q!BH+\202\307D\201z \267\202O\201{ \201| !\202R\201{ \201} !\202R\332\333!\211F\204^\345\201~ !\210F\211\242\342=\262\203\307FG\201 >\203\301\201\200 FA@!\204\301\337F8\211\242\342=\262\203\301\334F8\203\237\334F8\211\242\342=\262\203\301\337F8I\334F8\206\266\342\201\201 FI\"BHFA@F\202\307\345\201\202 !\210I\204QF\242\201\203 = RH\203\372\351R\203\345
T\202\346
HGS\201\201 FH\"\306R%\210\202P\201\201 F\306\" S\306 T\355SG
R\203\337\202\333# @\306 U@\314W\203(\345\201\204 !\210@S\233\211UATU\306\241\210\351R\203C
T\202D
SGTSR%\210,)M\203f\201\205 \342FIHF\201\206 \"\210*\2020 \201\207 \210\2020 B\204\314\373 V\201\210 @\201\211 \201\212 E!\201\213 \"EFIG\337U\203\237IA@\202\241IHG\337U\203\260HA@\202\262H \257#\210V:\203\313\201\205 \201\214 V!\201\215 \"\210)A\205`A;\203\341\316\201\216 A\"\202`\373 W\201\217 \373!\210\201\220 \201\221 AA@!\201\221 \337A8!\"\210\201\222 \337 8!\204\201\223 \306\211\"\210\201\224 \306!\210AA@AA\201\225 \334\201\226 \201\227 A@AA\"\201\230 \"\201k \201\231 \201\232 A@AA\"\"FA\201\220 \201\221 A!\201\221 \201\233 \333!!\"\210\201\234 \306!). \207" [calc-aborted-prefix msgs homog data temp calc-curve-nvars nil calc-is-inverse calcFunc-xfit calc-is-hyperbolic calcFunc-efit calcFunc-fit 0 ("(Press ? for help)" "1 = linear or multilinear" "2-9 = polynomial fits; i = interpolating polynomial" "p = a x^b, ^ = a b^x" "e = a exp(b x), x = exp(a + b x), l = a + b ln(x)" "E = a 10^(b x), X = 10^(a + b x), L = a + b log10(x)" "q = a + b (x-c)^2" "g = (a/b sqrt(2 pi)) exp(-0.5*((x-c)/b)^2)" "s = a/(1 + exp(b (x - c)))" "b = a exp(b (x - c))/(1 + exp(b (x - c)))^2" "o = (y/x) = a (1 - x/b)" "h prefix = homogeneous model (no constant term)" "P prefix = plot result" "' = alg entry, $ = stack, u = Model1, U = Model2") message "Fit to model: %s:%s%s" "P" " " "h" "" read-char 7 keyboard-quit 63 (104 72) 80 calc-top 1 3 "Can't plot" 36 2 math-matrixp append vec prefix-numeric-value calc-top-list error "Bad prefix argument" "Data matrix is not a matrix!" 49 calc-get-fit-variables math-mul 50 57 - 48 -1 math-build-polynomial-expr 105 112 calcFunc-reduce (var mul var-mul) calcFunc-map (var pow var-pow) 94 (var mul var-mul) (var pow var-pow) 115 t require calc-nlfit calc-fit-s-shaped-logistic-curve 98 n plot nonlinear which key func calc-curve-model arg calc-curve-coefnames calc-curve-varnames v c defvars record-entry calc-arg-values calc-dollar-values calc-dollar-used calc-hashes-used with-y coefs vars p calc-fit-to-trail calc-graph-no-auto-view calc-fit-bell-shaped-logistic-curve 111 math-map-binop calcFunc-div calc-fit-hubbert-linear-curve #s(hash-table size 6 test eq rehash-size 1.5 rehash-threshold 0.8125 purecopy t data (101 661 69 661 120 722 88 722 108 779 76 779)) (var mul var-mul) 101 (var exp var-exp) (calcFunc-lambda (var a var-a) (^ 10 (var a var-a))) (var mul var-mul) 120 calcFunc-exp ^ 10 108 (var ln var-ln) (var log10 var-log10) 113 math-add * 103 math-read-expr "(AFit / BFit sqrt(2 pi)) exp(-0.5 * ((XFit - CFit) / BFit)^2)" (vec (var XFit var-XFit)) (vec (var AFit var-AFit) (var BFit var-BFit) (var CFit var-CFit)) (36 39 117 85) 39 calc-do-alg-entry "Model formula: " calc-curve-fit-history "Bad format" reverse calc-invent-args #s(hash-table size 2 test eq rehash-size 1.5 rehash-threshold 0.8125 purecopy t data (117 1083 85 1093)) calc-var-value var-Model1 var-Model2 "User model not yet defined" (3 4) math-objvecp math-all-vars-but "Incorrect model specifier" calcFunc-eq "Not enough variables in model" calc-record "modl" beep calc-enter-result substring symbol-name 9 calc-normalize "parm" "%s" calc-graph-delete calc-graph-add-curve calc-graph-lookup math-contains-sdev-p calc-graph-set-styles calc-graph-point-style intv math-sub math-min-list (float 5 -1) (float 5 -1) math-max-list calc-top-n calc-graph-plot] 12]] 3 nil "P"])
(defalias 'calc-invent-independent-variables #[(n &optional but) "\302 \303\304$\207" [n but calc-invent-variables (x y z t) "x"] 5])
(defalias 'calc-invent-parameter-variables #[(n &optional but) "\302 \303\304$\207" [n but calc-invent-variables (a b c d) "a"] 5])
(defalias 'calc-invent-variables #[(num but names base) "\306\307\306 \307V\203D
\203D \310
:\203
@\202& \311\nT\211!P!\312\313B \"\2049 \fB S
9\204
A\202 \307U\203O \f\237\202W \314\315$,\207" [num var nn n vars names nil 0 math-build-var-name int-to-string math-expr-contains vec calc-invent-variables t base but] 5])
(defalias 'calc-get-fit-variables #[(nv nc &optional defv defc with-y homog) " \203\n \nT\202 \nU\204 \306\307!\210 \250\203 \310\f\203#
S \204+ \311!&\2046 \312
\"&\313\314\315\316\317\320\321 \"\322#\316\317\320\323&\"\322##!\310 ' (\324\325(\"\203` \326(!\202g \326\327(\330Q!\211(\242\306=\203x \306\331\332(8\"\210(\211\242\333=\262\204\207 \306\334!\210(\335\232\203\234 \333 B(\333&B'\202\332 (A@\211\242\333=\262\203\324 (G\336U\203\315 \332(8\211\242\333=\262\203\315 \332(8'(A@(\202\332 \306\337!\210\202\332 \333&B'(GSU\204 \204\361 T(GSU\204\306\340\211\341U\203\376 \342\202\377 \343#\210
'GSU\204\306\344
\211\341U\203\342\202\343#\210\f\203'\333\f'ABB')\2038\345*)A(A#*+\203I\345*+A'A#*()'\211+*\207" [nv with-y calc-curve-nvars defv homog nc error "Wrong number of data vectors for this type of model" nil calc-invent-independent-variables calc-invent-parameter-variables read-string format "Fitting variables (default %s; %s): " mapconcat symbol-name mapcar #[(v) "A@\207" [v] 1] "," #[(v) "A@\207" [v] 1] string-match "\\[" math-read-expr "[" "]" "Bad format in expression: %s" 2 vec "Expected a variable or vector of variables" (vec) 3 "Expected independent variables vector, then parameters vector" "Expected %d independent variable%s" 1 "" "s" "Expected %d parameter variable%s" math-multi-subst defc coefs vars calc-curve-varnames calc-curve-model calc-curve-coefnames] 10])
(defalias 'math-newton-root #[(expr deriv guess orig-guess limit) "\306=\203 \307\310 \"\210 \311\211 \312
!\312! \211:?\206% \211@\313>\262\203\247 \f\211:?\2065 \211@\313>\262\203\247 \f\211:\203Y \211@\314=\203Q \211A@\315=\262\202^ \316!\262\202^ \211\315=\262\204\247 \317 \320 \f\"\"\321 \322 !\211\"\203 \323 \312
!E\202\253 \324\317 \"! \325\"\326U\266\202\203\240 \327
%\202\253 \330 \331\"\202\253 \330 \332\"+\207" [calc-display-working-message guess var-DUMMY next dval expr lots math-do-working "newton" nil math-evaluate-expr (frac float cplx polar) float 0 math-zerop math-sub math-div math-nearly-equal math-float vec math-abs-approx math-compare -1 math-newton-root math-reject-arg "*Newton's method failed to converge" "*Newton's method encountered a singularity" deriv orig-guess limit] 7])
(defalias 'math-newton-search-root #[(expr deriv guess vguess ostep oostep low vlow high vhigh) "\306\307\211\211\211 '\203& (\310=\203\236 \311\312\313\314)*F\"\210\202\236 (\310=\2038 \311\315\313\314)*F\"\210\316\317*)\"\320\"+\321)+\"\211'\322,!\211-\211:?\206^ \211@\323>\262\204\236 \316+\324\"+\321)+\"\211'\322,!\211-\211:?\206\205 \211@\323>\262\204\236 \316+\325\"+\321)+\"\211'\307--\204\251 \322,!--\211:?\206\265 \211@\323>\262\204\277 \326\327\"\210\330-!-.\211:\203\353 \211@\331>\203\343 \211A@\211\250\205\336 \211\314W\266\202\202\360 \332!\262\202\360 \211\314W\262-\211:\203\211@\333>\203\211A@\211\250\205 \211\314V\266\202\202\334!\262\202\211\314V\262\211=\203+*-/\202g/\211:\203Q\211@\331>\203I\211A@\211\250\205D\211\314W\266\202\202V\332!\262\202V\211\314W\262\f=\203e)-.\202g\307-\211:\203\207\211@\335=\203\211A@\314=\262\202\214\336!\262\202\214\211\314=\262\204\230\337)*\"\203\240\340-E\202b\3220!\211\211:?\206\260\211@\323>\262\203\204 \211:\203\324\211@\335=\203\314\211A@\314=\262\202\331\336!\262\202\331\211\314=\262\204\204\341-\330 !\"\317 \"\211\203\204\342*\n\"\204\204\342\n)\"\204\204\n'\322,!\211\211:\203$\211@\335=\203 \211A@\314=\262\202)\336!\262\202)\211\314=\262\2043\337\n\"\203:\340\n E\202b
\203m\342\3431\206K\317*)\"!\343\316\344 \"!\"\203m\345,0\307\211\211+).*/&\n\202b\345,0\n +).*/&\n\202b.\211:\203\252\211@\333>\203\242\211A@\211\250\205\235\211\314V\266\202\202\257\334!\262\202\257\211\314V\262\203\340/\211:\203\330\211@\333>\203\320\211A@\211\250\205\313\211\314V\266\202\202\335\334!\262\202\335\211\314V\262\204<.\211:\203\211@\331>\203\376\211A@\211\250\205\371\211\314W\266\202\202 \332!\262\202 \211\314W\262\203N/\211:\2034\211@\331>\203,\211A@\211\250\205'\211\314W\266\202\2029\332!\262\2029\211\314W\262\203N\346,0).*/&\202b\345,0\307\211\211+).*/&\n.\207" [guess vnext next step pos better t nil lots math-do-working "newton" intv 0 "bisect" math-mul-float math-sub-float (float 5 -1) math-add-float math-evaluate-expr (frac float) (float 6 -1) (float 123456 -5) math-reject-arg "*Newton's method encountered a singularity" math-float (frac float) math-negp (frac float) math-posp float math-zerop math-nearly-equal vec math-div-float math-lessp-float math-abs (float 2 0) math-newton-search-root math-search-root var-DUMMY calc-display-working-message low high ostep expr vguess vlow vhigh deriv oostep] 12])
(defalias 'math-search-root #[(expr deriv low vlow high vhigh) "\306 \203 \307 \310=\203 \n\311\\\202 \312 \310=\203 \313\202 \314\306\211\211\211 - . / 0\3151!\203< \3152!\204J \3161!\203\376 \3162!\203\376 3\317=\203\\ \320\321\322\30745F\"\2100T\2110/V\203s \323\322\30745F\324\"\2100\nU\203} \325.\326\32754\".\"\211\211:\203\247 \211@\330=\203\237 \211A@\307=\262\202\254 \331!\262\202\254 \211\307=\262\203\271 \3325\311\"5\202. \333\3341!\3342!\"\203\342 \3354-\3364 \"\211461\3377!1\202. \3065-\3405 \"\211562\3377!2\202. -\203\f\203-5
2\202-4
1\335.\202D1\211:?\206+\211@\341>\262\2046\3231\342\"\2102\211:?\206B\211@\341>\262\204M\3232\342\"\21045D12D1\211:\203}\211@\343>\203u\211A@\211\250\205p\211\307V\266\202\202\202\315!\262\202\202\211\307V\262\307\32754\"\306\211\211 6 8 9 : ; < = >;T\211;\344X\203B\204B>9=8\326:\345\":9A\203\234\204\2348@\211:?\206\321\211@\341>\262\203\3409@48@1\3469@:\"\21156\3377!23\317=\203\377\320\3475\"\2102\211:?\206 \211@\341>\262\203 \3162!<=\203 \335\202\272959AB\241\210828AB\241\2109AA98AA8\202\272.\203\3242\211:\203h\211@\330=\203`\211A@\307=\262\202m\331!\262\202m\211\307=\262\203y\35052E\202\3361\211:\203\231\211@\330=\203\221\211A@\307=\262\202\236\331!\262\202\236\211\307=\262\203\252\35041E\202\336?\203\305\3517?\306\211\211\2114152&\n\202\336\35274152%\202\336\323\322\35345F\354\")\207" [found math-root-widen calc-internal-prec diff waslow vprev nil 0 point 10 20 (float 9 0) (float 16 -1) math-posp math-negp lots math-do-working "widen" intv math-reject-arg "*Unable to bracket root" (float 16 -1) math-mul-float math-sub-float float math-zerop calcFunc-incr math-lessp-float math-abs t math-sub math-evaluate-expr math-add (frac float) realp (frac float) 5 (float 497 -3) math-add-float "search" vec math-newton-search-root math-bisect-root 3 "*Unable to find a sign change in this interval" prev factor iterlim iters vlow vhigh calc-display-working-message low high var-DUMMY expr yp xp step levels pos yvals xvals deriv] 12])
(defalias 'math-bisect-root #[(expr low vlow high vhigh) "\306 \"\n\211:\203) \211@\307>\203! \211A@\211\250\205 \211\310V\266\202\202. \311!\262\202. \211\310V\262\312\211\211 \313 \314\315\"\316 \"\211\"\204\305 \f\317!\211\211:\203r \211@\320=\203j \211A@\310=\262\202w \321!\262\202w \211\310=\262\204\305 \322=\203\206 \323\324\f\"\210 \211:\203\253 \211@\307>\203\243 \211A@\211\250\205\236 \211\310V\266\202\202\260 \311!\262\202\260 \211\310V\262=\203\275 \f \2028 \f \2028 \325\f E-\207" [high low vhigh vmid mid var-DUMMY math-sub-float (frac float) 0 math-posp nil math-nearly-equal math-mul-float (float 5 -1) math-add-float math-evaluate-expr float math-zerop lots math-do-working "bisect" vec pos step expr calc-display-working-message vlow] 6])
(defvar math-root-vars [(var DUMMY var-DUMMY)])
(defalias 'math-newton-multi #[(expr jacob n guess orig-guess limit) "\306\307\211\211\211
AT\211 W\203* \310 H8
@L\210\202
\311 !\311 !\312 !\203@ \312\n!\204E \313\314\"\210\315\316\317\320 !!\317\n!\"\" \321=\203c \322\323 \"\210
A\fA\211\203v \324
@\f@\"\204c
\203\246 \325\326 !\"!\"\327\"\306U\266\202\203\237 \330 !\"&\202\252 \313\307\331\"\202\252 \332 E.\207" [guess next jacob-val expr-val p2 p -1 nil 2 math-evaluate-expr math-constp math-reject-arg "*Newton's method encountered a singularity" math-add math-div math-float math-neg lots math-do-working "newton" math-nearly-equal math-abs-approx math-sub math-compare math-newton-multi "*Newton's method failed to converge" vec m n math-root-vars expr jacob calc-display-working-message orig-guess limit] 7])
(defalias 'math-find-root #[(expr var guess math-root-widen) "\242\306=\203rGS\307\211\306C\307\211\211\211 @ A B\242\306=\204+ \310\311\"\210GBTU\2048 \312 \210\313!BCGY\203a \314\315\316CG!P! D\317C\320\321D\211E!\"C)\202< \322\nT\211BW\203y \323C\nH8\307L\210\202c \322\nT\fA\211\203\257 \f@\242\321=\204\226 \310\324\"\210 A\211\203~ \325 @\f@C\nH#\240\210\202\230 E\242\306=\204\275 \310E\311\"\210EGBTU\204\312 \312 \210\313E!\211E\fA\211\203\373 E@\211:?\206\346 \211@\326>\262\204\361 \310E\327\"\210\f\330\f@!\240\210\202\322 \fA\211\203R\f@\242F\236\203\f\331\f@A@\323\f@8\"\240\210\f\332\f@!\240\210\306C\322\nT\211BW\203J \332\333\f@C\nH\307\334$\206C\310\335\"!C\244\210\202&
C\244\210\202\375 \336E!\337
BE\211\340\n!\203i\341\202m\342\n\343\"&.\207\242\321=\204\200\310\324\"\210\344\"\204\215\310\345\"\210@F\236\203\236\331A@\3238\"\346G\323\\ G\325\347#\307\211 A @\332!\333\350\307\334$ H\307\211 I J\307\211 K LH\203\325\332H!H\330E!\211E\211:?\206\346\211@\326>\262\203
H\203
\351HE\211\340E!\203\352\202 \342\336E!\353\"%\202>E\211:?\206\211@\354>\262\2038EIEJE@\332!\211KL\355M\202\246E@\356=\203\222\357E!\204M\310E\360\"\210\323E8I\361E8JEA@\362>\203k\363I\364J#IEA@\365>\203}\363J\322I#JI@\332!KJ@\332!L\202\246\366E!\203\240\310\367!\210\202\246\310E\370\"\210K\211:\203\306\211@\371=\203\276\211A@\372=\262\202\313\340!\262\202\313\211\372=\262\203\327\306IKE\202>L\211:\203\367\211@\371=\203\357\211A@\372=\262\202\374\340!\262\202\374\211\372=\262\203\306JLE\202>H\203DK\211:?\206\211@\326>\262\203DL\211:?\206*\211@\326>\262\203D\373H\307\211\211\211IKJL&\n\202>K\211:\203j\211@\374>\203b\211A@\211\250\205]\211\372V\266\202\202o\375!\262\202o\211\372V\262\203\240L\211:\203\230\211@\374>\203\220\211A@\211\250\205\213\211\372V\266\202\202\235\375!\262\202\235\211\372V\262\204 K\211:\203\306\211@\376>\203\276\211A@\211\250\205\271\211\372W\266\202\202\313\377!\262\202\313\211\372W\262\203\374L\211:\203\364\211@\376>\203\354\211A@\211\250\205\347\211\372W\266\202\202\371\377!\262\202\371\211\372W\262\204 K\211:?\206\211@\326>\262\203 L\211:?\206\211@\326>\262\2041\201N HIKJL&\202>\201O IKJL%. !\207" [expr row m p2 p jacob vec nil math-reject-arg vectorp math-dimension-error copy-sequence intern "math-root-v" int-to-string vconcat vector var -1 2 "*Expected a variable" math-expr-subst (frac float cplx polar) numberp math-float math-sub math-evaluate-expr calcFunc-deriv t "*Formulas must be differentiable" math-abs-approx math-newton-multi math-zerop (float 1 3) math-mul 10 math-expr-contains "*Formula does not contain specified variable" math-normalize (var DUMMY var-DUMMY) (var DUMMY var-DUMMY) math-newton-root (float 1 6) 100 (frac float) point intv math-constp constp 3 (0 1) calcFunc-incr 1 (0 2) math-complexp "*Complex root finder must have derivative" realp float 0 math-newton-search-root (frac float) math-posp (frac float) math-negp var-DUMMY calc-symbolic-mode n math-root-vars symb guess calc-tweak-eqn-table calc-internal-prec deriv low high vlow vhigh math-root-widen math-search-root math-bisect-root] 13])
(defalias 'calcFunc-root #[(expr var guess) "\303 \n\304$\207" [expr var guess math-find-root nil] 5])
(defalias 'calcFunc-wroot #[(expr var guess) "\303 \n\304$\207" [expr var guess math-find-root t] 5])
(defvar math-min-vars [(var DUMMY var-DUMMY)])
(defalias 'math-min-eval #[(expr a) "\211\242\305=\262\203& \306 TA\211\203\" \307\n H8@L\210\202\f )\202( \310\f!\211\211:?\2067 \211@\311=\262\203@ \312!\207@\313=\203I \207\314\315\"\207" [a m math-min-vars var-DUMMY expr vec -1 2 math-evaluate-expr frac math-float float math-reject-arg realp] 4])
(defvar math-min-or-max "minimum")
(defalias 'math-widen-min #[(expr a b) "\306\307\306\211\211\211\211\211\211\211\211\211\211\211\211 $ % & ' ( ) * + ,-\2043 \310.\311\"-\312/.\"(\312/-\"'\313('\"\203d .%-.%-($'($'\314-\315\316\317-.\"\"\")\312/)\"&,\204\313&'\"\2030\320=\203\233 \321\322\323\324.)F\"\210+S\211+\324U\203\257 \325\306\326\3271\"\"\210\317-)\"\317-.\"\315 \317'&\"\"\315\n\317'(\"\"\317\f
\"\313\330!\331\"\203\353 \332!\203\351 \333\202\352 \334\317-\335\317\315\n\f\"\315
\"\"\315\336\"\"\"%\314-\315\337\n\"\"\332\n!\211*\203 \313-%\"\202\"\313%-\"\203\353*\2033\313%)\"\2029\313)%\"\203\216\313\312/%\"\211$&\"\203b-.'(%-$'\340,\202\375\313'$\"\203y%)$&\340,\202\375\314)\315\341\n\"\"%\312/%\"$\202\375*\203\233\313% \"\202\240\313 %\"\203\335\313\312/%\"\211$&\"\203\375)-&'%)$&\314)\315\342\317-)\"\"\"%\312/%\"$\202\375 %\312/%\"$\202\375\314)\315\343\n\"\"%\312/%\"$-.'()-&'%)$&\202{ \313.)\"\2032.(-')&\257\202@)&-'.(\257.\207" [qr ba bc ulim q r nil 30 math-mul (float 101 -2) math-min-eval math-lessp-float math-add-float math-mul-float (float 161803 -5) math-sub-float lots math-do-working "widen" intv 0 math-reject-arg format "*Unable to find a %s near the interval" math-abs (float 1 -20) math-negp (float -1 -20) (float 1 -20) math-div-float (float 2 0) (float -1 2) t (float -161803 -5) (float -161803 -5) (float -161803 -5) vu u vc vb va c incr iters done b a expr calc-display-working-message math-min-or-max] 16])
(defalias 'math-narrow-min #[(expr a c intv) " D\306\n\"\306\n \"D\307\310 \"\311\211\211\211 T\211 \312X\203\270 \204\270 \313\314\"
A\203\" \204\" \315
@\" \316=\203\\ \317\320 \"\210
AB\241\210\f\306\n \"\fAB\241\210\321\fA@\f@\"\203\210 \321\fA@\322\f8\"\203\210 \323\202?
A\fA\211AA\203\257 \321\fA@\f@\"\203\257 \321\fA@\322\f8\"\203\257 \323\202?
A\fA\202? \203\322
@\f@
A@\fA@\322
8\322\f8\257\202S\321 @ A@\"\203 A@\324>\203K @ ! A\211 \203\377 \321! @\"\204\354 ?\205\f\322 8!D)\202H A@\325>\203K \237\211 \203K @ ! A\211 \203:\321! @\"\204' ?\205G\326 8!D)\206S\327\311\330\331\"\"\".\207" [a c expr b yp xp math-min-eval 0 math-sub-float nil 5 math-mul-float (float 497 -3) math-add-float lots math-do-working "search" math-lessp-float 2 t (2 3) (1 3) 3 math-reject-arg format "*Unable to find a %s in the interval" found step levels yvals xvals calc-display-working-message intv min math-min-or-max] 9])
(defalias 'math-brent-min #[(expr prec a va x vx b vb) "\306\307_\\ \n \n\310\311\312ZE\310\311\313ZE\314\315\211\211\211\211\211\211\211\211\211\211\211 ' ( ) * + , - . / 0 1 2 3 4 5 6 7\316\317\32089\"\",\3201\3162\321 !\"\"+\316+\322\"*\323\324*\316\325\32498\"\"\"\321\324 ,\"!\"\203f7S\2117\326U\203\223 \327\315\330\331:\"\"\210;\332=\203\237 \333\334 \"\210\323\3210!+\"\203\313 \323 ,\"\203\272 \3249 \"\202\277 \3248 \"0\316\3350\"/\202\260\324 6\"\316 \324\n3\"\"\324 4\"\316\f\324\n5\"\"'\324\316\f'\"\316
\"\"(\316\336\324'
\"\"'\337'!\203\340(!(\202\340'!'0)/0\323\321(!\321\316\341\316')\"\"!\"\203\222\323\316'\3248 \"\"(\"\203\222\323(\316'\3249 \"\"\"\203\222\342('\"/\320 /\".\323\324.8\"*\"\204|\323\3249.\"*\"\203\260\323, \"\203\213\340+!\202\215+/\202\260\323 ,\"\203\242\3249 \"\202\247\3248 \"0\316\3430\"/\320 \323\321/!+\"\203\320\344/!\203\313\340+!\202\322+\202\322/\".\345<.\"-\323\n-\"\203>\323. \"\203\364.8\202\370.96 \232\204\3235-\"\2046453.6-5\202G 4 \232\204346\232\2043\3233-\"\204G .4-3\202G \323. \"\203L 9\202O 86453 6\n5.-\202G \346 \nE.\207" [prec x vx xw xv r 20 5 float 1 -1 -5 (float 0 0) nil math-mul-float (float 5 -1) math-add-float math-abs (float 2 0) math-lessp-float math-sub-float (float 5 -1) 0 math-reject-arg format "*Unable to converge on a %s" lots math-do-working "brent" (float 381966 -6) (float 2 0) math-posp math-neg-float (float 5 -1) math-div-float (float 381966 -6) math-negp math-min-eval vec q p etemp tol2 tol1 xm vu u d e zeps tol vv v vw w iters a b math-min-or-max calc-display-working-message expr] 20])
(defalias 'math-powell-min #[(expr n guesses prec) "\306 \"\307\310 \"\311\312\313\f\"B\211 !\314\310\"[E #\315
\" $\316\211 % &\316\211 ' (\316\211 ) *\316\211 + ,\316 -$%\317*\320+\317))T\211) X\203\226 $'\321\n
\322 )\" \"%-\" .\323-@!)\324-8$\325\326'$\"!,\327+,\"\203J ,+)*\202J \327\330#\331\325%!\325$!\"\"\330\332\325\326%$\"!\"\"\203<\333\334\335
\"!\"&\333
!\"(
!\315&\"'\327'%\"\203= \327\330\211\336\331\326%\330\337$\"\"'\"\"\340\326\211%$\"+\"!\"\330+\340\326%'\"!\"\"\203= \321\n
( \"%\211-@\324-8$\317))T\211) X\203= *) 8\233)-A@8\240\210\202 \311
$E.\207" [expr n f1dim xi guesses p math-line-min-func calcFunc-idn 1 vec mapcar car float math-min-eval nil 0 (float 0 0) math-line-min math-mat-col math-normalize 2 math-abs math-sub-float math-lessp-float math-mul-float math-add-float (float 2 0) math-sub math-mul (float 2 0) (float 2 0) (float 2 0) math-sqr-float pt prec ftol fret fp ptt fptt xit i ibig del diff res calc-internal-prec] 11])
(defalias 'math-line-min-func #[(expr n) "\304T\211 W\203$ \305\nH8\306\307\310\311\312TEE\311\313TEEL\210\202 \314 !)\207" [m n math-min-vars expr -1 2 + * (var DUMMY var-DUMMY) calcFunc-mrow (var line-xi line-xi) (var line-p line-p) math-evaluate-expr] 7])
(defalias 'math-line-min #[(f1dim line-p line-xi n prec) "\306\307 !\310\n\311\312#\313\314\n\f $ \315
A@\" \316\"\317
8E-\207" [var-DUMMY f1dim expr params prec res nil math-evaluate-expr math-widen-min (float 0 0) (float 1 0) apply math-brent-min math-mul math-add 2 line-xi xi line-p] 5])
(defalias 'math-find-minimum #[(expr var guess min-widen) "\306\307\306 \211\242\310=\262 \306\211 - \211\242\310=\262\204! \310 D.\211\242\310=\262\2042 \310.D. G.GU\204> \311 \210 A.A\211.\203O @\242\303=\204X \312 @\313\"\210\314/ @\"\204g \312 @\315\"\210 T0GY\203\214 \316\317\3200G!P! 1\3210\322\3031\211E!\"0)\202g \3230 H8\306L\210\3230 TH8\306L\210\324.@!\203\340 \325/ @\3260 H\3270 TH\330EE#/\331\332.@!! \323
8\306\211E
A@\306\211\333F-)BB- \323\\\202> \325/ @0 H#/.@\211:?\206\371 \211@\334>\262\203 \331.@!\306\211E\202D.@\242\335=\203>\336.@!\203>\337\340\323.@8\341.@8\"\342\"\331\323.@8!\331\341.@8!.@F\202D\312.@\343\"-B- T\202> -\237-\344/!/ \345U\203\340-@A@\203\2242\203~\346/-@A@\323-@8#\202\235\347/-@A@\323-@8\341-@8$\202\235\346/-@@\306# 34 53AA\203\301\3504\211\323\\\\ 4\351\352/53$)!\202\305\3103B 6\f\203\332\310\2116A@D\32368E\202\3346+\202[4 5\3504\211\323\\\\ 4\353/ -5$)!\211 6A@ 7\310C 7A\2117\203A\341-@8\2031\350\3547@7A@E!C\244\2107A7-A-\20297@C\244\210-A-\202\f\203P\310\32368E\202Z\310A@\32368E,.\207" [calc-symbolic-mode n var-DUMMY var isvec g nil 0 vec math-dimension-error math-reject-arg "*Expected a variable" math-expr-contains "*Formula does not contain specified variable" intern "math-min-v" int-to-string vconcat vector 2 math-complexp math-expr-subst + * (cplx 0 1) math-float math-complex t (frac float) intv math-constp math-mul math-add 3 (float 5 -1) realp math-evaluate-expr 1 math-widen-min math-narrow-min math-normalize apply math-brent-min math-powell-min cplx guesses guess expr math-min-vars symb min-widen params calc-internal-prec prec res p] 9])
(defalias 'calcFunc-minimize #[(expr var guess) "\305\245\306]\307\310\311\n!\311 !\311\f!\312$*\207" [calc-internal-prec math-min-or-max expr var guess 2 3 "minimum" math-find-minimum math-normalize nil] 5])
(defalias 'calcFunc-wminimize #[(expr var guess) "\305\245\306]\307\310\311\n!\311 !\311\f!\312$*\207" [calc-internal-prec math-min-or-max expr var guess 2 3 "minimum" math-find-minimum math-normalize t] 5])
(defalias 'calcFunc-maximize #[(expr var guess) "\306\245\307]\310\311\312\313\n!!\312 !\312\f!\314$ \315
A@\313\306
8!E+\207" [calc-internal-prec math-min-or-max expr var guess res 2 3 "maximum" math-find-minimum math-normalize math-neg nil vec] 5])
(defalias 'calcFunc-wmaximize #[(expr var guess) "\306\245\307]\310\311\312\313\n!!\312 !\312\f!\314$ \315
A@\313\306
8!E+\207" [calc-internal-prec math-min-or-max expr var guess res 2 3 "maximum" math-find-minimum math-normalize math-neg t vec] 5])
(defalias 'calcFunc-polint #[(data x) "\304!\204 \305\306\"\210G\307U\204 \305\310\"\210A@G\311V\204% \305\312\"\210 \211\242\313=\262\203A \314 !\2049 \n\203A \313\315\316 A\"B\207\317 !\204P \n\204P \305 \320\"\210\321 \311\\\313\322A@A\3118A \323$)B!\207" [data x math-expand-formulas calc-internal-prec math-matrixp math-reject-arg matrixp 3 "*Wrong number of data rows" 2 "*Too few data points" vec math-constp mapcar #[(x) "\302 \"\207" [data x calcFunc-polint] 3] math-objectp objectp math-normalize math-poly-interp nil] 7])
(put 'calcFunc-polint 'math-expandable t)
(defalias 'calcFunc-ratint #[(data x) "\304!\204 \305\306\"\210G\307U\204 \305\310\"\210A@G\311V\204% \305\312\"\210 \211\242\313=\262\203A \314 !\2049 \n\203A \313\315\316 A\"B\207\317 !\204P \n\204P \305 \320\"\210\321 \311\\\313\322A@A\3118A A@AAA$)B!\207" [data x math-expand-formulas calc-internal-prec math-matrixp math-reject-arg matrixp 3 "*Wrong number of data rows" 2 "*Too few data points" vec math-constp mapcar #[(x) "\302 \"\207" [data x calcFunc-ratint] 3] math-objectp objectp math-normalize math-poly-interp] 7])
(put 'calcFunc-ratint 'math-expandable t)
(defalias 'math-poly-interp #[(xa ya x ratp) "G\306\211\211\307 !\307 !\310\211\306\211\306\211\211\211 T\211 X\203h \311@ \"BA\312@! \203^ \313\n \"\203. \n \202. \237 S\211 8\314 !\203\205 \310D\202aT\211 W\203\\\310\233T\211 ZX\203/ \203\353 \315\316@ @\"
@\" !\315\311\fA@ @\"\311!\fA@\"\" \316\fA@\n\"\240\210\f\316!\n\"\240\210)\202 \317@
@\"\203\374 \320\321B\322\"\210\315\311\fA@ @\"\311@
@\"\" \316
@\n\"\240\210\f\316@\n\"\240\210\fA AA
A\202\244 \211\\ ZW\203F 8\202Q S\211 8\323\"\202\205 D.\207" [xa ya temp dp cp xpm nil copy-sequence 0 math-sub math-abs math-lessp math-zerop math-div math-mul math-equal math-reject-arg vec "*Duplicate X values" math-add xp dy y m i d c xax ns dif n x ratp t2] 16])
(defalias 'calcFunc-ninteg #[(expr var lo hi) "\306!\306 !\211:?\206 \211@\307>\262\204# \310!\204# \311\312\"\210 \211:?\206. \211@\307>\262\204> \310 !\204> \311 \312\"\210\313 \"\203N \314\315\n $!\207\316\n \317#\320\211\211\321 \306\n!\322\232\203\207 \313 \323\"\203s \202t \324 \325\326\n\327!\327 !\330% ) \331\232\203\257 \313\332\"\203\230 \202\231 \333 !\334\f\325\326\n\327!!\327 !\330%\"!)\335 \"\204\305 \334\f\325\326\n\327!\327 !\320%\"\f,\207" [lo hi expr var sum calc-prefer-frac math-evaluate-expr (frac float cplx polar) math-infinitep math-reject-arg numberp math-lessp math-neg calcFunc-ninteg math-expr-subst (var DUMMY var-DUMMY) nil 0 (neg (var inf var-inf)) (float -2 0) (float -2 0) math-ninteg-romberg math-ninteg-midpoint math-float inf (var inf var-inf) (float 2 0) (float 2 0) math-add math-equal calc-symbolic-mode var-DUMMY thi tlo] 8])
(defalias 'math-ninteg-romberg #[(func expr lo hi mode) "\306\307\211\310\307\307 \311\312\\ \204y \fT\211\313X\203y
$C\244C\244\f\314Y\203a \315
\316\307$ \317\320 A@!\321\320 @!\n[\"\"\203` @)\f\322Y\203o
AA\323\324\"\202 \206\203 \325\307\326\327!\")!.\207" [calc-internal-prec math-ninteg-temp prec ss j s (float 1 0) nil 0 math-normalize 2 8 3 math-poly-interp (float 0 0) math-lessp math-abs calcFunc-scf 5 math-div-float (float 9 0) math-reject-arg format "*Integral failed to converge" h curh func expr lo hi mode res] 7])
(defalias 'math-ninteg-evaluate #[(expr x mode) "\305=\203 \306\307 \" \310 !\211 \211:?\206 \211@\311>\262\204& \312\f\313\"\210\305=\2033 \314\f\315 !\"\f*\207" [mode x var-DUMMY expr res inf math-div (float 1 0) math-evaluate-expr (frac float cplx polar) math-reject-arg "*Integrand does not evaluate to a number" math-mul math-sqr] 5])
(defalias 'math-ninteg-midpoint #[(expr lo hi mode) "\306=\203 \307\310\311\312 \"\311\312\f\" *
\203\257 \313
@_ \314
@_ \315 \316\f \" \311\317!\" \320\211\" \320 \" \320 \321\322\"\" \323 \315 \310 T\211
@X\203\230 T\324! #T\320 \320 \324!\320 \"#\"\" \320 \" \202[ \320\311
A@\325\"\321 \"\"D. \202\304 \312\321\316\f \"\324!\321\320 \f\"\326\"#\"D
A@\207" [mode temp math-infinite-mode lo hi math-ninteg-temp inf t nil math-div 1 3 2 0 math-sub math-float math-add math-mul (float 5 -1) (float 0 0) math-ninteg-evaluate (float 3 0) (float 5 -1) it3 math-working-step-2 math-working-step range del del2 del3 x sum j expr] 9])
(defvar math-dummy-vars [(var DUMMY var-DUMMY)])
(defvar math-dummy-counter 0)
(defalias 'math-dummy-variable #[nil " GU\203 \303\304\305\"!\306 \307\310\n\211E!\")\311 H8\312L\210 HT\207" [math-dummy-counter math-dummy-vars symb intern format "math-dummy-%d" vconcat vector var 2 nil] 6])
(defvar math-in-fit 0)
(defvar calc-fit-to-trail nil)
(defalias 'calcFunc-fit #[(expr vars &optional coefs data) "\306\307 \310\\\311\n \f
\312%)!)\207" [math-in-fit calc-internal-prec expr vars coefs data 10 math-normalize 2 math-general-fit nil] 7])
(defalias 'calcFunc-efit #[(expr vars &optional coefs data) "\306\307 \310\\\311\n \f
\312%)!)\207" [math-in-fit calc-internal-prec expr vars coefs data 10 math-normalize 2 math-general-fit sdev] 7])
(defalias 'calcFunc-xfit #[(expr vars &optional coefs data) "\306\307 \310\\\311\n \f
\312%)!)\207" [math-in-fit calc-internal-prec expr vars coefs data 10 math-normalize 2 math-general-fit full] 7])
(defalias 'math-general-fit #[(expr vars coefs data mode) "\306\307 \310=\306\n\306\211\211\211\211\211\211\211\211\211\211\211\211\211\211\211\211\211\211\211\211\211 > ? @ A B C D E F G H I J K L M N O P Q R S T U VW\204\237 X\203i XW\306X\202\237 \n\211\242\311=\262\203\225 \nG\312>\203\217 YW\313\n8Y\314\n8X\nA@\202\237 \315 \210\202\237 YW\306\211YX\316W!\204\254 \317W\320\"\210WGSWA@GSAY\211\242\311=\262\204\317 Y\203\317 \311YDYX\211\242\311=\262\204\345 X\203\345 \311XDXX\204\363 \311\321\nY\"BXY\2041XGS
X\203\n\317X\322\"\210\2021\323X!X\324XG
\n\242\325=\203 \307\202!\326#X\233\311 ABY \306\241\210)YGS
U\204HYG
U\204H\317Y\327\"\210XGS\211?\307W\203Z\317X\330\"\210X A\211\203\204 @\242\331=\204r\317 @\332\"\210\333 \334\n @\335SZD#\202]RY A\211\203\261 @\242\331=\204\237\317 @\332\"\210\333 \334\n @\336RZD#\202\212R
\\W\203\275\333 \fJ\nP\n\242\325=\204\321\325\336
D\nE\306 Z\337\n!\340\341\nD\342\"\313U\337\n!\211\242\343=\203 \nG\344U\203 \313\n8\211\242\311=\262\203 \314\n8\211\242\311=\262\203 \313\n8G\307V\203 \314\n8G?TU\204&\317Q\345\"\210\nA@K\346\313\n8A!L\314\n8ILG>KJ\232\203H\306K\311\347>\326\"BB\311\350\323\347>B\"\"BC\346WA! [\307 \\\351>\326\" ]\326 ^\306\211 _ `\306\211 a b\306\211 c d\306\211 e f\306\211 g h\306 i^T\211^AX\203I\326_\306c_
W\203.[_[_HAI\210[_H@a_
SU\203\371c\203a\242\352=\203\347\353\354\313a8!c\"caA@aK\203\355a\356c!\"a\202a\242\352=\203\353\354\313a8!c\206\326\"caA@a\313jR_\\H8aL\210_T\211_\202\247K\203<\337K!b\202C\313J8Jbb\242\352=\203Y\354\313b8!cbA@b^\307U\203tcOT\206r \352=\205rO?NO\203\235c\203\224\357\307c\"\\N\203\250\\FBF\202\250\317b\360\"\210\202\250c\203\250\317b\360\"\210a\242\352=\203\272\313J8aA@L\210\326_CeBh_>W\2035\337L_H!d]_dI\210\361d\\\"dhAheA\211e@f\326\211`_X\203 fA\211f\353f@\361d]`H\"\"\240\210`T\211`\202\370h\353h@\361db\"\"\240\210_T\211_\202\307N\203\223\362]bC\"GBG\202\223\326_Ce_>SW\203\227_`_T\211_eA\211e@\233f_C\233g`T\211`>W\203PfAfgAgf_g@8\240\210\202t. \203\343\363C!\211C\203\262\361CB\"B\202\315\364\365B!!\203\307\366\326BGS\"C\202\315\317P\367\"\210C\211\242\311=\262\204\353\311\211CDDC\202\353\357BC\"BN\203\\\306\211\211 k l m\326MG\203[BmG@l\326kmA\211m\203,\353k\361m@l@\"\"klAl\202
\354\370l@k\"!kF\203E\361kF@\"k\353Mk\"MFAFGA\211G\204+\323B!\211H mC n\307 o\314 U\326 _ \203\222O\204\222A>X\203\207\326\202\220\357MA>Z\"o \203\277mA\211m\203\277m\355m@\356\361_T\211_nA\211n@8o\"!\"\240\210\202\226\337I!Hp\203\373HX\306 n mmAmnA\211n\203\363\325n@m@E B \202\325\311 \237Bp+.\371P
\306 \211 A\307Y\203 \336AD B AS\211A\202?\211A\307Y\203:\335AD B AS\211A\202# *\362HAYA\"#\372T\203\244\311\nBCI\326 AATA A\211\203w @\242\373=\203w @A@A=\204Y \203\200I\202\202\311C*MO\203\236A>V\203\236\374MA>ZE\202\237\375\257\202\245\n!. \207" [math-dummy-counter mode expr p dummy v nil 1 full vec (3 4) 2 3 math-dimension-error math-matrixp math-reject-arg matrixp math-all-vars-but "*Not enough variables in model" copy-sequence - calcFunc-eq 0 "*Number of variables does not match data" "*Need at least one parameter" var "*Expected a variable" math-dummy-variable math-expr-subst calcFunc-fitparam calcFunc-fitvar math-evaluate-expr math-rewrite calcFunc-fitmodel (var FitRules var-FitRules) calcFunc-fitsystem 4 "*Model expression is too complex" vconcat make-list mapcar make-vector sdev math-add math-sqr math-make-sdev math-sqrt math-div "*Mixed error forms and plain numbers" math-mul append math-matrix-inv-raw math-zerop math-abs calcFunc-diag "*Singular matrix" math-sub math-multi-subst math-normalize calcFunc-fitdummy calcFunc-utpc (var nan var-nan) mm m nn n beta covar var-YVALX var-YVAL weights xy-values math-fit-new-coefs coef-filters y-dummy y-filter x-funcs chisq need-chisq have-sdevs orig-expr plain-expr math-fit-first-var math-fit-first-coef extended math-in-fit calc-simplify-mode data coefs vars calc-symbolic-mode ptrs isigsq xvals i j k xval yval sigmasqr wt covj covjk covk betaj lud math-dummy-vars sum xp bp cp sigdat calc-fit-to-trail] 30])
(defalias 'calcFunc-fitvar #[(x) "\304Y\203 \305\n \306#H\307\304 8!\206 \207\310 !\207" [math-in-fit math-dummy-vars math-fit-first-var x 2 + -1 calc-var-value math-reject-arg] 5])
(defalias 'calcFunc-fitparam #[(x) "\304Y\203 \305\n \306#H\307\304 8!\206 \207\310 !\207" [math-in-fit math-dummy-vars math-fit-first-coef x 2 + -1 calc-var-value math-reject-arg] 5])
(defalias 'calcFunc-fitdummy #[(x) "\303U\203\n \n8\207\304 !\207" [math-in-fit x math-fit-new-coefs 3 math-reject-arg] 2])
(defalias 'calcFunc-hasfitvars #[(expr) "\211:?\206 \211@\301>\262\203 \302\207@\303=\203 A@\207\304\305\306\307A\"\"\207" [expr (frac float cplx polar hms date mod var) 0 calcFunc-fitvar apply max mapcar calcFunc-hasfitvars] 5])
(defalias 'calcFunc-hasfitparams #[(expr) "\211:?\206 \211@\301>\262\203 \302\207@\303=\203 A@\207\304\305\306\307A\"\"\207" [expr (frac float cplx polar hms date mod var) 0 calcFunc-fitparam apply max mapcar calcFunc-hasfitparams] 5])
(defalias 'math-all-vars-but #[(expr but) "\304!\n\211\203 \305\306 \242 \" \" A\211\204\n \307\310\311 \"\312\"*\207" [expr vars but p math-all-vars-in delq assoc sort mapcar car #[(x y) "A@ A@\231\207" [x y] 2]] 5])
(defalias 'math-all-vars-in #[(expr) "\303\211\304\n!\210 *\207" [math-all-vars-found math-all-vars-vars expr nil math-all-vars-rec] 2])
(defalias 'math-all-vars-rec #[(expr) "\211:?\206 \211@\303>\262\2034 \242\304=\205C \305!\2063 \306 \"\211\203, \n\211AT\241\207\307B B\211\207A\211\205C \310@!\210\2024 \207" [expr math-all-vars-vars math-all-vars-found (frac float cplx polar hms date mod var) var math-const-var assoc 1 math-all-vars-rec] 3])
(provide 'calcalg3)