| Current Path : /usr/share/emacs/24.3/lisp/calc/ |
Linux gator3171.hostgator.com 4.19.286-203.ELK.el7.x86_64 #1 SMP Wed Jun 14 04:33:55 CDT 2023 x86_64 |
| Current File : //usr/share/emacs/24.3/lisp/calc/calcalg3.elc |
;ELC
;;; Compiled by mockbuild@buildfarm06-new.corp.cloudlinux.com on Fri Oct 11 10:10:48 2024
;;; from file /builddir/build/BUILD/emacs-24.3/lisp/calc/calcalg3.el
;;; in Emacs version 24.3.1
;;; 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)
#@58 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 (#$ . 628)])
(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 (#$ . 2635))
(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\204\207\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] >\203\322\351
\211T\n\205\236\333#\210\352HA@\363\201^ \365D\201_ =\203\271\201` \202\274\201a \365\201b \342HAABI#\"\"\"\211F\2022 D\201c >\203\351
\211T\n\205\344\314#\210\352H\342\333IABB\"FD\201d =\203\201e FD\202\201f \201g FE\211F\2022 D\201h >\203M\351
\211T\n\205&\314#\210\352H\342\333\365D\201i =\203=\201j \202@\201k I\"ABB\"\211F\2022 D\201l U\203\243\351
\211\337_T\n\205a\314#\210HI J\211 KA@FJAJKAA\211K\203\237\201m F\201n K@\201f \355J@KA@E\337EE\"F\202p*\2020 D\201o U\203\317\201p \201q !F\201r I\201s H\351\333HGS\n\205\312\333#\210\2020 D\201t >\203\306\211 L MD\201u =\203DN O\314\211 P Q\201v \323\201w \306\201x $\211FG\333U\204\f\345\201y !\210F@F\373MP\314V\203/\342NGPZ\201z N!\233BH\202@Q\314V\203@\342\201{ Q!BH+\202\334D\201| =\203W\201} \201~ !\202mD\201 =\203j\201} \201\200 !\202m\332\333!\211F\204y\345\201\201 !\210\201\202 F!\203\334FG\201\203 >\203\326\201\204 FA@!\204\326\201\202 \337F8!\203\326\334F8\203\264\201\202 \334F8!\203\326\337F8I\334F8\206\313\342\201\205 FI\"BHFA@F\202\334\345\201\206 !\210I\204fF\242\201\207 = RH\203\351R\203\372
T\202\373
HGS\201\205 FH\"\306R%\210\202e\201\205 F\306\" S\306 T\355SG
R\203)\337\202*\333# @\306 U@\314W\203=\345\201\210 !\210@S\233\211UATU\306\241\210\351R\203X
T\202Y
SGTSR%\210,)M\203{\201\211 \342FIHF\201\212 \"\210*\2020 \201\213 \210\2020 B\204\337\373 V\201\214 @\201\215 E!\201\216 \306OEFIG\337U\203\262IA@\202\264IHG\337U\203\303HA@\202\305H \257#\210V:\203\336\201\211 \201\217 V!\201\220 \"\210)A\205sA;\203\364\316\201\221 A\"\202s\373 W\201\222 \373!\210\201\223 \201\224 AA@!\201\224 \337A8!\"\210\201\225 \337 8!\204)\201\226 \306\211\"\210\201\227 \306!\210AA@AA\201\230 \334\201\231 \201\232 A@AA\"\201\233 \"\201m \201\234 \201\235 A@AA\"\"FA\201\223 \201\224 A!\201\224 \201\236 \333!!\"\210\201\237 \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 (101 69) (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 88) 120 calcFunc-exp ^ 10 (108 76) 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 117 calc-var-value var-Model1 85 var-Model2 "User model not yet defined" math-vectorp (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 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\333)!\204\203 \306\334!\210)\335\232\203\230 \336 B)\336'B(\202\316 \333)A@!\203\310 )G\337U\203\301 \333\332)8!\203\301 \332)8()A@)\202\316 \306\340!\210\202\316 \336'B()GSU\204\365 \204\345 T)GSU\204\365 \306\341\211\342U\203\362 \343\202\363 \344#\210
(GSU\204\306\345
\211\342U\203 \343\202\f\344#\210\f\203\336\f(ABB(*\203,\346+*A)A#+,\203=\346+,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 math-vectorp "Expected a variable or vector of variables" (vec) 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( @\313>)\203\272 \f\211 :?\206: @\313>)\203\272 \f\211 :\203g @\314>?\205k @\315=\203` A@\316=\202k \317 !\202k \316=)\204\272 \320 \321 \f\"\"\322 \323 !\211\"\203\215 \324 \312
!E\202\276 \325\320 \"! \326 \"\327U*\203\263 \330
%\202\276 \331 \332\"\202\276 \331 \333\"+\207" [calc-display-working-message guess var-DUMMY next dval expr lots math-do-working "newton" nil math-evaluate-expr (bigpos bigneg frac float cplx polar) (bigpos bigneg) 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 a orig-guess limit b] 7])
(defalias 'math-newton-search-root #[(expr deriv guess vguess ostep oostep low vlow high vhigh) "\306\307\211\211\211 +\203& ,\310=\203\227 \311\312\313\314-.F\"\210\202\227 ,\310=\2038 \311\315\313\314-.F\"\210\316\317.-\"\320\"/\321-/\"\211+\3220!\2111\211 2:?\206a 2@\323>)\204\227 \316/\324\"/\321-/\"\211+\3220!1\3251!\204\227 \316/\326\"/\321-/\"\211+\30711\204\242 \3220!11\211 2:?\206\261 2@\323>)\204\272 \327\330\"\210\3311!13\211 2:\2032@\332=\2062@\333=?\2052@\334>\203\375 2A@\211 2:\203\365 2@\332=\202\371 2\314W)\202\3352!\2022\314W)1\211 2:\203M2@\333=\206Q2@\332=?\205Q2@\336>\203F2A@\211 2:\203>2@\333=\202B2\314V)\202Q\3372!\202Q2\314V)\211=\203b.14\202\2744\211 2:\203\2462@\332=\206\2522@\333=?\205\2522@\334>\203\2372A@\211 2:\203\2272@\332=\202\2332\314W)\202\252\3352!\202\2522\314W)\f=\203\272-13\202\274\3071\211 2:\203\3462@\340>?\205\3522@\341=\203\3372A@\314=\202\352\3422!\202\3522\314=)\204\367\343-.\"\203\377\3441E\202O\3225!\211\211 2:?\2062@\323>)\203\371 \211 2:\203?2@\340>?\205C2@\341=\20382A@\314=\202C\3422!\202C2\314=)\204\371\3451\331 !\"\317 \"\211\203\371\346.\n\"\204\371\346\n-\"\204\371\n+\3220!\211\211 2:\203\2312@\340>?\205\2352@\341=\203\2222A@\314=\202\235\3422!\202\2352\314=)\204\250\343\n\"\203\257\344\n E\202O
\203\342\346\3476\206\300\317.-\"!\347\316\350 \"!\"\203\342\35105\307\211\211/-3.4&\n\202O\35105\n /-3.4&\n\202O3\211 2:\203=2@\333=\206A2@\332=?\205A2@\336>\20362A@\211 2:\203.2@\333=\20222\314V)\202A\3372!\202A2\314V)\203\2214\211 2:\203\2112@\333=\206\2152@\332=?\205\2152@\336>\203\2022A@\211 2:\203z2@\333=\202~2\314V)\202\215\3372!\202\2152\314V)\204)3\211 2:\203\3252@\332=\206\3312@\333=?\205\3312@\334>\203\3162A@\211 2:\203\3062@\332=\202\3122\314W)\202\331\3352!\202\3312\314W)\203;4\211 2:\203!2@\332=\206%2@\333=?\205%2@\334>\2032A@\211 2:\2032@\332=\2022\314W)\202%\3352!\202%2\314W)\203;\35205-3.4&\202O\35105\307\211\211/-3.4&\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 (bigpos bigneg frac float) (float 6 -1) math-realp (float 123456 -5) math-reject-arg "*Newton's method encountered a singularity" math-float bigneg bigpos (frac float) math-negp (frac float) math-posp (bigpos bigneg) 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 a 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 1 2 3\3154!\203< \3155!\204J \3164!\203\3165!\2036\317=\203\\ \320\321\322\30778F\"\2103T\21132V\203s \323\322\30778F\324\"\2103\nU\203} \3251\326\32787\"1\"\211\211 9:\203\261 9@\330>?\205\265 9@\331=\203\252 9A@\307=\202\265 \3329!\202\265 9\307=)\203\303 \3338\311\"8\202. \334\3354!\3355!\"\203\354 \33670\3377 \"\2117:4\340;!4\202. \30680\3418 \"\2118:5\340;!5\202. 0\203\"\f\20308
5\202\"07
4\336.\202t4\211 9:?\20689@\342>)\204B\3234\343\"\2105\211 9:?\206Q9@\342>)\204[\3235\343\"\21078D45D4\211 9:\203\2519@\344=\206\2559@\345=?\205\2559@\346>\203\2429A@\211 9:\203\2329@\344=\202\2369\307V)\202\255\3159!\202\2559\307V)\307\32787\"\306\211\211 : < = > ? @ A B?T\211?\347X\203r\204rB=A<\326>\350\">=A\203\310\204\310<@\211 9:?\206 9@\342>)\203=@7<@4\351=@>\"\2118:\340;!56\317=\203-\320\3528\"\2105\211 9:?\206<9@\342>)\203O\3165!@=\203O\336\202\346=8=AB\241\210<5<AB\241\210=AA=<AA<\202\346.\2035\211 9:\203\2429@\330>?\205\2469@\331=\203\2339A@\307=\202\246\3329!\202\2469\307=)\203\263\35385E\202\"4\211 9:\203\3359@\330>?\205\3419@\331=\203\3269A@\307=\202\341\3329!\202\3419\307=)\203\356\35374E\202\"C\203 \354;C\306\211\211\2117485&\n\202\"\355;7485%\202\"\323\322\35678F\357\")\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 (bigpos bigneg) float math-zerop calcFunc-incr math-lessp-float math-abs t math-sub math-evaluate-expr math-add (bigpos bigneg frac float) realp bigpos bigneg (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 a 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> @\307=\206A @\310=?\205A @\311>\2038 A@\211:\2031 @\307=\2024 \312V)\202A \313 !\202A \312V)\314\211\211 \315 \316\317\"\320 \"\211\"\204\363
\321 !\211\211:\203\215 @\322>?\205\220 @\323=\203\207 A@\312=\202\220 \324 !\202\220 \312=)\204\363 \325=\203\240 \326\327
\"\210\f\211:\203\332 @\307=\206\335 @\310=?\205\335 @\311>\203\324 A@\211:\203\315 @\307=\202\320 \312V)\202\335 \313 !\202\335 \312V)=\203\353
\f\202M
\f \202M \330
\fE-\207" [high low vhigh a vmid mid math-sub-float bigpos bigneg (frac float) 0 math-posp nil math-nearly-equal math-mul-float (float 5 -1) math-add-float math-evaluate-expr (bigpos bigneg) float math-zerop lots math-do-working "bisect" vec var-DUMMY 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\253 \325\326 !\"!\" # $\327$#\"\306U*\203\244 \330 !\"&\202\257 \313\307\331\"\202\257 \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 b a] 7])
(defalias 'math-find-root #[(expr var guess math-root-widen) "\242\306=\203tGS\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\375 E@\211 F:?\206\351 F@\326>)\204\363 \310E\327\"\210\f\330\f@!\240\210\202\322 \fA\211\203T\f@\242G\236\203 \f\331\f@A@\323\f@8\"\240\210\f\332\f@!\240\210\306C\322\nT\211BW\203L \332\333\f@C\nH\307\334$\206E\310\335\"!C\244\210\202(
C\244\210\202\377 \336E!\337
BE\211\340\n!\203k\341\202o\342\n\343\"&.\207\242\321=\204\202\310\324\"\210\344\"\204\217\310\345\"\210@G\236\203\240\331A@\3238\"\346H\323\\ H\325\347#\307\211 A @\332!\333\350\307\334$ I\307\211 J K\307\211 L MI\203\327\332I!I\330E!E\351E!\203I\203\352IE\211\340E!\203\373\353\202\342\336E!\354\"%\202\335E\211 F:?\206F@\355>)\2033EJEKE@\332!\211LM\356N\202\241E@\357=\203\215\360E!\204H\310E\361\"\210\323E8J\362E8KEA@\363>\203f\364J\365K#JEA@\366>\203x\364K\322J#KJ@\332!LK@\332!M\202\241\367E!\203\233\310\370!\210\202\241\310E\371\"\210L\211 F:\203\313F@\372>?\205\317F@\373=\203\304FA@\374=\202\317\340F!\202\317F\374=)\203\334\306JLE\202\335M\211 F:\203F@\372>?\205\nF@\373=\203\377FA@\374=\202\n\340F!\202\nF\374=)\203\306KME\202\335I\203WL\211 F:?\206+F@\326>)\203WM\211 F:?\206>F@\326>)\203W\375I\307\211\211\211JLKM&\n\202\335L\211 F:\203\237F@\376=\206\243F@\377=?\205\243F@\201O >\203\226FA@\211 F:\203\216F@\376=\202\222F\374V)\202\243\201P F!\202\243F\374V)\203\367M\211 F:\203\357F@\376=\206\363F@\377=?\205\363F@\201O >\203\346FA@\211 F:\203\336F@\376=\202\342F\374V)\202\363\201P F!\202\363F\374V)\204\275L\211 F:\203?F@\377=\206CF@\376=?\205CF@\201Q >\2036FA@\211 F:\203.F@\377=\2022F\374W)\202C\201R F!\202CF\374W)\203\227M\211 F:\203\217F@\377=\206\223F@\376=?\205\223F@\201Q >\203\206FA@\211 F:\203~F@\377=\202\202F\374W)\202\223\201R F!\202\223F\374W)\204\275L\211 F:?\206\246F@\326>)\203\275M\211 F:?\206\271F@\326>)\204\320\201S IJLKM&\202\335\201T JLKM%. !\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 (bigpos bigneg 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-numberp math-newton-root (float 1 6) 100 (bigpos bigneg frac float) point intv math-constp constp 3 (0 1) calcFunc-incr 1 (0 2) math-complexp "*Complex root finder must have derivative" realp (bigpos bigneg) float 0 math-newton-search-root bigpos bigneg var-DUMMY calc-symbolic-mode n math-root-vars symb guess a calc-tweak-eqn-table calc-internal-prec deriv low high vlow vhigh math-root-widen (frac float) math-posp (frac float) math-negp 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:\205 @\305=)\203+ \306 TA\211\203' \307\n H8@L\210\202 )\202- \310\f!\211\211:?\206= @\311>)\203E \312!\207@\313=\203N \207\314\315\"\207" [a m math-min-vars var-DUMMY expr vec -1 2 math-evaluate-expr (bigpos bigneg 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\310 ! \306\211 .\310 !\204 \311 D\310/!\204& \311/D/ G/GU\2042 \312 \210 A/A\211/\2039 @\242\303=\204L \313 @\314\"\210\3150 @\"\204[ \313 @\316\"\210 T1GY\203\200 \317\320\3211G!P! 2\3221\323\3032\211E!\"1)\202[ \3241 H8\306L\210\3241 TH8\306L\210\325/@!\203\324 \3260 @\3271 H\3301 TH\331EE#0\332\333/@!! \324
8\306\211E
A@\306\211\334F.)BB. \324\\\2022 \3260 @1 H#0\335/@!\203\363 \332/@!\306\211E\202./@\242\336=\203(\337/@!\203(\340\341\324/@8\342/@8\"\343\"\332\324/@8!\332\342/@8!/@F\202.\313/@\344\".B. T\2022 .\237.\3450!0 \346U\203\312.@A@\203~3\203h\3470.@A@\324.@8#\202\207\3500.@A@\324.@8\342.@8$\202\207\3470.@@\306# 45 64AA\203\253\3515\211\324\\\\ 5\352\353064$)!\202\257\3114B 7\f\203\304\311\2117A@D\32478E\202\3067+\202E5 6\3515\211\324\\\\ 5\3540 .6$)!\211 7A@ 8\311C 8A\2118\203+\342.@8\203 \351\3558@8A@E!C\244\2108A8.A.\202# 8@C\244\210.A.\202\355\f\203:\311 \32478E\202D\311 A@\32478E,.\207" [calc-symbolic-mode n var-DUMMY var isvec g nil 0 math-vectorp 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 math-realp 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\313 !\203= \314 !\2045 \n\203= \315\316\317 A\"B\207\320 !\204L \n\204L \305 \321\"\210\322 \311\\\315\323A@A\3118A \324$)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" math-vectorp math-constp vec 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\313 !\203= \314 !\2045 \n\203= \315\316\317 A\"B\207\320 !\204L \n\204L \305 \321\"\210\322 \311\\\315\323A@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" math-vectorp math-constp vec 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 !\307!\204 \310!\204 \311\312\"\210\307 !\204* \310 !\204* \311 \312\"\210\313 \"\203: \314\315\n $!\207\316\n \317#\320\211\211\321 \306\n!\322\232\203s \313 \323\"\203_ \202` \324 \325\326\n\327!\327 !\330% ) \331\232\203\233 \313\332\"\203\204 \202\205 \333 !\334\f\325\326\n\327!!\327 !\330%\"!)\335 \"\204\261 \334\f\325\326\n\327!\327 !\320%\"\f,\207" [lo hi expr var sum calc-prefer-frac math-evaluate-expr math-numberp 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) "\306=\203 \307\310 \" \311 !\211 \211 :?\206
@\312>)\204& \313\f\314\"\210\306=\2033 \315\f\316 !\"\f*\207" [mode x var-DUMMY expr res a inf math-div (float 1 0) math-evaluate-expr (bigpos bigneg 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 V WX\204\233 Y\203i YX\306Y\202\233 \311\n!\203\221 \nG\312>\203\213 ZX\313\n8Z\314\n8Y\nA@\202\233 \315 \210\202\233 ZX\306\211ZY\316X!\204\250 \317X\320\"\210XGSXA@GSB\311Z!\204\307 Z\203\307 \321ZDZ\311Y!\204\331 Y\203\331 \321YDYY\204\347 \321\322\nZ\"BYZ\204%YGS
X\203\376 \317Y\323\"\210\202%\324Y!Y\325YG
\n\242\326=\203\307\202\327#Y\233\321 ABZ \306\241\210)ZGS
U\204<ZG
U\204<\317Z\330\"\210YGS\211@\307W\203N\317Y\331\"\210Y A\211\203x @\242\332=\204f\317 @\333\"\210\334 \335\n @\336TZD#\202QSZ A\211\203\245 @\242\332=\204\223\317 @\333\"\210\334 \335\n @\337SZD#\202~S
\\W\203\261\334 \fK\nQ\n\242\326=\204\305\326\337
D\nE\306 [\340\n!\341\342\nD\343\"\313V\340\n!\211\242\344=\203\f\nG\345U\203\f\311\313\n8!\203\f\311\314\n8!\203\f\313\n8G\307V\203\f\314\n8G@TU\204\317R\346\"\210\nA@L\347\313\n8A!M\314\n8JMG?LK\232\2034\306L\321\350?\327\"BC\321\351\324\350?C\"\"BD\347XA! \\\307 ]\352?\327\" ^\327 _\306\211 ` a\306\211 b c\306\211 d e\306\211 f g\306\211 h i\306 j_T\211_BX\2035\327`\306d`
W\203\\`\\`HAI\210\\`H@b`
SU\203\345d\203b\242\353=\203\323\354\355\313b8!d\"dbA@bL\203\356b\357d!\"b\202b\242\353=\203\354\355\313b8!d\206\372\327\"dbA@b\313kS`\\H8bL\210`T\211`\202\223L\203(\340L!c\202/\313K8Jcc\242\353=\203E\355\313c8!dcA@c_\307U\203`dPU\206^ \353=\205^P?OP\203\211d\203\200\360\307d\"]O\203\224]GBG\202\224\317c\361\"\210\202\224d\203\224\317c\361\"\210b\242\353=\203\246\313K8bA@L\210\327`DfCi`?W\203!\340M`H!e^`eI\210\362e]\"eiAifA\211f@g\327\211a`X\203 gA\211g\354g@\362e^aH\"\"\240\210aT\211a\202\344i\354i@\362ec\"\"\240\210`T\211`\202\263O\203\363^cC\"HBH\202\327`Df`?SW\203\203`a`T\211`fA\211f@\233g`D\233haT\211a?W\203<gAghAhg`h@8\240\210\202`. \203\313\364D!\211D\203\236\362DC\"C\202\271\365\366C!!\203\263\367\327CGS\"D\202\271\317Q\370\"\210\311D!\204\323\321\211DDDD\202\323\360CD\"CO\203D\306\211\211 l m n\327NH\203CCnH@m\327lnA\211n\203\354l\362n@m@\"\"lmAm\202\365\355\371m@l\"!lG\203-\362lG@\"l\354Nl\"NGAGHA\211H\204\351+\324C!\211I nD o\307 p\314 V\327 ` \203zP\204zB?X\203o\327\202x\360NB?Z\"p \203\247nA\211n\203\247n\356n@\357\362`T\211`oA\211o@8p\"!\"\240\210\202~\340J!Iq\203\343IY\306 o nnAnoA\211o\203\333\326o@n@EB\202\275\321\237Bq+.\372Q
\306 \211 B\307Y\203\337BDBBS\211B\202\357@\211B\307Y\203\"\336BDBBS\211B\202 *\363IAZA\"#\373U\203\214\321\nCDJ\327 BBTB A\211\203_ @\242\374=\203_ @A@B=\204A \203hJ\202j\321C*NP\203\206B?V\203\206\375NB?ZE\202\207\376\257\202\215\n!. \207" [math-dummy-counter mode expr p dummy v nil 1 full math-vectorp (3 4) 2 3 math-dimension-error math-matrixp math-reject-arg matrixp vec 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\f @\302>)\203 \303\207@\304=\203 A@\207\305\306\307\310A\"\"\207" [expr a (bigpos bigneg frac float cplx polar hms date mod var) 0 calcFunc-fitvar apply max mapcar calcFunc-hasfitvars] 6])
(defalias 'calcFunc-hasfitparams #[(expr) "\211:?\206\f @\302>)\203 \303\207@\304=\203 A@\207\305\306\307\310A\"\"\207" [expr a (bigpos bigneg frac float cplx polar hms date mod var) 0 calcFunc-fitparam apply max mapcar calcFunc-hasfitparams] 6])
(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\f @\304>)\2034 \242\305=\205C \306!\2063 \307\n\"\211\203, \211AT\241\207\310B\nB\211\207A\211\205C \311@!\210\2024 \207" [expr a math-all-vars-vars math-all-vars-found (bigpos bigneg frac float cplx polar hms date mod var) var math-const-var assoc 1 math-all-vars-rec] 4])
(provide 'calcalg3)