@mathzetta/inlinemath 0.0.0-stage → 0.1.0

package/dist/cas.js

@@ -0,0 +1,2 @@
+/* Inline Math 0.1.0 — MIT licence; includes mathzetta (and @mathzetta/cas, for the CAS) */
+import{m as e,o as n,q as t,r,u as s,G as o,v as i,w as a,x as u,y as l,i as c,z as f,A as p,f as d,B as h,C as m}from"./chunks/lite-Bj_L6w2g.js";import{c as g}from"./chunks/calculate-BX9lJ7Ff.js";import{a as v}from"./chunks/html-DWtfytjA.js";const b=new Map,y=-1,w="\\:",x=e=>"string"==typeof e&&e.startsWith(w)&&e.length>2;function q(e){if(!e||!e.symbol)throw new Error("operator needs a symbol");if(!Number.isInteger(e.arity)||e.arity<0&&e.arity!==y)throw new Error(`operator "${e.symbol}" needs an integer arity >= 0, or VARIADIC`);if(e.commutative&&2!==e.arity)throw new Error(`only binary operators can be commutative ("${e.symbol}")`);const n={symbol:e.symbol,arity:e.arity,commutative:!!e.commutative,fold:"function"==typeof e.fold?e.fold:null,label:e.label||e.symbol};return b.set(n.symbol,n),n}const k=Object.freeze({symbol:w,arity:y,commutative:!1,fold:null,label:"function"});function $(e){const n=b.get(e);return n||(x(e)?{...k,symbol:e,label:e.slice(2)}:null)}function S(e){return b.has(e)||x(e)}const O=e=>Number.isFinite(e)&&Number.isInteger(e)?e:null,N=e=>Number.isFinite(e)?e:null;!function(){q({symbol:"+",arity:2,commutative:!0,label:"add",fold:(e,n)=>N(e+n)}),q({symbol:"-",arity:2,label:"subtract",fold:(e,n)=>N(e-n)}),q({symbol:"*",arity:2,commutative:!0,label:"multiply",fold:(e,n)=>N(e*n)}),q({symbol:"/",arity:2,label:"divide",fold:(e,n)=>0===n?null:O(e/n)}),q({symbol:"^",arity:2,label:"power",fold:(e,n)=>Number.isInteger(n)&&Math.abs(n)<=64?N(e**n):null}),q({symbol:"neg",arity:1,label:"negate",fold:e=>N(-e)}),q({symbol:"pm",arity:1,label:"plus or minus"}),q({symbol:"±",arity:2,label:"plus or minus"}),q({symbol:"mod",arity:2,label:"modulo",fold:(e,n)=>0===n?null:O((e%n+n)%n)}),q({symbol:"!",arity:1,label:"factorial",fold:e=>Number.isInteger(e)&&e>=0&&e<=18?[...Array(e).keys()].reduce((e,n)=>e*(n+1),1):null}),q({symbol:"_",arity:2,label:"subscript"}),q({symbol:"der",arity:2,label:"derivative"}),q({symbol:"int",arity:y,label:"integral"}),q({symbol:"lim",arity:y,label:"limit"}),q({symbol:"\\sumidx",arity:y,label:"sum over an index"}),q({symbol:"\\prodidx",arity:y,label:"product over an index"}),q({symbol:"abs",arity:1,label:"absolute value",fold:e=>N(Math.abs(e))}),q({symbol:"floor",arity:1,label:"floor",fold:e=>N(Math.floor(e))}),q({symbol:"ceil",arity:1,label:"ceiling",fold:e=>N(Math.ceil(e))}),q({symbol:"round",arity:2,label:"round",fold:(e,n)=>Number.isInteger(n)?N(Math.round(e*10**n)/10**n):null}),q({symbol:"fixed",arity:2,label:"fixed",fold:()=>null}),q({symbol:"sqrt",arity:1,label:"square root",fold:e=>O(Math.sqrt(e))}),q({symbol:"root",arity:2,label:"nth root",fold:(e,n)=>O(Math.round(n**(1/e))**e===n?Math.round(n**(1/e)):NaN)});for(const e of["sin","cos","tan","csc","sec","cot","sinh","cosh","tanh","arcsin","arccos","arctan","arcsinh","arccosh","arctanh","ln","loglog","logloglog","expm1"])q({symbol:e,arity:1,label:e});q({symbol:"log",arity:2,label:"log to a base"}),q({symbol:"atan2",arity:2,label:"atan2"}),q({symbol:"perm",arity:2,label:"permutations"}),q({symbol:"comb",arity:2,label:"combinations"}),q({symbol:"random",arity:0,label:"random"});for(const e of["max","min","gcd","mean","sum","prod"])q({symbol:e,arity:y,label:e});for(const e of["()","[]","{}","(]","[)"])q({symbol:e,arity:y,label:"group"});for(const e of["=","===","≠","≠=","<",">","≤","≥","∈","∉","∪","∩","≈","≡","→","pmod","&","|","\\xor","\\nand","\\nor"])q({symbol:e,arity:2,label:e});q({symbol:"\\not",arity:1,label:"not"}),q({symbol:"?:",arity:3,label:"if then else"})}();const A="num",I="sym",M="op",E=Object.freeze([]);class T extends Error{constructor(e){super(e),this.name="ExprError"}}function j(e){if("number"!=typeof e||!Number.isFinite(e))throw new T(`not a finite number: ${e}`);const n=(t=e,Object.is(t,-0)?"0":String(t));var t;return Object.freeze({type:A,token:n,value:e,args:E,s:n})}function Z(e){if("string"!=typeof e||""===e)throw new T("symbol needs a non-empty token");return Object.freeze({type:I,token:e,value:null,args:E,s:e})}function z(e){const n=/^(-?)0*(\d+?)(?:\.(\d*?)0*)?$/.exec(String(e));if(!n)throw new T(`not an exact decimal: ${e}`);const[,t,r,s=""]=n,o=s?`${r}.${s}`:r,i=/^0(\.0*)?$/.test(o)?"0":t+o;return Object.freeze({type:A,token:i,value:Number(i),args:E,s:i})}const D=/^-?\d+(\.\d+)?$/;function B(e){const n=String(e);if(D.test(n))return z(n);if(""!==n&&!/\s/.test(n)){const e=Number(n);if(Number.isFinite(e))return j(e)}return Z(n)}function C(e,n){const t=$(e);if(!t)throw new T(`unknown operator "${e}"`);if(t.arity!==y&&n.length!==t.arity)throw new T(`"${e}" takes ${t.arity} operand(s), got ${n.length}`);const r=Object.freeze(n.slice()),s=t.arity===y?` ${r.length}`:"",o=`${r.length?r.map(e=>e.s).join(" "):""}${r.length,""}${s} ${e}`.trimStart();return Object.freeze({type:M,op:e,token:e,value:null,args:r,s:o})}const F=e=>e.type===A,V=e=>e.type===M;function R(e){return U(function(e){const n=String(e).trim();return""===n?[]:n.split(/\s+/)}(e))}function U(e){if(!Array.isArray(e))throw new T("fromPostfix needs an array of tokens");if(0===e.length)throw new T("empty expression");const n=[];for(let t=0;t<e.length;t++){const r=String(e[t]);if(S(r)){let e=$(r).arity;if(e===y){const s=n.pop();if(e=s&&s.type===A&&Number.isInteger(s.value)&&s.value>=0?s.value:-1,e<0)throw new T(`"${r}" at position ${t+1} takes any number of operands, so a count must come before it`)}if(n.length<e)throw new T(`"${r}" at position ${t+1} needs ${e} operand(s), only ${n.length} on the stack`);const s=n.splice(n.length-e,e);n.push(C(r,s))}else n.push(B(r))}if(1!==n.length){const t=[...new Set(e.filter(e=>!S(e)))];throw new T(`expression leaves ${n.length} values on the stack; it should leave exactly 1. Read as operands: ${t.join(", ")} — is one of those meant to be an operator?`)}return n[0]}function P(e){if(e.type!==M)return 0;let n=0;for(const t of e.args)n=Math.max(n,P(t));return n+1}function W(e){if(e.type!==M)return 1;let n=1;for(const t of e.args)n+=W(t);return n}function L(e,n=[]){if(e.type!==M)n.push(e.token);else for(const t of e.args)L(t,n);return n}function _(e){return e.type===A||e.type!==I&&e.args.every(_)}function H(e,n="pre",t=[],r=[]){"pre"===n&&r.push({node:e,path:t});for(let s=0;s<e.args.length;s++)H(e.args[s],n,t.concat(s),r);return"post"===n&&r.push({node:e,path:t}),r}function J(e,n,t){if(0===n.length)return t;const[r,...s]=n,o=e.args.slice();return o[r]=J(o[r],s,t),C(e.op,o)}class Q extends Error{constructor(e){super(e),this.name="CasError"}}const Y=1e5,G=e=>e<0n?-e:e,K=e=>0n===e?0:G(e).toString(2).length;function X(e,n){for(e=G(e),n=G(n);n;)[e,n]=[n,e%n];return e}class ee{constructor(e,n=1n){if(0n===n)throw new Q("division by zero");n<0n&&(e=-e,n=-n);const t=X(e,n)||1n;if(this.n=e/t,this.d=n/t,K(this.n)>Y||K(this.d)>Y)throw new Q("a number too large to keep exactly")}static of(e){return new ee(BigInt(e))}add(e){return new ee(this.n*e.d+e.n*this.d,this.d*e.d)}sub(e){return new ee(this.n*e.d-e.n*this.d,this.d*e.d)}mul(e){return new ee(this.n*e.n,this.d*e.d)}div(e){if(0n===e.n)throw new Q("division by zero");return new ee(this.n*e.d,this.d*e.n)}neg(){return new ee(-this.n,this.d)}abs(){return this.n<0n?this.neg():this}inv(){return ee.ONE.div(this)}pow(e){if((e=BigInt(e))<0n)return this.inv().pow(-e);if((0n===this.n||1n===G(this.n))&&1n===this.d)return 0n===e||-1n===this.n&&e%2n==0n?ee.ONE:this;if(BigInt(Math.max(K(this.n),K(this.d)))*e>BigInt(Y))throw new Q("a power too large to work out exactly");return new ee(this.n**e,this.d**e)}get isZero(){return 0n===this.n}get isOne(){return 1n===this.n&&1n===this.d}get isInteger(){return 1n===this.d}get sign(){return this.n<0n?-1:this.n>0n?1:0}equals(e){return this.n===e.n&&this.d===e.d}compare(e){const n=this.n*e.d,t=e.n*this.d;return n<t?-1:n>t?1:0}floor(){const e=this.n/this.d;return new ee(this.n<0n&&e*this.d!==this.n?e-1n:e)}ceil(){return this.neg().floor().neg()}toNumber(){const e=Number(this.n),n=Number(this.d);if(Number.isFinite(e)&&Number.isFinite(n))return e/n;const t=BigInt(Math.max(K(this.n),K(this.d))-1e3);return Number(this.n>>t)/Number(this.d>>t)}toString(){return 1n===this.d?String(this.n):`${this.n}/${this.d}`}}function ne(e){const n=/^(-?)(\d*)(?:\.(\d*))?(?:[eE]([+-]?\d+))?$/.exec(String(e));if(!n||""===n[2]&&!n[3])return null;const[,t,r,s="",o="0"]=n;let i=new ee(BigInt((r||"0")+s),10n**BigInt(s.length));const a=Number(o);return a&&(i=i.mul(ee.of(10).pow(a))),t?i.neg():i}function te(e,n){if(e<2n)return e;if(n>=K(e))return 1n;const t=BigInt(n);let r=1n<<BigInt(Math.ceil(K(e)/n));for(;;){const n=((t-1n)*r+e/r**(t-1n))/t;if(n>=r)break;r=n}for(;r**t>e;)r--;for(;(r+1n)**t<=e;)r++;return r}ee.ZERO=new ee(0n),ee.ONE=new ee(1n),ee.MINUS_ONE=new ee(-1n),ee.HALF=new ee(1n,2n);const re=(()=>{const e=[],n=new Uint8Array(1e4);for(let t=2;t<n.length;t++)if(!n[t]){e.push(BigInt(t));for(let e=t*t;e<n.length;e+=t)n[e]=1}return e})();function se(e,n){const t=BigInt(n);let r=1n,s=1n,o=e;for(const e of re){if(e*e>o)break;if(o%e)continue;let n=0n;for(;o%e===0n;)o/=e,n++;r*=e**(n/t),s*=e**(n%t)}if(o>1n){const e=function(e,n){if(e<0n)return null;if(e<2n)return e;const t=te(e,n);return t**BigInt(n)===e?t:null}(o,n);null!==e?r*=e:s*=o}return{outside:r,inside:s}}const oe=new Set(["pi","e","π"]),ie={complex:!1};function ae(e,n){const t=ie.complex;ie.complex=Boolean(e);try{return n()}finally{ie.complex=t}}const ue=e=>oe.has(e)||ie.complex&&"i"===e,le=e=>ie.complex&&"atom"===e.k&&"sym"===e.node.type&&"i"===e.node.token,ce=e=>({k:"num",q:e}),fe=ce(ee.ONE),pe=ce(ee.ZERO),de=ce(ee.HALF),he=e=>"num"===e.k,me=e=>({k:"atom",node:e}),ge=new WeakMap;function ve(e){let n=ge.get(e);return void 0===n&&ge.set(e,n=Qe(e).s),n}function be(e){if("num"===e.k)return 0;if("atom"===e.k)return"sym"===e.node.type&&ue(e.node.token)?0:1;if("add"===e.k)return Math.max(...e.terms.map(e=>be(e.x)));let n=0;for(const t of e.fs){const e=be(t.b);e&&(n+=e*(he(t.e)?t.e.q.toNumber():1))}return n}function ye(e,n){return"num"===e.k?n&&!he(n)?3:0:le(e)?1.5:"atom"===e.k&&"sym"===e.node.type?ue(e.node.token)?n&&!he(n)?3:1:2:"atom"===e.k?3:4}const we=(e,n)=>e<n?-1:e>n?1:0;function xe(e){const n=new Map,t=(e,t)=>{"atom"!==e.k||"sym"!==e.node.type||ue(e.node.token)||n.set(e.node.token,(n.get(e.node.token)??0)+t)};if("mul"===e.k)for(const n of e.fs)t(n.b,he(n.e)?n.e.q.toNumber():1);else t(e,1);return n}function qe(e,n){const t=be(n.x)-be(e.x);if(t)return t;const r=xe(e.x),s=xe(n.x);for(const e of[...new Set([...r.keys(),...s.keys()])].sort()){const n=(s.get(e)??0)-(r.get(e)??0);if(n)return n}return we(ve(e.x),ve(n.x))}function ke(e,n){return e.isOne?n:e.isZero?pe:Ne([ce(e),n])}function $e(e){if("num"===e.k)return{coef:e.q,rest:null};if("mul"!==e.k||e.c.isOne)return{coef:ee.ONE,rest:e};const n=1===e.fs.length&&he(e.fs[0].e)&&e.fs[0].e.q.isOne?e.fs[0].b:{k:"mul",c:ee.ONE,fs:e.fs};return{coef:e.c,rest:n}}function Se(e){let n=ee.ZERO;const t=new Map,r=(e,n)=>{const r=ve(n),s=t.get(r);s?s.coef=s.coef.add(e):t.set(r,{coef:e,x:n})};for(const t of e)if("num"===t.k)n=n.add(t.q);else if("add"===t.k){n=n.add(t.c);for(const e of t.terms)r(e.coef,e.x)}else{const{coef:e,rest:n}=$e(t);r(e,n)}const s=[...t.values()].filter(e=>!e.coef.isZero);return s.length?n.isZero&&1===s.length?ke(s[0].coef,s[0].x):(s.sort(qe),{k:"add",c:n,terms:s}):ce(n)}function Oe(e,n){if(!he(n))return e.isOne?{coef:ee.ONE,fs:[]}:{coef:ee.ONE,fs:[{b:ce(e),e:n}]};const t=n.q;if(t.isInteger)return{coef:e.pow(t.n),fs:[]};if(e.isZero){if(t.sign<0)throw new Q("division by zero");return{coef:ee.ZERO,fs:[]}}if(e.isOne)return{coef:ee.ONE,fs:[]};if(e.sign<0){if(t.d%2n==1n){const{coef:r,fs:s}=Oe(e.neg(),n);return{coef:t.n%2n==0n?r:r.neg(),fs:s}}if(ie.complex&&2n===t.d){const{coef:r,fs:s}=Oe(e.neg(),n);return{coef:r,fs:[...s,{b:Me,e:ce(new ee(t.n))}]}}return{coef:ee.ONE,fs:[{b:ce(e),e:n}]}}if(!e.isInteger){const r=Oe(new ee(e.n),n),s=Oe(new ee(e.d),ce(t.neg()));return{coef:r.coef.mul(s.coef),fs:[...r.fs,...s.fs]}}const r=t.floor(),s=t.sub(r);if(BigInt(K(e.n))*s.n>BigInt(Y)||s.d>1000n)return{coef:ee.ONE,fs:[{b:ce(e),e:n}]};let o=r.isZero?ee.ONE:e.pow(r.n);const{outside:i,inside:a}=se(e.n**s.n,Number(s.d));o=o.mul(new ee(i));return{coef:o,fs:a>1n?[{b:ce(new ee(a)),e:ce(new ee(1n,s.d))}]:[]}}function Ne(e){let n=ee.ONE;const t=new Map,r=(e,n)=>{const r=ve(e),s=t.get(r);s?s.e=Se([s.e,n]):t.set(r,{b:e,e:n})};for(const t of e)if("num"===t.k)n=n.mul(t.q);else if("mul"===t.k){n=n.mul(t.c);for(const e of t.fs)r(e.b,e.e)}else r(t,fe);if(n.isZero)return pe;const s=[],o=new Map;for(const e of t.values())if(!he(e.e)||!e.e.q.isZero){if("num"===e.b.k){if(he(e.e)&&!e.e.q.isInteger&&e.b.q.sign>0&&e.b.q.isInteger){const n=e.e.q.toString(),t=o.get(n);t?t.n*=e.b.q.n:o.set(n,{n:e.b.q.n,e:e.e});continue}const t=Oe(e.b.q,e.e);n=n.mul(t.coef),s.push(...t.fs);continue}if("mul"===e.b.k&&he(e.e)&&e.e.q.isInteger){const t=Ae(e.b,e.e);"mul"===t.k?(n=n.mul(t.c),s.push(...t.fs)):"num"===t.k?n=n.mul(t.q):s.push({b:t,e:fe});continue}s.push(e)}for(const{n:e,e:t}of o.values()){const r=Oe(new ee(e),t);n=n.mul(r.coef),s.push(...r.fs)}if(n.isZero)return pe;const i=new Map;for(const e of s){const n=ve(e.b),t=i.get(n);t?t.e=Se([t.e,e.e]):i.set(n,{...e})}let a=[...i.values()].filter(e=>!(he(e.e)&&e.e.q.isZero));if(ie.complex){a=a.flatMap(e=>{if(!le(e.b)||!he(e.e)||!e.e.q.isInteger)return[e];const t=(e.e.q.n%4n+4n)%4n;return t>=2n&&(n=n.neg()),t%2n?[{b:e.b,e:fe}]:[]});const e=a.filter(e=>"add"===e.b.k&&he(e.e)&&e.e.q.isInteger&&Te(e.b)),t=a.some(e=>le(e.b)),r=e.every(e=>(e.e.q.sign<0?-e.e.q.n:e.e.q.n)<=60n);if(e.length&&r&&(e.length>1||t||!e[0].e.q.isOne)){const r=a.filter(n=>!e.includes(n)&&!le(n.b));let s=t?[Me]:[fe];for(const n of e){const e=n.e.q.sign<0?ze(n.b):n.b;for(let t=0n;t<(n.e.q.sign<0?-n.e.q.n:n.e.q.n);t++)s=je(s,Xe(e))}return Ne([ce(n),Se(s),...r.map(e=>Ie(e.b,e.e))])}}if(a.some(e=>"num"===e.b.k&&he(e.e)&&e.e.q.compare(ee.ONE)>=0))return Ne([ce(n),...a.map(e=>"num"===e.b.k?Ae(e.b,e.e):{k:"mul",c:ee.ONE,fs:[e]})]);if(!a.length)return ce(n);if(1===a.length&&!n.isOne&&"add"===a[0].b.k&&he(a[0].e)&&a[0].e.q.isOne){const e=a[0].b;return Se([ce(n.mul(e.c)),...e.terms.map(e=>ke(n.mul(e.coef),e.x))])}return n.isOne&&1===a.length&&he(a[0].e)&&a[0].e.q.isOne?a[0].b:(a=a.sort((e,n)=>ye(e.b,e.e)-ye(n.b,n.e)||we(ve(e.b),ve(n.b))),{k:"mul",c:n,fs:a})}function Ae(e,n){if(ie.complex&&he(n)&&n.q.isInteger){if(le(e))return Ne([{k:"mul",c:ee.ONE,fs:[{b:e,e:n}]}]);if("add"===e.k&&Te(e)&&!n.q.isOne)return Ne([{k:"mul",c:ee.ONE,fs:[{b:e,e:n}]}])}if(he(n)){if(n.q.isZero){if(he(e)&&e.q.isZero)throw new Q("0^0 is undefined");return fe}if(n.q.isOne)return e}if(he(e)){const t=Oe(e.q,n);if(!t.fs.length)return ce(t.coef);if(t.fs.some(e=>le(e.b)))return Ne([ce(t.coef),...t.fs.map(e=>({k:"mul",c:ee.ONE,fs:[e]}))]);const r={k:"mul",c:t.coef,fs:t.fs};return t.coef.isOne&&1===t.fs.length&&he(t.fs[0].e)&&t.fs[0].e.q.isOne?t.fs[0].b:r}if(ie.complex&&he(n)&&2n===n.q.d&&Te(e)){const t=Ze(e);if(t&&he(t.im)&&!t.im.q.isZero){const e=e=>Ae(e,ce(ee.of(2))),r=Ae(Se([e(t.re),e(t.im)]),de),s=Ae(Ne([Se([r,t.re]),ce(ee.HALF)]),de),o=Ae(Ne([Se([r,Ce(t.re)]),ce(ee.HALF)]),de);return Ae(Se([s,Ne([ce(ee.of(t.im.q.sign)),o,Me])]),ce(new ee(n.q.n)))}}if("mul"===e.k){if(he(n)&&n.q.isInteger)return Ne([ce(e.c.pow(n.q.n)),...e.fs.map(e=>Ie(e.b,Ne([e.e,n])))]);if(!e.c.isOne&&e.c.sign>0)return Ne([Ae(ce(e.c),n),Ae({k:"mul",c:ee.ONE,fs:e.fs},n)]);if(e.c.isOne&&1===e.fs.length){const[t]=e.fs;if(he(t.e)&&he(n)&&1n===t.e.q.n&&t.e.q.d%2n==1n)return Ie(t.b,Ne([t.e,n]))}}return{k:"mul",c:ee.ONE,fs:[{b:e,e:n}]}}function Ie(e,n){return he(n)&&n.q.isZero?fe:he(e)?Ae(e,n):he(n)&&n.q.isOne?e:{k:"mul",c:ee.ONE,fs:[{b:e,e:n}]}}const Me=me(Z("i")),Ee=e=>ie.complex&&Qe(e).s.split(" ").includes("i");function Te(e){return ie.complex&&Ee(e)&&0===Ke(Qe(e)).size}function je(e,n){const t=[];for(const r of e)for(const e of n)t.push(Ne([r,e]));return Xe(Se(t))}function Ze(e){const n=[],t=[];for(const r of Xe(e)){const{coef:e,rest:s}=$e(r),o=s?en(s):[],i=o.filter(e=>le(e)),a=o.filter(e=>!le(e));if(a.some(e=>Ee(e)))return null;(i.length?t:n).push(Ne([ce(e),...a]))}return{re:Se(n),im:Se(t)}}function ze(e){const n=Ze(e);if(!n)return Ae(e,ce(ee.MINUS_ONE));const t=Se([Ae(n.re,ce(ee.of(2))),Ae(n.im,ce(ee.of(2)))]);return Ne([Se([n.re,Ne([ce(ee.MINUS_ONE),n.im,Me])]),Ae(t,ce(ee.MINUS_ONE))])}const De={re:([e])=>Ze(e)?.re??null,im:([e])=>Ze(e)?.im??null,conj:([e])=>{const n=Ze(e);return n?Se([n.re,Ne([ce(ee.MINUS_ONE),n.im,Me])]):null},arg:([e])=>{const n=Ze(e);return n?Ue("atan2",[n.im,n.re]):null},abs:([e])=>{const n=Ze(e);return!n||he(n.im)&&n.im.q.isZero?null:Ae(Se([Ae(n.re,ce(ee.of(2))),Ae(n.im,ce(ee.of(2)))]),de)}},Be={Re:"re",re:"re",Im:"im",im:"im",conj:"conj",overline:"conj",bar:"conj",arg:"arg"},Ce=e=>Ne([ce(ee.MINUS_ONE),e]),Fe=(e,n)=>Ne([e,Ae(n,ce(ee.MINUS_ONE))]),Ve=3000n,Re={abs:([e])=>he(e)?ce(e.q.abs()):"mul"!==e.k||e.c.isOne?null:Ne([ce(e.c.abs()),Ue("abs",[$e(e).rest])]),floor:([e])=>he(e)?ce(e.q.floor()):null,ceil:([e])=>he(e)?ce(e.q.ceil()):null,"!":([e])=>{if(!he(e)||!e.q.isInteger||e.q.sign<0||e.q.n>Ve)return null;let n=1n;for(let t=2n;t<=e.q.n;t++)n*=t;return ce(new ee(n))},mod:([e,n])=>he(e)&&he(n)&&!n.q.isZero?ce(e.q.sub(n.q.mul(e.q.div(n.q).floor()))):null,max:e=>e.every(he)&&e.length?e.reduce((e,n)=>n.q.compare(e.q)>0?n:e):null,min:e=>e.every(he)&&e.length?e.reduce((e,n)=>n.q.compare(e.q)<0?n:e):null,gcd:e=>e.length&&e.every(e=>he(e)&&e.q.isInteger)?ce(new ee(e.reduce((e,n)=>X(e,n.q.n),0n))):null,log:([e,n])=>{if(!he(e)||!he(n)||e.q.sign<=0||n.q.sign<=0||e.q.isOne)return null;if(n.q.isOne)return pe;const t=Math.log(n.q.toNumber())/Math.log(e.q.toNumber());if(!Number.isFinite(t))return null;for(let r=1;r<=12;r++){const s=Math.round(t*r);if(!(Math.abs(t*r-s)>1e-9||0===s)&&(Math.abs(s)<=256&&e.q.pow(BigInt(s)).equals(n.q.pow(BigInt(r)))))return ce(new ee(BigInt(s),BigInt(r)))}return null},comb:([e,n])=>he(e)&&he(n)&&e.q.isInteger&&n.q.isInteger&&e.q.sign>=0&&n.q.sign>=0&&e.q.n<=Ve?ce(function(e,n){if(n>e)return ee.ZERO;n>e-n&&(n=e-n);let t=1n;for(let r=1n;r<=n;r++)t=t*(e-n+r)/r;return new ee(t)}(e.q.n,n.q.n)):null,perm:([e,n])=>{if(!(he(e)&&he(n)&&e.q.isInteger&&n.q.isInteger&&e.q.sign>=0&&n.q.sign>=0&&e.q.n<=Ve))return null;if(n.q.n>e.q.n)return pe;let t=1n;for(let r=0n;r<n.q.n;r++)t*=e.q.n-r;return ce(new ee(t))}};function Ue(e,n){if(ie.complex&&De[e]){const t=De[e](n);if(t)return t}const t=Re[e];if(t){const e=t(n);if(e)return e}return me(C(e,n.map(Qe)))}const Pe=new Map;function We(e){const n=ie.complex?`ℂ ${e.s}`:e.s,t=Pe.get(n);if(t)return t;const r=function(e){if("num"===e.type){const n=ne(e.token);if(!n)throw new Q(`not an exact number: ${e.token}`);return ce(n)}if("sym"===e.type)return"π"===e.token?me(Z("pi")):me(e);const n=e.args.map(We);switch(e.op){case"+":return Se(n);case"-":return((e,n)=>Se([e,Ce(n)]))(n[0],n[1]);case"neg":return Ce(n[0]);case"*":return Ne(n);case"/":return Fe(n[0],n[1]);case"^":return Ae(n[0],n[1]);case"sqrt":return Ae(n[0],de);case"root":return Ae(n[1],Ae(n[0],ce(ee.MINUS_ONE)));default:{const t=e.op.startsWith("\\:")?e.op.slice(2):null;return t&&ie.complex&&Be[t]&&1===n.length?Ue(Be[t],n):t&&S(t)&&$(t).arity===n.length?Ue(t,n):Ue(e.op,n)}}}(e);return Pe.size>2e4&&Pe.clear(),Pe.set(n,r),r}const Le=e=>Qe(We(e)),_e=e=>z(String(e));function He(e){return e.isInteger?_e(e.n):C("/",[_e(e.n),_e(e.d)])}const Je=new WeakMap;function Qe(e){let n=Je.get(e);return n||Je.set(e,n=function(e){switch(e.k){case"num":return He(e.q);case"atom":return e.node;case"mul":return function(e,n){const t=[],r=[];for(const e of n)he(e.e)&&e.e.q.sign<0?r.push(Ye(e.b,ce(e.e.q.neg()))):t.push(Ye(e.b,e.e));let s=1n!==e.n&&-1n!==e.n||!t.length?Ge([_e(e.n),...t]):Ge(t);-1n===e.n&&t.length&&(s=C("neg",[s]));1n!==e.d&&r.unshift(_e(e.d));return r.length?C("/",[s,Ge(r)]):s}(e.c,e.fs);case"add":return function(e){let n=null;const t=e=>ie.complex&&Ee(e.x)&&0===Ke(Qe(e.x)).size,r=[...e.terms.filter(e=>!t(e)),...e.terms.filter(t)],s=e.c,o=()=>{s.isZero||(n=n?s.sign<0?C("-",[n,He(s.neg())]):C("+",[n,He(s)]):He(s))};let i=!1;for(const e of r)!i&&t(e)&&(o(),i=!0),n=n?e.coef.sign<0?C("-",[n,Qe(ke(e.coef.neg(),e.x))]):C("+",[n,Qe(ke(e.coef,e.x))]):Qe(ke(e.coef,e.x));i||o();return n}(e);default:throw new Error(`bad form ${e.k}`)}}(e)),n}function Ye(e,n){const t=Qe(e);if(he(n)){const e=n.q;if(e.isOne)return t;if(1n===e.n&&2n===e.d)return C("sqrt",[t]);if(1n===e.n)return C("root",[_e(e.d),t])}return C("^",[t,Qe(n)])}const Ge=e=>e.reduce((e,n)=>e?C("*",[e,n]):n,null);function Ke(e,n=new Set){if("sym"===e.type)ue(e.token)||n.add(e.token);else if(V(e))for(const t of e.args)Ke(t,n);return n}function Xe(e){if("add"!==e.k)return[e];const n=e.terms.map(e=>ke(e.coef,e.x));return e.c.isZero||n.push(ce(e.c)),n}function en(e){if("num"===e.k)return[e];if("mul"!==e.k)return[e];const n=e.c.isOne?[]:[ce(e.c)];for(const t of e.fs)n.push(Ie(t.b,t.e));return n}const nn="_",tn=new WeakMap;function rn(e,n){return[...n,...$(e.op)?.arity===y?[String(e.args.length)]:[],e.op].join(" ")}function sn(e,n){let t=tn.get(e);t||tn.set(e,t=new Map);const r=t.get(n);if(r)return r;const s=function(e,n){const t={key:nn,slots:[e]};if(e.type!==M||n<=0)return[t];const r=e.args.map(e=>sn(e,n-1)),s=[t],o=new Array(r.length).fill(0);for(;;){const n=[],t=[];for(let e=0;e<r.length;e++){const s=r[e][o[e]];n.push(s.key);for(const e of s.slots)t.push(e)}if(s.push({key:rn(e,n),slots:t}),s.length>=4096)break;let i=r.length-1;for(;i>=0&&!(++o[i]<r[i].length);i--)o[i]=0;if(i<0)break}return s}(e,n);return t.set(n,s),s}function on(e){return e.type!==M?nn:rn(e,e.args.map(on))}function an(e,n=[]){if(e.type!==M)n.push(e);else for(const t of e.args)an(t,n);return n}const un=new Map,ln=new Map;class cn extends Error{constructor(e){super(e),this.name="RequirementError"}}function fn(e,n,t=""){un.set(e,{name:e,fn:n,doc:t,arity:Math.max(0,n.length-1)})}function pn(e,n,t=""){ln.set(e,{name:e,fn:n,doc:t})}function dn(e){return un.get(e)||null}function hn(e){return ln.get(e)||null}const mn=e=>e.type===A;function gn(e,n){return e.type===M&&(e.op===n||e.args.some(e=>gn(e,n)))}fn("always",()=>!0,"matches anything"),fn("never",()=>!1,"matches nothing"),fn("isLeaf",e=>e.type!==M,"a bare operand, not an operation"),fn("isCompound",e=>e.type===M,"an operation, not a bare operand"),fn("isSymbol",e=>e.type===I,"a non-numeric operand such as x"),fn("isOperation",(e,n)=>e.type===M&&e.op===n,"is this specific operation"),fn("notOperation",(e,n)=>!(e.type===M&&e.op===n),"is anything except this operation"),fn("isNumber",mn,"a numeric literal"),fn("isInteger",e=>mn(e)&&Number.isInteger(e.value),"a whole number"),fn("isPositive",e=>mn(e)&&e.value>0,"greater than zero"),fn("isNegative",e=>mn(e)&&e.value<0,"less than zero"),fn("isNonNegative",e=>mn(e)&&e.value>=0,"zero or greater"),fn("isZero",e=>mn(e)&&0===e.value,"exactly 0"),fn("isOne",e=>mn(e)&&1===e.value,"exactly 1"),fn("isNonZero",e=>mn(e)&&0!==e.value,"a number that is provably not 0"),fn("isNotZeroLiteral",e=>!(mn(e)&&0===e.value),"anything except a literal 0"),fn("isPositiveInteger",e=>mn(e)&&Number.isInteger(e.value)&&e.value>0,"a whole number above 0"),fn("isNegativeInteger",e=>mn(e)&&Number.isInteger(e.value)&&e.value<0,"a whole number below 0"),fn("isEven",e=>mn(e)&&Number.isInteger(e.value)&&e.value%2==0,"a whole even number"),fn("isOdd",e=>mn(e)&&Number.isInteger(e.value)&&1===Math.abs(e.value%2),"a whole odd number"),fn("greaterThan",(e,n)=>mn(e)&&e.value>Number(n),"above a given value"),fn("lessThan",(e,n)=>mn(e)&&e.value<Number(n),"below a given value"),fn("isValue",(e,n)=>mn(e)&&e.value===Number(n),"equals a given number"),fn("isNumeric",_,"contains no symbols: a constant (it need not fold to a number: 1 sin)"),fn("hasSymbol",e=>!_(e),"mentions a symbol somewhere, so it is not just a number"),fn("isVar",(e,n)=>e.type===I&&e.token===String(n),"is exactly this symbol"),fn("contains",(e,n)=>L(e).includes(String(n)),"mentions this operand somewhere"),fn("isFree",(e,n)=>!L(e).includes(String(n)),"never mentions this operand"),fn("hasOperation",(e,n)=>gn(e,String(n)),"uses this operation somewhere inside"),fn("equals",(e,n)=>e.s===String(n),"is exactly this postfix expression");const vn=e=>e.type===M?null:e.token;function bn(e){if("string"==typeof e)return{name:e,args:[]};if(Array.isArray(e)&&e.length>0&&"string"==typeof e[0])return{name:e[0],args:e.slice(1)};throw new cn(`bad requirement: ${JSON.stringify(e)}`)}pn("freeOf",(e,n)=>{const t=vn(n);return null!==t&&!L(e).includes(t)},"the first operand never mentions the second"),pn("dependsOn",(e,n)=>{const t=vn(n);return null!==t&&L(e).includes(t)},"the first operand does mention the second"),pn("sameAs",(...e)=>e.every(n=>n.s===e[0].s),"these operands are structurally identical"),pn("differs",(...e)=>new Set(e.map(e=>e.s)).size===e.length,"these operands are all different"),pn("sameValue",(...e)=>e.every(e=>mn(e))&&e.every(n=>n.value===e[0].value),"these operands are the same number");const yn="simplify";class wn extends Error{constructor(e,n){super(n?`rule "${n}": ${e}`:e),this.name="RuleError",this.ruleId=n}}const xn=/^[A-Za-z_][A-Za-z0-9_]*$/;function qn(e,n){if(F(e))return{literal:e.value};const t=e.token;if(t.startsWith("=")){const e=t.slice(1);if(""===e)throw new wn('"=" needs an operand after it',n);return{token:e}}if(!xn.test(t))throw new wn(`"${t}" is not a usable slot name; use letters, digits and _`,n);if(S(t))throw new wn(`slot name "${t}" collides with an operator`,n);return{name:t}}function kn(e){if(!V(e))return[e];const n=e.args.map(kn),t=$(e.op),r=[];for(const s of function(e){let n=[[]];for(const t of e){const e=[];for(const r of n)for(const n of t)e.push(r.concat([n]));if(n=e,n.length>64)return n.slice(0,64)}return n}(n))if(r.push(C(e.op,s)),t.commutative&&2===s.length&&s[0].s!==s[1].s&&r.push(C(e.op,[s[1],s[0]])),r.length>=64)break;return r}function $n(e,n=[]){return"op"===e.kind?e.args.forEach(e=>$n(e,n)):"const"===e.kind&&e.free&&n.push(e.node.token),n}function Sn(e){if(!e||"object"!=typeof e)throw new wn("rule must be a mapping");const n=e.id||e.name;if(!n)throw new wn("rule needs an id");if(!e.pattern)throw new wn("rule needs a pattern",n);if(void 0===e.output||null===e.output)throw new wn("rule needs an output",n);let t;try{t=R(String(e.pattern))}catch(e){if(e instanceof T)throw new wn(`bad pattern: ${e.message}`,n);throw e}const r=an(t).map(e=>qn(e,n)),s=new Set(r.filter(e=>e.name).map(e=>e.name)),o=function(e,n,t){const r=[],s=e.when||{};if("object"!=typeof s||Array.isArray(s))throw new wn('"when" must map slot names to lists of requirements',t);for(const[e,o]of Object.entries(s)){if(!n.has(e))throw new wn(`"when" mentions "${e}", which is not a slot in the pattern`,t);const s=Array.isArray(o)?o:[o];for(const n of s){const{name:s,args:o}=bn(n),i=dn(s);if(!i)throw new wn(`unknown requirement "${s}"`,t);r.push({slot:e,name:s,args:o,fn:i.fn})}}return r}(e,s,n),i=function(e,n,t){const r=[],s=e.require||[];if(!Array.isArray(s))throw new wn('"require" must be a list of relations',t);for(const e of s){const{name:s,args:o}=bn(e),i=hn(s);if(!i)throw new wn(`unknown relation "${s}"`,t);if(o.length<2)throw new wn(`relation "${s}" needs at least two operands`,t);for(const e of o)if("string"==typeof e&&!n.has(e)&&xn.test(e))throw new wn(`relation "${s}" mentions "${e}", which is not a slot in the pattern`,t);r.push({name:s,args:o,fn:i.fn})}return r}(e,s,n),a=function(e,n,t){let r;try{r=R(e)}catch(n){throw new wn(`bad output "${e}": ${n.message}`,t)}const s=e=>V(e)?{kind:"op",op:e.op,args:e.args.map(s)}:F(e)?{kind:"const",node:e}:n.has(e.token)?{kind:"slot",name:e.token}:e.token.startsWith("=")?{kind:"const",node:Z(e.token.slice(1))}:{kind:"const",node:e,free:!0};return s(r)}(String(e.output),s,n),u=e.commutative?kn(t):[t],l=[],c=new Set;for(const e of u){const t=on(e),r=an(e).map(e=>qn(e,n)),s=t+"\0"+r.map(e=>e.name??(void 0!==e.literal?`#${e.literal}`:`=${e.token}`)).join(",");c.has(s)||(c.add(s),l.push({key:t,slots:r}))}const f={id:String(n),name:e.name?String(e.name):String(n),pattern:String(e.pattern),patternTree:t,output:String(e.output),outputTemplate:a,explain:e.explain?String(e.explain).trim():"",mode:e.mode?String(e.mode):yn,tags:Array.isArray(e.tags)?e.tags.map(String):[],priority:Number.isFinite(e.priority)?e.priority:0,fold:!!e.fold,commutative:!!e.commutative,enabled:void 0===e.enabled||!!e.enabled,when:o,require:i,slots:[...s],depth:P(t),bindings:l,freeSymbols:$n(a),source:e};for(const e of l)e.rule=f;return f}function On(e){const n={};for(const t of e.when)(n[t.slot]||=[]).push(t.args.length?[t.name,...t.args]:t.name);const t={id:e.id};return e.name!==e.id&&(t.name=e.name),e.mode!==yn&&(t.mode=e.mode),t.pattern=e.pattern,Object.keys(n).length&&(t.when=n),e.require.length&&(t.require=e.require.map(e=>[e.name,...e.args])),t.output=e.output,e.explain&&(t.explain=e.explain),e.fold&&(t.fold=!0),e.commutative&&(t.commutative=!0),e.priority&&(t.priority=e.priority),e.tags.length&&(t.tags=e.tags),e.enabled||(t.enabled=!1),t}class Nn extends Error{constructor(e,n){super(null==n?e:`${e} (line ${n+1})`),this.name="YamlError",this.line=n}}class An{constructor(e){this.records=e,this.i=0,this.pushed=[]}peek(){if(this.pushed.length)return this.pushed[this.pushed.length-1];for(;this.i<this.records.length&&this.records[this.i].blank;)this.i++;return this.i<this.records.length?this.records[this.i]:null}next(){if(this.pushed.length)return this.pushed.pop();const e=this.peek();return e&&this.i++,e}push(e){this.pushed.push(e)}}const In=e=>"-"===e||e.startsWith("- ");function Mn(e){let n=null,t=0;for(let r=0;r<e.length;r++){const s=e[r];if(n)s===n?"'"===n&&"'"===e[r+1]?r++:n=null:'"'===n&&"\\"===s&&r++;else if('"'===s||"'"===s)n=s;else if("["===s||"{"===s)t++;else if("]"===s||"}"===s)t--;else if(":"===s&&0===t&&(r+1===e.length||" "===e[r+1])){const n=e.slice(0,r).trim();return""===n?null:{key:En(n),rest:e.slice(r+1).trim()}}}return null}function En(e){return e.length>=2&&'"'===e[0]&&'"'===e[e.length-1]?JSON.parse(e):e.length>=2&&"'"===e[0]&&"'"===e[e.length-1]?e.slice(1,-1).replace(/''/g,"'"):e}function Tn(e){const n=new An(function(e){return String(e).split(/\r?\n/).map((e,n)=>{const t=e.search(/\S/);if(-1===t)return{raw:e,n:n,indent:0,text:"",blank:!0};const r=function(e){let n=null;for(let t=0;t<e.length;t++){const r=e[t];if(n)r===n?"'"===n&&"'"===e[t+1]?t++:n=null:'"'===n&&"\\"===r&&t++;else if('"'===r||"'"===r)n=r;else if("#"===r&&(0===t||/\s/.test(e[t-1])))return e.slice(0,t)}return e}(e.slice(t)).trimEnd();return{raw:e,n:n,indent:t,text:r,blank:""===r}})}(e));let t=n.peek();if(t&&"---"===t.text&&(n.next(),t=n.peek()),!t)return null;const r=Zn(n,t.indent),s=n.peek();if(s&&"..."!==s.text&&"---"!==s.text)throw new Nn(`unexpected content "${s.text}"`,s.n);return r}const jn=/^[|>][-+]?\d*$/;function Zn(e,n){const t=e.peek();return!t||t.indent<n?null:In(t.text)?zn(e,t.indent):Mn(t.text)?function(e,n){const t={};for(;;){const r=e.peek();if(!r||r.indent!==n)break;const s=Mn(r.text);if(!s)break;e.next();const{key:o,rest:i}=s;if(""===i){const r=e.peek();r&&r.indent>n?t[o]=Zn(e,r.indent):r&&r.indent===n&&In(r.text)?t[o]=zn(e,n):t[o]=null}else jn.test(i)?t[o]=Dn(e,i,n):t[o]=Bn(i,r.n)}return t}(e,t.indent):(e.next(),jn.test(t.text)?Dn(e,t.text,t.indent-1):Bn(t.text,t.n))}function zn(e,n){const t=[];for(;;){const r=e.peek();if(!r||r.indent!==n||!In(r.text))break;if(e.next(),"-"===r.text){const r=e.peek();t.push(r&&r.indent>n?Zn(e,r.indent):null);continue}if(jn.test(r.text.slice(2).trim())){t.push(Dn(e,r.text.slice(2).trim(),n));continue}const s=r.text.slice(2),o=n+2+(s.length-s.trimStart().length);e.push({raw:r.raw,n:r.n,indent:o,text:s.trim(),blank:!1}),t.push(Zn(e,o))}return t}function Dn(e,n,t){const r=n[0],s=(n.match(/[-+]/)||[""])[0],o=(n.match(/\d+/)||[null])[0];let i=o?t+Number(o):null;const a=[];for(;e.i<e.records.length;){const n=e.records[e.i];if(""===n.raw.trim()){a.push(""),e.i++;continue}const r=n.raw.search(/\S/);if(r<=t)break;if(null===i&&(i=r),r<i)break;a.push(n.raw.slice(i)),e.i++}for(;a.length&&""===a[a.length-1];)a.pop();let u;if("|"===r)u=a.join("\n");else{const e=[];for(const n of a)""===n?e.push("\n"):e.length&&"\n"!==e[e.length-1]?e[e.length-1]+=" "+n:e.push(n);u=e.join("").replace(/\n(?!$)/g,"\n")}return"-"===s?u:"+"===s?u+"\n":""===u?"":u+"\n"}function Bn(e,n){const t=e.trim();return""===t||"~"===t||"null"===t||"Null"===t||"NULL"===t?null:"true"===t||"True"===t||"TRUE"===t||"false"!==t&&"False"!==t&&"FALSE"!==t&&("["===t[0]||"{"===t[0]?function(e,n){const t={text:e,i:0,line:n},r=Fn(t);if(Cn(t),t.i<t.text.length)throw new Nn(`unexpected "${t.text.slice(t.i)}" in flow value`,n);return r}(t,n):'"'===t[0]||"'"===t[0]?En(t):/^[-+]?(\d+\.?\d*|\.\d+)([eE][-+]?\d+)?$/.test(t)?Number(t):t)}const Cn=e=>{for(;e.i<e.text.length&&/\s/.test(e.text[e.i]);)e.i++};function Fn(e){Cn(e);const n=e.text[e.i];return"["===n?function(e){e.i++;const n=[];for(;;){if(Cn(e),e.i>=e.text.length)throw new Nn("unterminated [",e.line);if("]"===e.text[e.i])return e.i++,n;if(n.push(Fn(e)),Cn(e),","===e.text[e.i])e.i++;else if("]"!==e.text[e.i])throw new Nn(`expected , or ] near "${e.text.slice(e.i)}"`,e.line)}}(e):"{"===n?function(e){e.i++;const n={};for(;;){if(Cn(e),e.i>=e.text.length)throw new Nn("unterminated {",e.line);if("}"===e.text[e.i])return e.i++,n;const t=String(Fn(e));if(Cn(e),":"!==e.text[e.i])throw new Nn(`expected : after "${t}"`,e.line);if(e.i++,n[t]=Fn(e),Cn(e),","===e.text[e.i])e.i++;else if("}"!==e.text[e.i])throw new Nn(`expected , or } near "${e.text.slice(e.i)}"`,e.line)}}(e):Bn(function(e){const n=e.text[e.i];if('"'===n||"'"===n){const t=e.i;for(e.i++;e.i<e.text.length;)if("\\"===e.text[e.i]&&'"'===n)e.i+=2;else if(e.text[e.i]===n){if("'"!==n||"'"!==e.text[e.i+1])return e.i++,e.text.slice(t,e.i);e.i+=2}else e.i++;throw new Nn("unterminated quoted string",e.line)}const t=e.i;for(;e.i<e.text.length;){const n=e.text[e.i];if(","===n||"]"===n||"}"===n)break;if(":"===n){const n=e.text[e.i+1];if(void 0===n||/[\s,\]}]/.test(n))break}e.i++}return e.text.slice(t,e.i).trim()}(e),e.line)}const Vn=/^(true|false|null|yes|no|on|off|~)$/i;function Rn(e){if(null==e)return"null";if("boolean"==typeof e)return e?"true":"false";if("number"==typeof e){if(!Number.isFinite(e))throw new Nn(`cannot serialise ${e}`);return String(e)}const n=String(e);return function(e){return!!(""===e||/^\s|\s$/.test(e)||/[\u0000-\u001f\u007f]/.test(e)||/^[-?:,[\]{}#&*!|>'"%@`]/.test(e)||/: |\s#/.test(e)||Vn.test(e)||/^[-+]?(\d+\.?\d*|\.\d+)([eE][-+]?\d+)?$/.test(e))}(n)?JSON.stringify(n):n}const Un=e=>null==e||"object"!=typeof e;function Pn(e){return Array.isArray(e)?`[${e.map(Pn).join(", ")}]`:e&&"object"==typeof e?`{${Object.entries(e).map(([e,n])=>`${Rn(e)}: ${Pn(n)}`).join(", ")}}`:Rn(e)}function Wn(e){return Array.isArray(e)?!!e.every(e=>Un(e)||Array.isArray(e)&&e.every(Un))&&Pn(e).length<=72:!(!e||"object"!=typeof e)&&(Object.values(e).every(e=>Un(e)||Array.isArray(e)&&e.every(Un))&&Pn(e).length<=72)}function Ln(e){const n=[];return _n(e,0,n),n.length?n.join("\n")+"\n":""}function _n(e,n,t){const r=" ".repeat(n);if(Array.isArray(e)){if(0===e.length)return void t.push(r+"[]");for(const r of e)Hn(r,n,t)}else{if(e&&"object"==typeof e){const s=Object.entries(e).filter(([,e])=>void 0!==e);if(0===s.length)return void t.push(r+"{}");for(const[e,r]of s)Jn(e,r,n,t);return}!function(e,n,t){if("string"==typeof e&&e.includes("\n")){t.push(n+"|-");for(const r of e.replace(/\n+$/,"").split("\n"))t.push(n+"  "+r);return}t.push(n+Rn(e))}(e,r,t)}}function Hn(e,n,t){const r=" ".repeat(n);if(null!==e&&"object"==typeof e&&!Wn(e)){const s=[];if(_n(e,n+2,s),0===s.length)return void t.push(r+"-");t.push(r+"- "+s[0].slice(n+2));for(let e=1;e<s.length;e++)t.push(s[e]);return}if(null===e||"object"!=typeof e)if("string"==typeof e&&e.includes("\n")){t.push(r+"- |-");for(const r of e.replace(/\n+$/,"").split("\n"))t.push(" ".repeat(n+4)+r)}else t.push(r+"- "+Rn(e));else t.push(r+"- "+Pn(e))}function Jn(e,n,t,r){const s=" ".repeat(t),o=Rn(e);if("string"==typeof n&&n.includes("\n")){r.push(`${s}${o}: |-`);for(const e of n.replace(/\n+$/,"").split("\n"))r.push(" ".repeat(t+2)+e)}else{if(null!==n&&"object"==typeof n)return Wn(n)?void r.push(`${s}${o}: ${Pn(n)}`):(r.push(`${s}${o}:`),void _n(n,Array.isArray(n)?t:t+2,r));r.push(`${s}${o}: ${Rn(n)}`)}}class Qn{constructor(e=[]){this.index=new Map,this.rules=[],this.maxDepth=0,this.name="",this.views=new Map;for(const n of e)this.add(n)}add(e){const n=e&&e.bindings?e:Sn(e);if(this.rules.some(e=>e.id===n.id))throw new wn(`duplicate rule id "${n.id}"`);return this.rules.push(n),n.enabled&&this.link(n),n.depth>this.maxDepth&&(this.maxDepth=n.depth),this.views.clear(),n}view(e){if(!e)return this;let n=this.views.get(e);return n||(n=new Qn(this.rules.filter(n=>n.mode===e)),n.name=this.name,this.views.set(e,n)),n}modes(){return[...new Set(this.rules.map(e=>e.mode))].sort((e,n)=>e===yn?-1:n===yn?1:e.localeCompare(n))}link(e){for(const n of e.bindings){let e=this.index.get(n.key);e||this.index.set(n.key,e=[]),e.push(n),e.sort((e,n)=>n.rule.priority-e.rule.priority||n.rule.depth-e.rule.depth)}}remove(e){const n=this.rules.findIndex(n=>n.id===e);return-1!==n&&(this.rules.splice(n,1),this.reindex(),!0)}setEnabled(e,n){const t=this.get(e);return!!t&&(t.enabled=!!n,this.reindex(),!0)}reindex(){this.index=new Map,this.maxDepth=0,this.views.clear();for(const e of this.rules)e.enabled&&this.link(e),e.depth>this.maxDepth&&(this.maxDepth=e.depth)}get(e){return this.rules.find(n=>n.id===e)||null}candidates(e){return this.index.get(e)||null}get size(){return this.rules.length}static fromYaml(e,{strict:n=!0}={}){const t=Tn(e);return Qn.fromSpecs(t,{strict:n})}static fromSpecs(e,{strict:n=!0}={}){let t=e,r="";if(e&&!Array.isArray(e)&&"object"==typeof e&&(t=e.rules,r=e.name?String(e.name):""),!Array.isArray(t))throw new wn('expected a list of rules, or a mapping with a "rules" list');const s=new Qn;s.name=r;const o=[];for(let e=0;e<t.length;e++)try{s.add(t[e])}catch(n){o.push(`#${e+1}: ${n.message}`)}if(o.length&&n)throw new wn(`${o.length} rule(s) failed to compile:\n  ${o.join("\n  ")}`);return s.errors=o,s}toSpecs(){return this.rules.map(On)}toYaml(){const e=this.toSpecs();return this.name?Ln({name:this.name,rules:e}):Ln(e)}keyReport(){return[...this.index.entries()].map(([e,n])=>({key:e,rules:n.map(e=>e.rule.id)})).sort((e,n)=>e.key.localeCompare(n.key))}}function Yn(e,n){const t=Object.create(null),r=e.slots;for(let e=0;e<r.length;e++){const s=r[e],o=n[e];if(void 0!==s.literal){if(!F(o)||o.value!==s.literal)return null;continue}if(void 0!==s.token){if(o.s!==s.token)return null;continue}const i=t[s.name];if(void 0!==i){if(i.s!==o.s)return null}else t[s.name]=o}return t}function Gn(e,n){for(const t of e.when){const e=n[t.slot];if(void 0===e)return!1;if(!t.fn(e,...t.args))return!1}for(const t of e.require){const e=t.args.map(e=>void 0!==n[e]?n[e]:e);if(!t.fn(...e))return!1}return!0}function Kn(e,n){return n.rule.priority-e.rule.priority||n.rule.depth-e.rule.depth}function Xn(e,n){const t=[];if(0===n.maxDepth)return t;for(const{key:r,slots:s}of function(e,n){const t=sn(e,n),r=[];for(const e of t)e.key!==nn&&r.push(e);return r}(e,n.maxDepth)){const e=n.candidates(r);if(e)for(const n of e){const e=Yn(n,s);e&&(Gn(n.rule,e)&&t.push({rule:n.rule,binding:n,key:r,bound:e}))}}return t.sort(Kn)}function et(e,n){switch(e.kind){case"slot":{const t=n[e.name];if(!t)throw new Error(`output references unbound slot "${e.name}"`);return t}case"const":return e.node;case"op":return C(e.op,e.args.map(e=>et(e,n)));default:throw new Error(`bad output template: ${JSON.stringify(e)}`)}}const nt=4e5,tt=()=>new Q("a number too large to approximate"),rt=e=>e<0n?-e:e;class st{constructor(e){this.P=e,this.S=10n**BigInt(e),this.cache=new Map}fromQ(e){return e.n*this.S/e.d}mul(e,n){return e*n/this.S}div(e,n){if(0n===n)throw new Q("division by zero");return e*this.S/n}sqrt(e){if(e<0n)throw new Q("the square root of a negative number is not real");return te(e*this.S,2)}series(e,n){let t=0n,r=e;for(let e=1;0n!==r;e++)if(t+=r,r=n(r,e),e>1e5)throw new Q("a series did not converge");return t}atanhInv(e){const n=BigInt(e),t=n*n;let r=this.S/n,s=0n;for(let e=1n;0n!==r;e+=2n)s+=r/e,r/=t;return s}atanInv(e){const n=BigInt(e),t=n*n;let r=this.S/n,s=0n;for(let e=1n,n=1n;0n!==r;e+=2n,n=-n)s+=n*r/e,r/=t;return s}get pi(){return this.memo("pi",()=>16n*this.atanInv(5)-4n*this.atanInv(239))}get ln2(){return this.memo("ln2",()=>2n*this.atanhInv(3))}get e(){return this.memo("e",()=>this.exp(this.S))}memo(e,n){return this.cache.has(e)||this.cache.set(e,n()),this.cache.get(e)}exp(e){const n=(e+(e>=0n?this.ln2/2n:-this.ln2/2n))/this.ln2;if(n>BigInt(nt))throw tt();if(n<-BigInt(nt))return 0n;const t=e-n*this.ln2,r=this.series(this.S,(e,n)=>e*t/this.S/BigInt(n));return n>=0n?r<<n:r>>-n}ln(e){if(e<=0n)throw new Q("the log of a number that is not positive is not real");let n=0n,t=e;for(;t>=this.S;)t>>=1n,n++;for(;t<this.S/2n;)t<<=1n,n--;const r=this.div(t-this.S,t+this.S),s=this.mul(r,r);let o=r,i=0n;for(let e=1n;0n!==o;e+=2n)i+=o/e,o=this.mul(o,s);return 2n*i+n*this.ln2}sin(e){const n=2n*this.pi;let t=e%n;t>this.pi&&(t-=n),t<-this.pi&&(t+=n);const r=this.mul(t,t);return this.series(t,(e,n)=>-e*r/this.S/BigInt(2*n*(2*n+1)))}cos(e){return this.sin(e+this.pi/2n)}atan(e){if(rt(e)>this.S){const n=this.pi/2n;return(e>0n?n:-n)-this.atan(this.div(this.S,e))}let n=e;for(let e=0;e<2;e++)n=this.div(n,this.S+this.sqrt(this.S+this.mul(n,n)));const t=this.mul(n,n);let r=n,s=0n;for(let e=1n,n=1n;0n!==r;e+=2n,n=-n)s+=n*r/e,r=this.mul(r,t);return 4n*s}asin(e){if(rt(e)>this.S)throw new Q("arcsin of a number beyond ±1 is not real");return rt(e)===this.S?e>0n?this.pi/2n:-this.pi/2n:this.atan(this.div(e,this.sqrt(this.S-this.mul(e,e))))}}function ot(e,n,t){if(t>0n)return e.atan(e.div(n,t));if(t<0n)return e.atan(e.div(n,t))+(n>=0n?e.pi:-e.pi);if(0n===n)throw new Q("the argument of 0 is undefined");return n>0n?e.pi/2n:-e.pi/2n}const it={add:(e,n)=>[e[0]+n[0],e[1]+n[1]],sub:(e,n)=>[e[0]-n[0],e[1]-n[1]],mul:(e,n,t)=>[e.mul(n[0],t[0])-e.mul(n[1],t[1]),e.mul(n[0],t[1])+e.mul(n[1],t[0])],div(e,n,t){const r=e.mul(t[0],t[0])+e.mul(t[1],t[1]);if(0n===r)throw new Q("division by zero");return[e.div(e.mul(n[0],t[0])+e.mul(n[1],t[1]),r),e.div(e.mul(n[1],t[0])-e.mul(n[0],t[1]),r)]},abs:(e,n)=>e.sqrt(e.mul(n[0],n[0])+e.mul(n[1],n[1])),exp(e,n){const t=e.exp(n[0]);return[e.mul(t,e.cos(n[1])),e.mul(t,e.sin(n[1]))]},ln(e,n){const t=it.abs(e,n);if(0n===t)throw new Q("the log of 0 is undefined");return[e.ln(t),ot(e,n[1],n[0])]},pow(e,n,t,r){if(0n===n[0]&&0n===n[1])return[0n,0n];if(r&&r.isInteger){let t=r.n<0n?-r.n:r.n,s=[e.S,0n],o=n;for(;t>0n;)1n&t&&(s=it.mul(e,s,o)),o=it.mul(e,o,o),t>>=1n;return r.n<0n?it.div(e,[e.S,0n],s):s}return it.exp(e,it.mul(e,t,it.ln(e,n)))},sinh:(e,n)=>(e.exp(n)-e.exp(-n))/2n,cosh:(e,n)=>(e.exp(n)+e.exp(-n))/2n},at={sin:(e,[n,t])=>[e.mul(e.sin(n),it.cosh(e,t)),e.mul(e.cos(n),it.sinh(e,t))],cos:(e,[n,t])=>[e.mul(e.cos(n),it.cosh(e,t)),-e.mul(e.sin(n),it.sinh(e,t))],tan:(e,n)=>it.div(e,at.sin(e,n),at.cos(e,n)),ln:(e,n)=>it.ln(e,n),sqrt:(e,n)=>it.exp(e,it.mul(e,[e.S/2n,0n],it.ln(e,n))),sinh:(e,[n,t])=>[e.mul(it.sinh(e,n),e.cos(t)),e.mul(it.cosh(e,n),e.sin(t))],cosh:(e,[n,t])=>[e.mul(it.cosh(e,n),e.cos(t)),e.mul(it.sinh(e,n),e.sin(t))]},ut=e=>({q:e}),lt=([e,n])=>0n===n?{f:e}:{z:[e,n]};function ct(e,n,t){const r=e=>((e,n)=>n.z??[ft(e,n),0n])(t,n[e]);switch(e.op){case"+":return lt(it.add(r(0),r(1)));case"-":return lt(it.sub(r(0),r(1)));case"neg":return lt([-r(0)[0],-r(0)[1]]);case"*":return lt(it.mul(t,r(0),r(1)));case"/":return lt(it.div(t,r(0),r(1)));case"^":return lt(it.pow(t,r(0),r(1),n[1].q));case"sqrt":return lt(at.sqrt(t,r(0)));case"root":return lt(it.pow(t,r(1),it.div(t,[t.S,0n],r(0)),n[0].q?.inv()));case"abs":return{f:it.abs(t,r(0))};case"arg":return{f:ot(t,r(0)[1],r(0)[0])};case"re":return{f:r(0)[0]};case"im":return{f:r(0)[1]};case"conj":return lt([r(0)[0],-r(0)[1]]);case"log":return lt(it.div(t,it.ln(t,r(1)),it.ln(t,r(0))));default:{const s=at[e.op];if(!s||1!==n.length)throw new Q(`no complex value for ${e.op}`);return lt(s(t,r(0)))}}}function ft(e,n){return n.q?e.fromQ(n.q):n.f}const pt={sin:(e,n)=>e.sin(n),cos:(e,n)=>e.cos(n),tan:(e,n)=>e.div(e.sin(n),e.cos(n)),sec:(e,n)=>e.div(e.S,e.cos(n)),csc:(e,n)=>e.div(e.S,e.sin(n)),cot:(e,n)=>e.div(e.cos(n),e.sin(n)),arcsin:(e,n)=>e.asin(n),arccos:(e,n)=>e.pi/2n-e.asin(n),arctan:(e,n)=>e.atan(n),sinh:(e,n)=>(e.exp(n)-e.exp(-n))/2n,cosh:(e,n)=>(e.exp(n)+e.exp(-n))/2n,tanh:(e,n)=>{const t=e.exp(n),r=e.exp(-n);return e.div(t-r,t+r)},arcsinh:(e,n)=>e.ln(n+e.sqrt(e.mul(n,n)+e.S)),arccosh:(e,n)=>{if(n<e.S)throw new Q("arccosh of a number below 1 is not real");return e.ln(n+e.sqrt(e.mul(n,n)-e.S))},arctanh:(e,n)=>{if(rt(n)>=e.S)throw new Q("arctanh of a number beyond ±1 is not real");return e.ln(e.div(e.S+n,e.S-n))/2n},ln:(e,n)=>e.ln(n),loglog:(e,n)=>e.ln(e.ln(n)),logloglog:(e,n)=>e.ln(e.ln(e.ln(n))),expm1:(e,n)=>e.exp(n)-e.S};function dt(e,n,t){if("num"===e.type)return ut(ne(e.token));if("sym"===e.type){if("pi"===e.token||"π"===e.token)return{f:n.pi};if("e"===e.token)return{f:n.e};if("i"===e.token&&ie.complex)return{z:[0n,n.S]};if(t&&e.token in t){const n=t[e.token],r="number"==typeof n||"string"==typeof n||"bigint"==typeof n?ne(String(n)):null;if(!r)throw new Q(`the value of ${e.token} is not a number`);return ut(r)}throw new Q(`${e.token} has no value`)}ie.complex&&e.op.startsWith("\\:")&&Be[e.op.slice(2)]&&(e=C(Be[e.op.slice(2)],e.args));const r=e.args.map(e=>dt(e,n,t));if(r.some(e=>e.z)||["re","im","conj","arg"].includes(e.op)||function(e,n){if(!ie.complex)return!1;const t=e=>e.q?e.q.sign<0:void 0!==e.f&&e.f<0n;return"sqrt"===e.op||"ln"===e.op||"log"===e.op?t(n["log"===e.op?1:0])||"log"===e.op&&t(n[0]):"^"===e.op?t(n[0])&&!(n[1].q&&n[1].q.d%2n==1n):"root"===e.op&&t(n[1])&&!(n[0].q&&n[0].q.isInteger&&n[0].q.n%2n==1n)}(e,r))return ct(e,r,n);const s=r.every(e=>e.q),o=e=>ft(n,r[e]);switch(e.op){case"+":return s?ut(r[0].q.add(r[1].q)):{f:o(0)+o(1)};case"-":return s?ut(r[0].q.sub(r[1].q)):{f:o(0)-o(1)};case"neg":return r[0].q?ut(r[0].q.neg()):{f:-r[0].f};case"*":return s?ut(r[0].q.mul(r[1].q)):{f:n.mul(o(0),o(1))};case"/":return s?ut(r[0].q.div(r[1].q)):{f:n.div(o(0),o(1))};case"^":return mt(n,r[0],r[1]);case"sqrt":return mt(n,r[0],ut(ee.HALF));case"root":return mt(n,r[1],r[0].q?ut(r[0].q.inv()):{f:n.div(n.S,r[0].f)});case"abs":return r[0].q?ut(r[0].q.abs()):{f:rt(r[0].f)};case"floor":return ut(r[0].q?r[0].q.floor():new ee(ht(r[0].f,n.S)));case"ceil":return ut(r[0].q?r[0].q.ceil():new ee(-ht(-r[0].f,n.S)));case"log":return{f:n.div(n.ln(o(1)),n.ln(o(0)))};case"atan2":return{f:ot(n,o(0),o(1))};case"max":case"min":if(!r.length)throw new Q(`${e.op} of nothing`);return s?ut(r.map(e=>e.q).reduce((n,t)=>("max"===e.op?t.sub(n).sign>0:t.sub(n).sign<0)?t:n)):{f:r.map((e,n)=>o(n)).reduce((n,t)=>"max"===e.op?t>n?t:n:t<n?t:n)};case"sum":return s?ut(r.reduce((e,n)=>e.add(n.q),ee.ZERO)):{f:r.reduce((e,n,t)=>e+o(t),0n)};case"prod":return s?ut(r.reduce((e,n)=>e.mul(n.q),ee.ONE)):{f:r.reduce((e,t,r)=>n.mul(e,o(r)),n.S)};case"mean":if(!r.length)throw new Q("mean of nothing");return s?ut(r.reduce((e,n)=>e.add(n.q),ee.ZERO).div(new ee(BigInt(r.length)))):{f:r.reduce((e,n,t)=>e+o(t),0n)/BigInt(r.length)};case"mod":{if(s){if(r[1].q.isZero)throw new Q("mod 0 is undefined");return ut(r[0].q.sub(r[1].q.mul(r[0].q.div(r[1].q).floor())))}const e=o(1);if(0n===e)throw new Q("mod 0 is undefined");return{f:o(0)-n.mul(e,ht(n.div(o(0),e),n.S)*n.S)}}case"round":{if(!r[1].q?.isInteger)throw new Q("round needs a whole number of places");const e=10n**rt(r[1].q.n),t=r[0].q??new ee(o(0),n.S),s=r[1].q.sign>=0,i=(s?t.mul(new ee(e)):t.div(new ee(e))).add(ee.HALF).floor();return ut(s?i.div(new ee(e)):i.mul(new ee(e)))}case"!":{if(!r[0].q?.isInteger||r[0].q.sign<0)throw new Q("factorial is only worked out for whole numbers");let e=1n;for(let n=2n;n<=r[0].q.n;n++)if(e*=n,0n==(255n&n)&&K(e)>nt)throw tt();return ut(new ee(e))}default:{const t=pt[e.op];if(!t||1!==r.length)throw new Q(`no numeric value for ${e.op}`);return{f:t(n,o(0))}}}}const ht=(e,n)=>{const t=e/n;return e<0n&&t*n!==e?t-1n:t};function mt(e,n,t){if(t.q&&t.q.isInteger){if(n.q)return ut(n.q.pow(t.q.n));let r=t.q.n;const s=r<0n;s&&(r=-r);let o=e.S,i=n.f;for(;r>0n;){if(1n&r&&(o=e.mul(o,i)),i=e.mul(i,i),K(o)>nt+4*e.P||K(i)>8e5){if(s)return{f:0n};throw tt()}r>>=1n}return{f:s?e.div(e.S,o):o}}const r=ft(e,n);if(0n===r){if(t.q?t.q.sign>0:t.f>0n)return ut(ee.ZERO);throw new Q("0 to a power that is not positive is undefined")}if(r<0n){if(t.q&&t.q.d%2n==1n){const n=mt(e,{f:-r},t).f;return{f:t.q.n%2n==0n?n:-n}}throw new Q("a negative number to a fractional power is not real")}if(t.q&&1n===rt(t.q.n)&&t.q.d<=64n){const n=te(r*e.S**(t.q.d-1n),Number(t.q.d));return{f:t.q.n<0n?e.div(e.S,n):n}}return{f:e.exp(e.mul(ft(e,t),e.ln(r)))}}function gt(e,n,t){const r=e<0n?"-":"",s=rt(e);if(0n===s)return{sign:"",digits:"0",point:1};const o=s.toString(),i=o.length-n.P;let a=o.slice(0,t);o.length>t&&o[t]>="5"&&(a=(BigInt(a)+1n).toString());let u=i;return a.length>t&&(a=a.slice(0,t),u++),{sign:r,digits:a.replace(/0+$/,"")||"0",point:u}}function vt({sign:e,digits:n,point:t}){if("0"===n)return"0";if(t>21||t<-6){return`${e}${n.length>1?`${n[0]}.${n.slice(1)}`:n}e${t-1}`}return t<=0?`${e}0.${"0".repeat(-t)}${n}`:t>=n.length?`${e}${n}${"0".repeat(t-n.length)}`:`${e}${n.slice(0,t)}.${n.slice(t)}`}function bt(e){let n=0n,t=0n,r=e.d;for(;r%2n==0n;)r/=2n,n++;for(;r%5n==0n;)r/=5n,t++;const s=n>t?n:t;let o=rt(e.n)*10n**s/e.d;if(0n===o)return 1;for(;o%10n==0n;)o/=10n;return o.toString().length}function yt(e){let n=e.d;for(;n%2n==0n;)n/=2n;for(;n%5n==0n;)n/=5n;return 1n===n}function wt(e){const n=/^(-?[\d.]+)e(-?\d+)$/.exec(e);return n?C("*",[z(n[1]),C("^",[z("10"),z(n[2])])]):z(e)}function xt(e,{digits:n=15,scope:t}={}){if(!Number.isInteger(n)||n<1||n>1e5)throw new Q("digits must be a whole number from 1 to 100000");let r=n+25,s=null;for(let o=0;o<4;o++){const i=new st(r),a=dt(e,i,t);if(a.z){const e=kt(a.z,i,n);if(e.retry&&o<3){r+=e.retry;continue}return void 0===e.real?e:qt(e.real,i,n,null)}if(a.q&&a.q.isInteger&&rt(a.q.n).toString().length<=n)return{text:a.q.n.toString(),node:z(a.q.n.toString()),number:Number(a.q.n),re:Number(a.q.n),im:0,exact:!0};const u=ft(i,a);if(0n===u&&(!a.q||!a.q.isZero)&&o<3){s=null,r*=4;continue}const l=gt(u,i,n);if(s&&l.point<-10&&s.point-l.point>=(r-s.P)/2)return{text:"0",node:z("0"),number:0,re:0,im:0,exact:!1};if(!(l.point<-10&&o<3)||s&&s.point===l.point)return qt(u,i,n,a.q);s={point:l.point,P:r},r+=10-l.point}throw new Q("could not approximate")}function qt(e,n,t,r){const s=vt(gt(e,n,t));return{text:s,node:wt(s),number:Number(s),re:Number(s),im:0,exact:Boolean(r)&&yt(r)&&bt(r)<=t}}function kt([e,n],t,r){const s=10n**BigInt(Math.max(0,t.P-r-12)),o=e=>rt(e)>s;if(!o(n))return{real:e};const i=gt(n,t,r);if(i.point<-10&&o(n))return{retry:10-i.point};const a=vt(i),u="1"===a?"i":"-1"===a?"-i":`${a}i`,l="1"===a?Z("i"):C("*",[wt(a.replace(/^-/,"")),Z("i")]);if(!o(e))return{text:u,node:a.startsWith("-")?C("neg",[l]):l,number:NaN,re:0,im:Number(a),exact:!1};const c=vt(gt(e,t,r)),f=a.startsWith("-"),p="1"===a.replace(/^-/,"")?Z("i"):C("*",[wt(a.replace(/^-/,"")),Z("i")]);return{text:`${c} ${f?"-":"+"} ${u.replace(/^-/,"")}`,node:C(f?"-":"+",[wt(c),p]),number:NaN,re:Number(c),im:Number(a),exact:!1}}function $t(e,n={}){const t=e=>(e=>0===Ke(e).size)(e)?"num"===e.type&&/^-?\d+$/.test(e.token)?e:xt(e,n).node:V(e)?C(e.op,e.args.map(t)):e;return t(e)}const St={variable:null},Ot=(e,n)=>null!==n&&Ke(e).has(n);function Nt(e){if("num"===e.type)return ne(e.token);if(V(e)&&"/"===e.op&&e.args.every(e=>"num"===e.type)){const n=ne(e.args[0].token),t=ne(e.args[1].token);return n&&t&&!t.isZero?n.div(t):null}return null}const At=e=>0===Ke(e).size;function It(e,n,t){const r=function(e){const n=We(e);return"num"===n.k?n.q.isZero?n.q:null:"atom"===n.k&&"sym"===n.node.type&&"pi"===n.node.token?ne("1"):"mul"===n.k&&1===n.fs.length&&"atom"===n.fs[0].b.k&&"pi"===n.fs[0].b.node.token&&"num"===n.fs[0].e.k&&n.fs[0].e.q.isOne?n.c:null}(e);if(!r)return!1;const[s,o="1"]=String(n).split("/");return t(r.compare(ne(s).div(ne(o))))}function Mt(e){if(!At(e))return NaN;try{const{number:n}=xt(e,{digits:30});return Math.sign(n)}catch{return NaN}}function Et(){dn("hasVariable")||(fn("hasVariable",e=>Ot(e,St.variable),"mentions the variable being solved for"),fn("freeOfVariable",e=>null!==St.variable&&!Ot(e,St.variable),"does not mention the variable being solved for"),fn("isRational",e=>null!==Nt(e),"an exact number: 3, -2, 1/2"),fn("isWhole",e=>Nt(e)?.isInteger??!1,"a whole number (any size)"),fn("isEvenWhole",e=>{const n=Nt(e);return Boolean(n?.isInteger&&n.n%2n==0n)},"an even whole number"),fn("isOddWhole",e=>{const n=Nt(e);return Boolean(n?.isInteger&&n.n%2n!=0n)},"an odd whole number"),fn("isPositiveRational",e=>(Nt(e)?.sign??0)>0,"an exact number above 0"),fn("isNegativeRational",e=>(Nt(e)?.sign??0)<0,"an exact number below 0"),fn("isConstant",At,"has no unknowns (numbers, π and e only)"),fn("notOne",e=>{const n=Nt(e);return!n||!n.isOne},"anything except the number 1"),fn("notMinusOne",e=>{const n=Nt(e);return!n||!n.equals(ne("-1"))},"anything except the number -1"),fn("hasMinusSign",e=>function(e){const n=We(e);return"num"===n.k?n.q.sign<0:"mul"===n.k&&n.c.sign<0}(e),"written with a leading minus: -x, -2x, -π/6"),fn("piMultipleAtLeast",(e,n)=>It(e,n,e=>e>=0),"a multiple of π, at least k·π"),fn("piMultipleAbove",(e,n)=>It(e,n,e=>e>0),"a multiple of π, above k·π"),fn("piMultipleBelow",(e,n)=>It(e,n,e=>e<0),"a multiple of π, below k·π"),fn("piMultipleAtMost",(e,n)=>It(e,n,e=>e<=0),"a multiple of π, at most k·π"),fn("isPositiveConstant",e=>Mt(e)>0,"no unknowns, and above 0 (e, π, √2, 3)"),fn("isNegativeConstant",e=>Mt(e)<0,"no unknowns, and below 0"),fn("hasImaginaryPart",e=>{if(!ie.complex)return!1;const n=Ze(We(e));return Boolean(n)&&!("num"===n.im.k&&n.im.q.isZero)},"complex mode: has a nonzero imaginary part (unknowns count as real)"),fn("isSymbolNode",e=>"sym"===e.type&&!oe.has(e.token),"a single unknown such as x"))}function Tt(e){const n=new Map;for(const t of e){const e=t.id??t.name;if(!e)throw new Error(`a rule needs an id (pattern ${t.pattern})`);n.set(String(e),{...n.get(String(e))??{},...t})}const t=[...n.values()].map(e=>{const n={...e};try{void 0!==n.pattern&&(n.pattern=function(e){const n=R(e),t=e=>"num"===e.type||("sym"===e.type?e.token.startsWith("="):e.args.every(t)),r=e=>"sym"===e.type?Z(e.token.slice(1)):"num"===e.type?e:C(e.op,e.args.map(r)),s=e=>"sym"===e.type?Z(`=${e.token}`):"num"===e.type?e:C(e.op,e.args.map(s)),o=e=>"op"!==e.type?e:t(e)?s(Le(r(e))):C(e.op,e.args.map(o));return o(n).s}(String(n.pattern)))}catch(e){throw new Error(`rule "${n.id}": bad pattern "${n.pattern}": ${e.message}`)}return"hidden"===n.step&&(n.tags=[...n.tags??[],"hidden"]),n});return Qn.fromSpecs(t,{strict:!0})}function jt(e,n,t,r){return e.explain?e.explain.replace(/\{(\w+)\}/g,(e,s)=>"result"===s?t?r(t):e:n[s]?r(n[s]):e).trim():""}const Zt=e=>V(e)&&("+"===e.op||"-"===e.op),zt=e=>V(e)&&("*"===e.op||"/"===e.op);function*Dt(e,n,{order:t="pre"}={}){for(const r of function(e,n,{order:t="post"}={}){const r=[];for(const{node:s,path:o}of H(e,t))for(const e of Xn(s,n))r.push({...e,node:s,path:o});return r}(e,n,{order:t})){const n=et(r.rule.outputTemplate,r.bound);yield{rule:r.rule,bound:r.bound,path:r.path,before:r.node,after:n,next:()=>Le(J(e,r.path,n))}}for(const{node:r,path:s}of H(e,t)){const t=Zt(r);if(!t&&!zt(r))continue;const o=t?Xe(We(r)):en(We(r));if(o.length<3)continue;const i=t?Se:Ne;for(let t=0;t<o.length;t++)for(let r=t+1;r<o.length;r++){const a=Qe(i([o[t],o[r]]));for(const u of Xn(a,n)){const n=et(u.rule.outputTemplate,u.bound),l=o.filter((e,n)=>n!==t&&n!==r);yield{rule:u.rule,bound:u.bound,path:s,before:a,after:n,next:()=>Le(J(e,s,Qe(i([We(n),...l]))))}}}}}function Bt(e,n,{format:t,order:r="pre",maxSteps:s=200,variable:o=null}={}){const i=St.variable;St.variable=o;try{let o=Le(e);const i=new Set([o.s]),a=[];for(let e=0;e<s;e++){let e=null;for(const t of Dt(o,n,{order:r})){const n=t.next();if(n.s!==o.s&&!i.has(n.s)){e={c:t,next:n};break}}if(!e)break;const{c:s,next:u}=e;a.push({id:s.rule.id,name:s.rule.name,tags:s.rule.tags,from:o.s,to:u.s,path:s.path,before:s.before.s,after:s.after.s,explanation:jt(s.rule,s.bound,Le(s.after),t)}),o=u,i.add(o.s)}return{node:o,steps:a}}finally{St.variable=i}}const Ct={full:[],normal:["hidden"],brief:["hidden","tidy","identity","value","split","base","arithmetic"]},Ft=(e,n)=>"all"===n||"*"===n||Array.isArray(n)&&n.some(n=>n===e.id||e.tags?.includes(n));function Vt(e){if(!1===e)return null;if(void 0===e||!0===e)return{show:"all",hide:Ct.normal,explain:!0};if("string"==typeof e){if(!Ct[e])throw new Error(`steps: "${e}" is not one of ${Object.keys(Ct).join(", ")}`);return{show:"all",hide:Ct[e],explain:!0}}const n=e.detail??"normal";if(!Ct[n])throw new Error(`steps.detail: "${n}" is not one of ${Object.keys(Ct).join(", ")}`);return{show:e.show??"all",hide:[...Ct[n],...e.hide??[]],explain:e.explain??!0}}function Rt(e,n){if(!n)return[];const t=[];let r=null;for(const s of e){if(!(Ft(s,n.show)&&!Ft(s,n.hide))){r=r??{from:s.from,ids:[]},r.ids.push(s.id);continue}const e={...s};r&&(e.from=r.from,e.also=r.ids,e.path=null,r=null);!0===n.explain||!1!==n.explain&&Ft(s,n.explain)||delete e.explanation,t.push(e)}if(r&&t.length){const n=t[t.length-1];n.to=e[e.length-1].to,n.also=[...n.also??[],...r.ids]}return t}const Ut=ce(ee.ONE),Pt=(e,n)=>"atom"===e.k&&"sym"===e.node.type&&e.node.token===n;function Wt(e){switch(e.k){case"num":default:return e;case"atom":return"op"===e.node.type?Ue(e.node.op,e.node.args.map(e=>Wt(We(e)))):e;case"add":return Se(Xe(e).map(Wt));case"mul":{let n=[ce(e.c)];for(const t of e.fs){const e=Wt(t.b);if("add"===e.k&&he(t.e)&&t.e.q.isInteger&&t.e.q.sign>0&&t.e.q.n<=BigInt(60))for(let r=0n;r<t.e.q.n;r++)n=Lt(n,Xe(e));else n=Lt(n,[Ie(e,t.e)])}return Se(n)}}}function Lt(e,n){const t=[];for(const r of e)for(const e of n)t.push(Ne([r,e]));return Xe(Se(t))}function _t(e,n){const t=new Map;for(const r of Xe(Wt(e))){const{coef:e,rest:s}=$e(r);let o=0;const i=[ce(e)];for(const e of s?en(s):[]){const{b:t,e:r}="mul"===e.k&&1===e.fs.length&&e.c.isOne?e.fs[0]:{b:e,e:Ut};if(Pt(t,n)){if(!he(r)||!r.q.isInteger||r.q.sign<0)return null;o+=Number(r.q.n)}else{if(Ht(t,n)||Ht(r,n))return null;i.push(e)}}t.set(o,Se([t.get(o)??ce(ee.ZERO),Ne(i)]))}for(const[e,n]of t)he(n)&&n.q.isZero&&t.delete(e);return t}const Ht=(e,n)=>Ke(Qe(e)).has(n);function Jt(e){if(!e)return null;const n=Math.max(0,...e.keys());if(n>256)return null;const t=[];for(let r=0;r<=n;r++){const n=e.get(r)??ce(ee.ZERO);if(!he(n))return null;t.push(n.q)}return t}function Qt(e,n){const t=We(Z(n));return Se(e.map((e,n)=>Ne([ce(e),Ae(t,ce(ee.of(n)))])))}const Yt=e=>{const n=e.slice();for(;n.length>1&&n[n.length-1].isZero;)n.pop();return n},Gt=(e,n)=>e.reduceRight((e,t)=>e.mul(n).add(t),ee.ZERO);function Kt(e,n){const t=e.length-1,r=new Array(t);let s=ee.ZERO;for(let o=t;o>=1;o--)s=e[o].add(s.mul(n)),r[o-1]=s;return r}function Xt(e){if(0n===(e=(e=>e<0n?-e:e)(e)))return null;if(e>10n**14n)return null;const n=[];let t=e;for(let e=2n;e*e<=t;e+=2n===e?1n:2n){if(t%e)continue;let r=0;for(;t%e===0n;)t/=e,r++;n.push([e,r])}t>1n&&n.push([t,1]);let r=[1n];for(const[e,t]of n){const n=[];for(const s of r)for(let r=0,o=1n;r<=t;r++,o*=e)n.push(s*o);if(r=n,r.length>2e4)return null}return r.sort((e,n)=>e<n?-1:e>n?1:0)}function er(e){let n=Yt(e);const t=[];for(;n.length>1&&n[0].isZero;)t.push(ee.ZERO),n=n.slice(1);if(n.length<=1)return{roots:t,rest:n};const r=n.reduce((e,n)=>e*n.d/X(e,n.d),1n),s=n.map(e=>e.n*r/e.d),o=Xt(s[0]),i=Xt(s[s.length-1]);if(!o||!i)return{roots:t,rest:n};const a=new Map;for(const e of o)for(const n of i)for(const t of[1n,-1n]){const r=new ee(t*e,n);a.set(r.toString(),r)}const u=[...a.values()].sort((e,n)=>e.compare(n));for(const e of u)for(;n.length>1&&Gt(n,e).isZero;)t.push(e),n=Kt(n,e);return{roots:t.sort((e,n)=>e.compare(n)),rest:n}}function nr(e){if("num"===e.k)return{coef:e.q,factors:[],notes:[]};if("atom"===e.k)return{coef:ee.ONE,factors:[{base:e,exp:Ut}],notes:[]};if("mul"===e.k){let n=e.c;const t=[],r=[];for(const s of e.fs)if("add"===s.b.k&&he(s.e)&&s.e.q.isInteger){const e=nr(s.b);n=n.mul(e.coef.pow(s.e.q.n));for(const n of e.factors)t.push({base:n.base,exp:Ne([n.exp,s.e])});r.push(...e.notes)}else t.push({base:s.b,exp:s.e});return tr(n,t,r)}return function(e){const n=Xe(e),t=[],r=n.map(e=>{const{coef:n,rest:t}=$e(e);return{coef:n,fs:t?en(t).map(sr):[]}}),s=r.reduce((e,n)=>({n:X(e.n,n.coef.n),d:e.d*n.coef.d/X(e.d,n.coef.d)}),{n:0n,d:1n});let o=new ee(s.n||1n,s.d);rr(e)<0&&(o=o.neg());const i=[];for(const e of r[0].fs){if(!he(e.e)||e.e.q.sign<=0)continue;let n=e.e.q;const t=ve(e.b);let s=!0;for(const e of r.slice(1)){const r=e.fs.find(e=>ve(e.b)===t&&he(e.e)&&e.e.q.sign>0);if(!r){s=!1;break}r.e.q.compare(n)<0&&(n=r.e.q)}s&&i.push({base:e.b,exp:ce(n)})}let a=e;if(!o.isOne||i.length){const e=Ne([ce(o),...i.map(e=>Ie(e.base,e.exp))]);a=Se(n.map(n=>Ne([n,Ae(e,ce(ee.MINUS_ONE))]))),t.push({id:"factor.common",values:{common:e}})}const u="add"===a.k?function(e,n){const t=[...or(e)];for(const r of t){const s=Jt(_t(e,r));if(!s||s.length<3)continue;if(1!==t.length)break;const{roots:o,rest:i}=er(s);if(!o.length)break;const a=[];let u=i[i.length-1]??ee.ONE;const l=i.map(e=>e.div(u)),c=[...o].sort((e,n)=>e.abs().compare(n.abs())||n.sign-e.sign);for(const e of c){const t=Qt([new ee(-e.n),new ee(e.d)],r);a.push({base:t,exp:Ut}),u=u.div(new ee(e.d)),n.push({id:"factor.root",values:{root:ce(e),factor:t}})}if(l.length>1){const e=l.reduce((e,n)=>e*n.d/X(e,n.d),1n),n=l.map(n=>n.mul(new ee(e))),t=n.reduce((e,n)=>X(e,n.n),0n)||1n;u=u.mul(new ee(t,e)),a.push({base:Qt(n.map(e=>e.div(new ee(t))),r),exp:Ut})}return tr(u,a,[])}const r=Xe(e);if(2===r.length){const[e,t]=r;for(const r of[2,3]){const s=ir(e,r),o=ir(Ce(t),r),i=3===r?ir(t,r):null;if(s&&o&&2===r){n.push({id:"factor.squares",values:{a:s,b:o}});const e=nr(Se([s,Ce(o)])),t=nr(Se([s,o]));return tr(e.coef.mul(t.coef),[...e.factors,...t.factors],[...e.notes,...t.notes])}if(3===r&&s&&(o||i)){const e=o??i,t=o?-1:1;n.push({id:t<0?"factor.cubes-difference":"factor.cubes-sum",values:{a:s,b:e}});const r=Se(t<0?[s,Ce(e)]:[s,e]),a=Se([Ae(s,ce(ee.of(2))),t<0?Ne([s,e]):Ce(Ne([s,e])),Ae(e,ce(ee.of(2)))]),u=nr(r);return tr(u.coef,[...u.factors,{base:a,exp:Ut}],u.notes)}}}return{coef:ee.ONE,factors:[{base:e,exp:Ut}],notes:[]}}(a,t):nr(a);return tr(o.mul(u.coef),[...i,...u.factors],[...t,...u.notes.filter(e=>!t.includes(e))])}(e)}function tr(e,n,t){const r=new Map;for(const e of n){const n=ve(e.base),t=r.get(n);t?t.exp=Se([t.exp,e.exp]):r.set(n,{...e})}return{coef:e,factors:[...r.values()].filter(e=>!(he(e.exp)&&e.exp.q.isZero)),notes:t}}const rr=e=>$e(Xe(e)[0]).coef.sign;const sr=e=>"mul"===e.k&&e.c.isOne&&1===e.fs.length?{b:e.fs[0].b,e:e.fs[0].e}:{b:e,e:Ut},or=e=>Ke(Qe(e));function ir(e,n){const{coef:t,rest:r}=$e(e);if(t.sign<0)return null;const s=function(e,n){const t=e=>{if(e<0n)return null;let t=0n,r=1n;for(;r**BigInt(n)<=e;)r*=2n;for(;t<r-1n;){const s=(t+r)/2n;s**BigInt(n)<=e?t=s:r=s}return t**BigInt(n)===e?t:null},r=t(e.n),s=t(e.d);return null!==r&&null!==s?new ee(r,s):null}(t,n);if(!s)return null;const o=[ce(s)];for(const e of r?en(r).map(sr):[]){if(!he(e.e))return null;const t=e.e.q.div(ee.of(n));if(!t.isInteger)return null;o.push(Ie(e.b,ce(t)))}return Ne(o)}function ar({coef:e,factors:n}){const t=[],r=[];for(const e of n){const n=he(e.exp)&&e.exp.q.sign<0,s=n?ce(e.exp.q.neg()):e.exp,o=Qe(Ie(e.base,s));(n?r:t).push(o)}const s=e=>e.reduce((e,n)=>e?C("*",[e,n]):n,null),o=e.n<0n?-e.n:e.n;let i=t.length?s(1n===o?t:[Qe(ce(new ee(o))),...t]):Qe(ce(new ee(o)));return e.n<0n&&(i=C("neg",[i])),1n!==e.d&&r.unshift(Qe(ce(new ee(e.d)))),r.length?C("/",[i,s(r)]):i}function ur(e){if("mul"!==e.k)return{num:e,den:Ut};const n=[ce(new ee(e.c.n))],t=[ce(new ee(e.c.d))];for(const r of e.fs)he(r.e)&&r.e.q.sign<0?t.push(Ie(r.b,ce(r.e.q.neg()))):n.push(Ie(r.b,r.e));return{num:Ne(n),den:Ne(t)}}function lr(e){if("add"!==e.k)return ur(e);const n=Xe(e).map(ur),t=new Map;let r=1n;for(const e of n){const{coef:n,rest:s}=$e(e.den);r=r*n.n/X(r,n.n);for(const e of s?en(s).map(sr):[]){const n=ve(e.b),r=t.get(n);(!r||he(e.e)&&he(r.e)&&e.e.q.compare(r.e.q)>0)&&t.set(n,e)}}const s=Ne([ce(new ee(r)),...[...t.values()].map(e=>Ie(e.b,e.e))]);if(he(s)&&s.q.isOne)return{num:e,den:Ut};return{num:Wt(Se(n.map(e=>Ne([e.num,s,Ae(e.den,ce(ee.MINUS_ONE))])))),den:s}}function cr(e){const{num:n,den:t}=lr(e);if(he(t)&&t.q.isOne)return null;const r=nr(n),s=nr(t),o=new Set(r.factors.map(e=>ve(e.base))),i=s.factors.filter(e=>o.has(ve(e.base))),a=[...r.factors,...s.factors.map(e=>({base:e.base,exp:Ce(e.exp)}))];return{...tr(r.coef.div(s.coef),a,[...r.notes,...s.notes]),cancelled:i.map(e=>e.base)}}const fr=(e,n)=>C("=",[e,n]),pr=Le(C("+",[Z("x"),C("neg",[Z("x")])])),dr=e=>V(e)&&"="===e.op;function hr(e,n,t,r={}){const s=[],o=(n,t,r,o,i,a={})=>s.push(e.step(n,t,r,o,i,a));o.raw=e=>s.push(e);const i=dr(n)?n:fr(n,pr),a=fr(Le(i.args[0]),Le(i.args[1]));a.s!==i.s&&o("normal","Work out each side",["tidy","arithmetic"],i,a);const u=[],l=[],c=new Set;u.numeric=c,wr(e,a,t,r,o,u,l,0);const f=[],p=new Set;for(const n of u.map(n=>c.has(n)?n:(n=>Bt(n,e.set.view("simplify"),{format:e.format}).node)(n))){const e=n.s;if(p.has(e))continue;p.add(e);const r=c.has(n)?null:vr(i,t,n);r?o("solve.reject","Check the answer",["check"],fr(Z(t),n),fr(Z(t),n),{value:n,reason:r}):f.push(n)}return f.sort(mr),{solutions:f,unsolved:l,steps:s,always:l.some(e=>e.always)}}function mr(e,n){const t=gr(e),r=gr(n);return null===t||null===r?e.s<n.s?-1:1:(Math.abs(t.im)>0)-(Math.abs(r.im)>0)||t.re-r.re||t.im-r.im}function gr(e){if(Ke(e).size)return null;try{const{re:n,im:t}=xt(e,{digits:20});return{re:n,im:t}}catch{return null}}function vr(e,n,t){const[r,s]=e.args,o=new Set([...Ke(r),...Ke(s),...Ke(t)]);if(o.delete(n),o.size)return null;const i=e=>br(e,n,t);try{const e=xt(i(r),{digits:40}),n=xt(i(s),{digits:40}),t=Math.max(1,Math.hypot(e.re,e.im),Math.hypot(n.re,n.im));return Math.hypot(e.re-n.re,e.im-n.im)>1e-12*t?"it does not make the two sides equal":null}catch(e){return`the original equation has no value there (${e.message})`}}function br(e,n,t){return"sym"===e.type?e.token===n?t:e:V(e)?C(e.op,e.args.map(e=>br(e,n,t))):e}const yr=24;function wr(e,n,t,r,s,o,i,a){if(a>yr)return Nr(e,n,t,r,s,o,i);const u=function(e){const n=[],t=e=>{if(V(e)&&"|"===e.op)return void e.args.forEach(t);if(!dr(e))return void n.push(e);const r=e.args[1],s=H(r,"pre").find(({node:e})=>V(e)&&"pm"===e.op);if(!s)return void n.push(e);const o=s.node.args[0],i=Le(J(r,s.path,o)),a=Le(J(r,s.path,C("neg",[o])));t(fr(e.args[0],i)),a.s!==i.s&&t(fr(e.args[0],a))};return t(e),n}(n);if(u.length>1){s("solve.split","Two cases",["split"],n,u.reduce((e,n)=>e?C("|",[e,n]):n,null));for(const n of u)wr(e,n,t,r,s,o,i,a+1);return}if(1===u.length&&u[0].s!==n.s)return wr(e,u[0],t,r,s,o,i,a+1);const[l,c]=n.args;if("sym"===l.type&&l.token===t&&!Ht(We(c),t))return void o.push(c);const f=Se([We(l),Ce(We(c))]);if(!Ht(f,t))return he(f)&&f.q.isZero?(s("solve.always","Always true",["check"],n,n,{v:t}),void i.push(Object.assign(Object.create(n),{always:!0}))):void s("solve.never","No solution",["check"],n,n,{v:t});const p=Xe(f).filter(e=>Ht(e,t)),d=Xe(f).filter(e=>!Ht(e,t)),h=Se(p),m=Ce(Se(d));if(ie.complex){const u=_t(f,t);if(u&&[...u.values()].every(e=>!Ht(e,t))&&Math.max(...u.keys())>=2){const l=fr(Qe(f),pr);return l.s!==n.s&&s("solve.one-side","Everything on one side",["move"],n,l,{v:t}),Sr(e,f,u,t,r,s,o,i,a,l)}}if(1===p.length)return xr(e,n,h,m,t,r,s,o,i,a);const g=fr(Qe(f),pr),v=_t(f,t);if(v&&[...v.values()].every(e=>!Ht(e,t))&&Math.max(...v.keys())>=2)return g.s!==n.s&&s("solve.one-side","Everything on one side",["move"],n,g,{v:t}),Sr(e,f,v,t,r,s,o,i,a,g);const b=nr(f),y=b.factors.filter(e=>Ht(e.base,t)&&!(he(e.exp)&&e.exp.q.sign<0));if(y.length>=2||1===y.length&&b.factors.length>1){const u=fr(ar(b),pr);g.s!==n.s&&s("solve.one-side","Everything on one side",["move"],n,g,{v:t}),u.s!==g.s&&s("solve.factor","Factor",["factor"],g,u,{});const l=y.map(e=>fr(Qe(e.base),pr));l.length>1&&s("solve.zero-product","A product is zero when a factor is",["split"],u,l.reduce((e,n)=>e?C("|",[e,n]):n,null),{});for(const n of l)wr(e,n,t,r,s,o,i,a+1);return}const w=lr(f);if(Ht(w.den,t)){const u=fr(Qe(w.num),pr);return s("solve.clear-denominator","Clear the denominator",["move"],n,u,{den:Qe(w.den)}),wr(e,u,t,r,s,o,i,a+1)}const x=p.filter(e=>!$r(e,t));if(1===x.length&&/\b(sqrt|root|abs)\b/.test(Qe(x[0]).s))return xr(e,n,x[0],Ce(Se(Xe(f).filter(e=>e!==x[0]))),t,r,s,o,i,a);const q=function(e,n){const t=Qe(e),r=new Map;for(const{node:e}of H(t,"pre"))if(V(e)&&Ke(e).has(n))if("^"!==e.op||Ke(e.args[0]).has(n))["sin","cos","tan","ln","sqrt","abs"].includes(e.op)&&e.args[0].s===n&&r.set(e.s,e);else{const t=We(e.args[1]),{coef:s,rest:o}=$e(t);if(o&&s.isInteger&&Qe(o).s===n){const t=Qe(Ae(We(e.args[0]),We(Z(n))));r.set(t.s,t)}}const s="u"===n?"w":"u";for(const e of r.values()){const r=Or(t,e,n,s);if(!r||Ke(r).has(n))continue;const o=We(r),i=_t(o,s);if(i&&Math.max(...i.keys())>=1)return{kernel:We(e),poly:o,symbol:s}}return null}(f,t);if(q){const n=q.symbol;s("solve.substitute","Substitute",["substitute"],g,fr(Qe(q.poly),pr),{u:Z(n),kernel:Qe(q.kernel)});const u=[];wr(e,fr(Qe(q.poly),pr),n,r,s,u,i,a+1);for(const n of u)wr(e,fr(Qe(q.kernel),n),t,r,s,o,i,a+1);return}return Nr(e,g,t,r,s,o,i)}function xr(e,n,t,r,s,o,i,a,u,l){let c=fr(Qe(t),Qe(r));c.s!==n.s&&i("solve.move","Move the terms",["move"],n,c,{v:s});const{coef:f,rest:p}=$e(t),d=p?en(p):[],h=Ne([ce(f),...d.filter(e=>!Ht(e,s))]);if(!he(h)||!h.q.isOne){const e=Ne(d.filter(e=>Ht(e,s))),n=fr(Qe(e),Qe(Fe(r,h))),t=he(h)&&!h.q.isInteger&&(1n===h.q.n||-1n===h.q.n);i(t?"solve.multiply":"solve.divide",t?"Multiply both sides":"Divide both sides",["move"],c,n,{by:Qe(t?ce(h.q.inv()):h),v:s}),c=n}if("sym"===c.args[0].type&&c.args[0].token===s)return void a.push(c.args[1]);const m=We(c.args[1]);if(he(m)&&m.q.isZero){const n=We(c.args[0]),t=en(n).filter(e=>Ht(e,s)&&!qr(e));if(t.length>1||t.length<en(n).filter(e=>Ht(e,s)).length){const n=t.map(e=>fr(Qe(kr(e)),pr));i("solve.zero-product","A product is zero when a factor is",["split"],c,n.reduce((e,n)=>e?C("|",[e,n]):n,null),{});for(const t of n)wr(e,t,s,o,i,a,u,l+1);return}}const g=Bt(c,e.set.view("isolate"),{format:e.format,maxSteps:1,variable:s});return g.steps.length?(function(e,n){for(const t of n)e.raw(t)}(i,g.steps),wr(e,g.node,s,o,i,a,u,l+1)):Nr(e,c,s,o,i,a,u)}const qr=e=>"mul"===e.k&&e.c.isOne&&1===e.fs.length&&he(e.fs[0].e)&&e.fs[0].e.q.sign<0,kr=e=>"mul"===e.k&&e.c.isOne&&1===e.fs.length&&he(e.fs[0].e)&&e.fs[0].e.q.sign>0?e.fs[0].b:e,$r=(e,n)=>{const t=_t(e,n);return Boolean(t)&&[...t.values()].every(e=>!Ht(e,n))};function Sr(e,n,t,r,s,o,i,a,u,l){const c=Math.max(...t.keys()),f=Jt(t);if(ie.complex&&c>=3&&c<=64&&f&&f.every((e,n)=>0===n||n===c||e.isZero)&&!f[0].isZero){const n=f[0].neg().div(f[c]),t=Ae(ce(n.abs()),ce(new ee(1n,BigInt(c)))),s=[];for(let r=0;r<c;r++){const o=Ne([ce(new ee(BigInt((n.sign<0?1:0)+2*r),BigInt(c))),We(Z("pi"))]),i=C("*",[Qe(t),C("+",[C("cos",[Qe(o)]),C("*",[Z("i"),C("sin",[Qe(o)])])])]);s.push(Bt(i,e.set.view("simplify"),{format:e.format}).node)}return o("solve.roots","All n roots",["complex"],l,s.map(e=>fr(Z(r),e)).reduce((e,n)=>e?C("|",[e,n]):n),{v:r,n:String(c),last:String(c-1),r:Qe(ce(n))}),void i.push(...s)}if(f){const{roots:t}=er(f);if(t.length){const t=nr(n),c=fr(ar(t),pr);if(t.factors.length>1||t.factors.some(e=>!(he(e.exp)&&e.exp.q.isOne))){o("solve.factor","Factor",["factor"],l,c,{});const n=t.factors.filter(e=>Ht(e.base,r)).map(e=>fr(Qe(e.base),pr));n.length>1&&o("solve.zero-product","A product is zero when a factor is",["split"],c,n.reduce((e,n)=>e?C("|",[e,n]):n,null),{});for(const t of n)wr(e,t,r,s,o,i,a,u+1);return}}}if(2===c){const e=t.get(2),n=t.get(1)??ce(ee.ZERO),s=t.get(0)??ce(ee.ZERO),a=Se([Ae(n,ce(ee.of(2))),Ne([ce(ee.of(-4)),e,s])]);if(he(a)&&a.q.sign<0&&!ie.complex)return void o("solve.no-real","No real solutions",["quadratic"],l,l,{D:Qe(a)});const u=Ne([ce(ee.of(2)),e]),c=Ae(a,ce(ee.HALF)),f=Qe(Fe(Se([Ce(n),c]),u)),p=Qe(Fe(Se([Ce(n),Ce(c)]),u)),d=f.s===p.s?fr(Z(r),f):C("|",[fr(Z(r),f),fr(Z(r),p)]);return o(ie.complex&&he(a)&&a.q.sign<0?"solve.quadratic-complex":"solve.quadratic","The quadratic formula",["quadratic"],l,d,{a:Qe(e),b:Qe(n),c:Qe(s),D:Qe(a),root:Qe(c),v:r}),i.push(f),void(p.s!==f.s&&i.push(p))}return Nr(e,l,r,s,o,i,a)}function Or(e,n,t,r){if(e.s===n.s)return Z(r);if(!V(e))return e;if("^"===e.op&&"^"===n.op&&e.args[0].s===n.args[0].s){const n=We(e.args[1]),{coef:s,rest:o}=$e(n);if(o&&Qe(o).s===t&&s.isInteger)return C("^",[Z(r),Qe(ce(s))])}return C(e.op,e.args.map(e=>Or(e,n,t,r)))}function Nr(e,n,t,r,s,o,i){if(r.approximate){const e="number"==typeof r.approximate?r.approximate:15,i=We(C("-",[n.args[0],n.args[1]])),a=ie.complex?Jt(_t(i,t)):null,u=a&&a.length>2?function(e,n){const t=e.length-1,r=e[t],s=e.map(e=>e/r),o=(e,n)=>{let r=1,o=0;for(let i=t-1;i>=0;i--)[r,o]=[r*e-o*n+s[i],r*n+o*e];return[r,o]};let i=Array.from({length:t},(e,n)=>[1.3*Math.cos(2*Math.PI*n/t+.4),1.3*Math.sin(2*Math.PI*n/t+.4)]);for(let e=0;e<500;e++){let e=0;if(i=i.map(([n,r],s)=>{let[a,u]=o(n,r);for(let e=0;e<t;e++){if(e===s)continue;const t=n-i[e][0],o=r-i[e][1],l=t*t+o*o;[a,u]=[(a*t+u*o)/l,(u*t-a*o)/l]}return e=Math.max(e,Math.hypot(a,u)),[n-a,r-u]}),e<1e-15)break}const a=e=>Math.abs(e)<1e-12?0:Number(e.toPrecision(n));return i.map(([e,n])=>[a(e),a(n)]).sort((e,n)=>(0!==e[1])-(0!==n[1])||e[0]-n[0]||e[1]-n[1]).map(([e,n])=>{const t=Ar(e);if(0===n)return t;const r=C("*",[Ar(Math.abs(n)),Z("i")]);return 0===e?n<0?C("neg",[r]):r:C(n<0?"-":"+",[t,r])})}(a.map(e=>e.toNumber()),Math.min(e,14)):function(e,n,t){const r=We(C("-",[e.args[0],e.args[1]])),s=Qe(r),o=e=>{try{return xt(br(s,n,Ir(e)),{digits:30}).number}catch{return NaN}},i=[],a=[];for(let e=-400;e<=400;e++)a.push(Math.sign(e)*(Math.abs(e)/20)**2);for(let e=0;e+1<a.length;e++){let n=a[e],r=a[e+1],s=o(n);const u=o(r);if(Number.isNaN(s)||Number.isNaN(u))continue;if(0===s){i.push(Ir(n));continue}if(s*u>0)continue;for(let e=0;e<200&&r-n>1e-16*Math.abs(n);e++){const e=(n+r)/2,t=o(e);if(Number.isNaN(t))break;s*t<=0?r=e:(n=e,s=t)}const l=(n+r)/2;Math.abs(o(l))<1e-6&&i.push(Ir(Number(l.toPrecision(Math.min(t,15)))))}return i}(n,t,e);for(const e of u)o.numeric?.add(e);return s("solve.numeric","Numeric solutions",["numeric"],n,u.length?u.map(e=>fr(Z(t),e)).reduce((e,n)=>e?C("|",[e,n]):n):n,{v:t,count:String(u.length)}),void o.push(...u)}s("solve.unsolved","No exact method",["unsolved"],n,n,{v:t}),i.push(n)}const Ar=e=>z(String(e).includes("e")?e.toFixed(20):String(e)),Ir=e=>z(Number.isInteger(e)?String(e):e.toFixed(17));const Mr=(e,n)=>C("int",[e,"string"==typeof n?Z(n):n]);function Er(e,n,t,r=0){if(r>6)throw new Q("integration went too deep");const s=[],o=e.set.view("integrate"),i=Bt(n,e.set.view("simplify"),{format:e.format});s.push(...i.steps);let a=Le(Mr(i.node,t));for(let n=0;n<40;n++){const n=Bt(a,o,{format:e.format});s.push(...n.steps),a=n.node;const t=H(a,"pre").find(({node:e})=>{return V(n=e)&&"int"===n.op&&2===n.args.length;var n});if(!t)return{node:a,steps:s};const i=Tr(e,t.node,r);if(!i)throw new Q(`no rule for ∫ ${e.format(t.node.args[0])} d${t.node.args[1].token}`);s.push(...i.steps),a=Le(J(a,t.path,i.node))}throw new Q("integration did not finish")}function Tr(e,n,t){const[r,s]=n.args,o=s.token,i=We(r),a=function(e,n){const{num:t,den:r}=lr(e),s=nr(t),o=nr(r),i=s.coef.div(o.coef);if(i.isOne||!Ht(e,n))return null;const a=Ne([...s.factors.map(e=>Ie(e.base,e.exp)),...o.factors.map(e=>Ie(e.base,Ce(e.exp)))]);return{c:i,rest:a}}(i,o);if(a){const t=C("*",[Qe(ce(a.c)),Mr(Qe(a.rest),s)]);return{node:t,steps:[e.step("integrate.constant","Take the number out",["split"],n,t,{c:Qe(ce(a.c))})]}}const u=Wt(i);if(Qe(u).s!==r.s&&"add"===u.k){const t=Mr(Qe(u),s);return{node:t,steps:[e.step("integrate.expand","Multiply out",["expand"],n,t,{})]}}const l=function(e,n){const t=new Map;for(const{node:r}of H(e,"pre")){if(!V(r))continue;const e="^"===r.op?r.args:"/"===r.op?[r.args[1]]:"der"===r.op||"int"===r.op?[]:1===r.args.length||"root"===r.op?[r.args[r.args.length-1]]:[];for(const r of e){if("sym"===r.type||!Ke(r).has(n))continue;const e=_t(We(r),n);if(!e||1!==Math.max(...e.keys()))continue;const s=e.get(1);Ht(s,n)||he(s)&&s.q.isZero||t.set(r.s,{node:r,a:s})}}for(const r of t.values()){if(!Ke(jr(e,r.node,Z("§"))).has(n))return r}return null}(r,o);if(l){const s="u"===o?"t":"u",i=jr(r,l.node,Z(s));if(!Ke(i).has(o)){const r=Er(e,Le(i),s,t+1),a=Le(br(r.node,s,l.node)),u=Qe(Fe(We(a),l.a));return{node:u,steps:[e.step("integrate.substitute","Substitute",["substitute"],n,Le(Mr(i,s)),{u:Z(s),inside:l.node,a:Qe(l.a),v:o}),...r.steps,e.step("integrate.back","Substitute back",["substitute"],r.node,u,{u:Z(s),inside:l.node})]}}}const c=function(e,n,t,r){const{num:s,den:o}=lr(t);if(!Ht(o,r))return null;const i=Jt(_t(s,r)),a=Jt(_t(o,r));if(!i||!a||a.length<2)return null;const u=We(Z(r)),l=e=>Qt(e,r);if(i.length>=a.length){const{q:t,r:s}=function(e,n){const t=Yt(n);let r=Yt(e);const s=t[t.length-1];if(s.isZero)throw new Q("division by the zero polynomial");const o=new Array(Math.max(1,r.length-t.length+1)).fill(ee.ZERO);for(;r.length>=t.length&&(1!==r.length||!r[0].isZero);){const e=r.length-t.length,n=r[r.length-1].div(s);if(o[e]=n,r=Yt(r.map((r,s)=>s>=e&&s-e<t.length?r.sub(n.mul(t[s-e])):r)),r.length-1<e+t.length-1&&r[r.length-1].isZero)break;if(r.length<t.length)break}return{q:Yt(o),r:r}}(i,a),o=Mr(Qe(Se([l(t),Fe(l(s),l(a))])),Z(r));return{node:o,steps:[e.step("integrate.divide","Divide the polynomials",["divide"],n,o,{})]}}const{roots:c,rest:f}=er(a),p=new Set(c.map(String)).size===c.length;if(!c.length||!p||1!==f.length)return null;const d=(e=>e.length<=1?[ee.ZERO]:e.slice(1).map((e,n)=>e.mul(ee.of(n+1))))(a),h=c.map(e=>Fe(ce(Gt(i,e).div(Gt(d,e))),Se([u,ce(e.neg())]))),m=Mr(Qe(Se(h)),Z(r));return{node:m,steps:[e.step("integrate.partial-fractions","Partial fractions",["partial-fractions"],n,m,{})]}}(e,n,i,o);return c||function(e,n,t,r,s){if(s>=6)return null;const o=en(t),i=o.filter(e=>Zr(e,r)||!Ht(e,r)),a=o.filter(e=>!i.includes(e));if(1!==a.length)return null;const u=Ne(i),l=a[0];if(!Ht(u,r))return null;const c=Qe(l),f=V(c)?c.op:"";let p,d,h;"ln"===f||"arctan"===f||"arcsin"===f?(p=l,d=u):(p=u,d=l);try{h=Er(e,Qe(d),r,s+1)}catch{return null}const m=We(e.differentiate(Qe(p),r)),g=Ne([p,We(h.node)]),v=Mr(Qe(Ne([m,We(h.node)])),Z(r)),b=Le(C("-",[Qe(g),v]));return{node:b,steps:[e.step("integrate.parts","Integration by parts",["parts"],n,b,{u:Qe(p),dv:Qe(d),V:h.node,du:Qe(m),v:r})]}}(e,n,i,o,t)}function jr(e,n,t){return e.s===n.s?t:V(e)?C(e.op,e.args.map(e=>jr(e,n,t))):e}const Zr=(e,n)=>{const t=_t(e,n);return Boolean(t)&&[...t.values()].every(e=>!Ht(e,n))&&Math.max(...t.keys())>=1};function zr(e,n,t,r,s,o){for(const e of[r,s])if(Ke(e).has("∞")||/∞|infinity/.test(e.s))throw new Q("integrals to infinity need limits, which this CAS does not do yet");if(o&&!Ke(r).size&&!Ke(s).size){const i=xt(r,{digits:20}).number,a=xt(s,{digits:20}).number,u=[...Ke(o)].filter(e=>e!==t);if(!u.length){const[r,s]=i<=a?[i,a]:[a,i];for(const i of function(e,n){const t=[],r={tan:"cos",sec:"cos",csc:"sin",cot:"sin"};return H(e,"pre").forEach(({node:e})=>{if(!V(e))return;let s=null;"/"===e.op?s=e.args[1]:"^"===e.op&&"num"===e.args[1].type&&e.args[1].value<0||"^"===e.op&&V(e.args[1])&&"neg"===e.args[1].op?s=e.args[0]:r[e.op]&&(s=C(r[e.op],[e.args[0]])),s&&Ke(s).has(n)&&!t.some(e=>e.s===s.s)&&t.push(s)}),t}(o,t)){let o=[];try{o=e.solve(C("=",[i,z("0")]),t,{steps:!1}).solutions}catch{continue}for(const i of o){const o=xt(i,{digits:20});if(!Number.isFinite(o.number)||o.im||o.number<r||o.number>s)continue;if(o.number>r&&o.number<s)throw new Q(`the integrand has no value at ${t} = ${e.format(i)}, between the limits, so F(b) - F(a) would be wrong`);let a=NaN;try{a=xt(Le(br(n,t,i)),{digits:15}).number}catch{}if(!Number.isFinite(a))throw new Q(`the integral diverges at ${t} = ${e.format(i)}: F has no value at that limit`)}}}for(let e=0;!u.length&&e<=64;e++){const n=i+(a-i)*(e+.5)/65;try{const e=xt(br(o,t,z(n.toFixed(15))),{digits:15}).number;if(!Number.isFinite(e))throw new Q("not finite")}catch(e){throw new Q(`the integrand has no value at ${t} = ${n.toPrecision(6)}, between the limits, so F(b) - F(a) would be wrong`)}}}const i=e=>Le(br(n,t,e)),a=Le(C("-",[i(s),i(r)])),u=C("-",[i(s),i(r)]);return{value:a,step:e.step("integrate.evaluate","Between the limits",["evaluate"],n,u,{a:r,b:s,F:n,v:t})}}const Dr={"=":1,"≠":1,"<":1,">":1,"≤":1,"≥":1,"+":2,"-":2,"±":2,"*":3,"/":3,neg:4,pm:4,"^":5},Br={pi:"π","≤":"≤","≥":"≥"};function Cr(e){return Pr(e).text}function Fr(e,n,t=!1){return e.prec<n||t&&e.prec===n?`(${e.text})`:e.text}const Vr=e=>"op"===e.type?e.op:null,Rr=e=>e.text.startsWith("-"),Ur={"()":["(",")"],"[]":["[","]"],"{}":["{","}"],"(]":["(","]"],"[)":["[",")"]};function Pr(e){if("num"===e.type)return{text:e.token,prec:e.token.startsWith("-")?4:9};if("sym"===e.type)return{text:Br[e.token]??e.token,prec:9};const n=e.args.map(Pr),t=Dr[e.op];switch(e.op){case"+":case"-":case"±":{const r="+"===e.op&&["+","-"].includes(Vr(e.args[1])),s=Rr(n[1])?`(${n[1].text})`:r?n[1].text:Fr(n[1],t,!0);return{text:`${Fr(n[0],t)} ${"±"===e.op?"±":e.op} ${s}`,prec:t}}case"*":{const r=Fr(n[0],t),s=Rr(n[1])?`(${n[1].text})`:"*"===Vr(e.args[1])?n[1].text:Fr(n[1],t,!0);return{text:"num"===e.args[0].type&&"num"!==e.args[1].type&&!/^[\d(-]/.test(s)||"sym"===e.args[1].type&&"i"===e.args[1].token&&!/[\dA-Za-z]$/.test(r)?`${r}${s}`:`${r}·${s}`,prec:t}}case"/":return{text:`${Fr(n[0],t)}/${Fr(n[1],t,!0)}`,prec:t};case"^":return{text:`${Fr(n[0],t+1)}^${Fr(n[1],t)}`,prec:t};case"neg":return{text:`-${Fr(n[0],3)}`,prec:4};case"pm":return{text:`±${Fr(n[0],4)}`,prec:4};case"!":return{text:`${Fr(n[0],9)}!`,prec:6};case"abs":return{text:`|${n[0].text}|`,prec:9};case"root":return{text:`root(${n[0].text}, ${n[1].text})`,prec:9};case"der":return{text:`d/d${n[1].text} (${n[0].text})`,prec:6};case"int":return{text:4===n.length?`∫ from ${n[2].text} to ${n[3].text} of ${n[0].text} d${n[1].text}`:`∫ ${n[0].text} d${n[1].text}`,prec:1};case"|":return{text:`${n[0].text} or ${n[1].text}`,prec:0};case"_":return{text:`${Fr(n[0],9)}_${n[1].prec>=9&&!/\s/.test(n[1].text)?n[1].text:`(${n[1].text})`}`,prec:9};case"()":case"[]":case"{}":case"(]":case"[)":return{text:`${Ur[e.op][0]}${n.map(e=>e.text).join(", ")}${Ur[e.op][1]}`,prec:9};case"&":return{text:`${n[0].text} and ${n[1].text}`,prec:0};default:return void 0!==t?{text:`${Fr(n[0],t)} ${e.op} ${Fr(n[1],t,!0)}`,prec:t}:{text:`${e.op.replace(/^\\:?/,"")}(${n.map(e=>e.text).join(", ")})`,prec:9}}}const Wr=['name: Differentiation\n\n# `u v der` means d(u)/dv (MathZetta\'s postfix: mathzetta/docs/postfix.md). Every rule here rewrites one derivative into smaller\n# derivatives, so differentiating always terminates: the expressions sitting\n# under a `D` get strictly smaller at every step.\n#\n# The variable is a slot, not a fixed letter, so these rules work for d/dx,\n# d/dt or anything else. That is what freeOf and dependsOn are for: they ask\n# whether one operand mentions whatever the variable slot bound to.\n#\n# Where two rules share a structure key, requirements keep them apart and\n# priority picks the one that reads better. The constant-multiple rule outranks\n# the product rule for exactly that reason.\n\nrules:\n  # ---- the two base cases ----------------------------------------------\n  - id: d-constant\n    name: Derivative of a constant\n    mode: differentiate\n    pattern: c v der\n    require:\n      - [freeOf, c, v]\n    output: "0"\n    priority: 10\n    explain: "{c} does not change when {v} changes, so its rate of change is 0."\n    tags: [base]\n\n  - id: d-variable\n    name: Derivative of the variable\n    mode: differentiate\n    pattern: u v der\n    require:\n      - [sameAs, u, v]\n    output: "1"\n    priority: 10\n    explain: "{v} changes exactly as fast as itself, so the rate is 1."\n    tags: [base]\n\n  # ---- splitting a derivative into smaller ones ------------------------\n  - id: d-sum\n    name: Sum rule\n    mode: differentiate\n    pattern: u w + v der\n    output: u v der w v der +\n    explain: "Differentiate {u} and {w} separately, then add the two results."\n    tags: [split]\n\n  - id: d-difference\n    name: Difference rule\n    mode: differentiate\n    pattern: u w - v der\n    output: u v der w v der -\n    explain: "Differentiate {u} and {w} separately, then subtract the second from the first."\n    tags: [split]\n\n  - id: d-negate\n    name: Derivative of a negative\n    mode: differentiate\n    pattern: u neg v der\n    output: u v der neg\n    explain: "The minus sign just carries through to the derivative of {u}."\n    tags: [split]\n\n  - id: d-constant-multiple\n    name: Constant multiple rule\n    mode: differentiate\n    pattern: c u * v der\n    commutative: true\n    require:\n      - [freeOf, c, v]\n      - [dependsOn, u, v]\n    output: c u v der *\n    priority: 5\n    explain: "{c} is only a scale factor, so it comes along unchanged while {u} is differentiated."\n    tags: [split]\n\n  - id: d-constant-divisor\n    name: Dividing by a constant\n    mode: differentiate\n    pattern: u c / v der\n    require:\n      - [freeOf, c, v]\n    output: u v der c /\n    priority: 5\n    explain: "Dividing by {c} is only a scale factor, so differentiate {u} and keep the division."\n    tags: [split]\n\n  - id: d-product\n    name: Product rule\n    mode: differentiate\n    pattern: u w * v der\n    require:\n      - [dependsOn, u, v]\n      - [dependsOn, w, v]\n    output: u v der w * u w v der * +\n    explain: "Differentiate {u} and keep {w}, then keep {u} and differentiate {w}, and add the two."\n    tags: [split]\n\n  - id: d-quotient\n    name: Quotient rule\n    mode: differentiate\n    pattern: u w / v der\n    output: u v der w * u w v der * - w 2 ^ /\n    explain: "Bottom times the derivative of the top, minus top times the derivative of the bottom, all over {w} squared."\n    tags: [split]\n\n  # ---- powers -----------------------------------------------------------\n  - id: d-power-of-variable\n    name: Power rule\n    mode: differentiate\n    pattern: u n ^ v der\n    require:\n      - [sameAs, u, v]\n      - [freeOf, n, v]\n    output: n u n 1 - ^ *\n    fold: true\n    priority: 6\n    explain: "Bring the {n} down in front and drop the power by one."\n    tags: [power]\n\n  - id: d-power-chain\n    name: Power rule with the chain rule\n    mode: differentiate\n    pattern: u n ^ v der\n    require:\n      - [freeOf, n, v]\n      - [dependsOn, u, v]\n    output: n u n 1 - ^ * u v der *\n    fold: true\n    explain: "Bring the {n} down and drop the power by one, then multiply by the derivative of {u}."\n    tags: [power, chain]\n\n  - id: d-exponential-base\n    name: Derivative of a constant to a power\n    mode: differentiate\n    pattern: a u ^ v der\n    require:\n      - [freeOf, a, v]\n      - [dependsOn, u, v]\n    output: a u ^ a ln * u v der *\n    explain: "A constant base gives itself back, multiplied by the natural log of {a} and the derivative of {u}."\n    tags: [power, chain]\n\n  # ---- the standard functions, each with its chain-rule factor ---------\n  # e^u is `=e u ^`: the base is the literal symbol e. Outranks the general\n  # a^u rule (which would give e^u · ln e · u\').\n  - id: d-exp\n    name: Derivative of e to a power\n    mode: differentiate\n    pattern: =e u ^ v der\n    require:\n      - [dependsOn, u, v]\n    output: =e u ^ u v der *\n    priority: 7\n    explain: "e to a power is its own derivative, times the derivative of {u}."\n    tags: [chain]\n\n  - id: d-ln\n    name: Derivative of a natural log\n    mode: differentiate\n    pattern: u ln v der\n    output: u v der u /\n    explain: "The derivative of {u} over {u} itself."\n    tags: [chain]\n\n  - id: d-sqrt\n    name: Derivative of a square root\n    mode: differentiate\n    pattern: u sqrt v der\n    output: u v der 2 u sqrt * /\n    explain: "The derivative of {u}, over twice the square root of {u}."\n    tags: [chain]\n\n  - id: d-sin\n    name: Derivative of sine\n    mode: differentiate\n    pattern: u sin v der\n    output: u cos u v der *\n    explain: "Sine differentiates to cosine, times the derivative of {u}."\n    tags: [chain, trig]\n\n  - id: d-cos\n    name: Derivative of cosine\n    mode: differentiate\n    pattern: u cos v der\n    output: u sin u v der * neg\n    explain: "Cosine differentiates to minus sine, times the derivative of {u}."\n    tags: [chain, trig]\n\n  - id: d-tan\n    name: Derivative of tangent\n    mode: differentiate\n    pattern: u tan v der\n    output: u v der u cos 2 ^ /\n    explain: "The derivative of {u}, over cosine squared of {u}."\n    tags: [chain, trig]\n\n  # ---- the rest of the standard functions -----------------------------\n  - id: d-csc\n    name: Derivative of cosecant\n    mode: differentiate\n    pattern: u csc v der\n    output: u csc u cot * u v der * neg\n    explain: "Cosecant differentiates to minus cosecant times cotangent, times the derivative of {u}."\n    tags: [chain, trig]\n\n  - id: d-sec\n    name: Derivative of secant\n    mode: differentiate\n    pattern: u sec v der\n    output: u sec u tan * u v der *\n    explain: "Secant differentiates to secant times tangent, times the derivative of {u}."\n    tags: [chain, trig]\n\n  - id: d-cot\n    name: Derivative of cotangent\n    mode: differentiate\n    pattern: u cot v der\n    output: u v der u sin 2 ^ / neg\n    explain: "Minus the derivative of {u}, over sine squared of {u}."\n    tags: [chain, trig]\n\n  - id: d-arcsin\n    name: Derivative of arcsine\n    mode: differentiate\n    pattern: u arcsin v der\n    output: u v der 1 u 2 ^ - sqrt /\n    explain: "The derivative of {u}, over the square root of 1 minus {u} squared."\n    tags: [chain, trig]\n\n  - id: d-arccos\n    name: Derivative of arccosine\n    mode: differentiate\n    pattern: u arccos v der\n    output: u v der 1 u 2 ^ - sqrt / neg\n    explain: "Minus the derivative of {u}, over the square root of 1 minus {u} squared."\n    tags: [chain, trig]\n\n  - id: d-arctan\n    name: Derivative of arctangent\n    mode: differentiate\n    pattern: u arctan v der\n    output: u v der 1 u 2 ^ + /\n    explain: "The derivative of {u}, over 1 plus {u} squared."\n    tags: [chain, trig]\n\n  - id: d-sinh\n    name: Derivative of hyperbolic sine\n    mode: differentiate\n    pattern: u sinh v der\n    output: u cosh u v der *\n    explain: "Hyperbolic sine differentiates to hyperbolic cosine, times the derivative of {u}."\n    tags: [chain, hyperbolic]\n\n  - id: d-cosh\n    name: Derivative of hyperbolic cosine\n    mode: differentiate\n    pattern: u cosh v der\n    output: u sinh u v der *\n    explain: "Hyperbolic cosine differentiates to hyperbolic sine — no minus sign, unlike cosine — times the derivative of {u}."\n    tags: [chain, hyperbolic]\n\n  - id: d-tanh\n    name: Derivative of hyperbolic tangent\n    mode: differentiate\n    pattern: u tanh v der\n    output: u v der u cosh 2 ^ /\n    explain: "The derivative of {u}, over hyperbolic cosine squared of {u}."\n    tags: [chain, hyperbolic]\n\n  - id: d-arcsinh\n    name: Derivative of inverse hyperbolic sine\n    mode: differentiate\n    pattern: u arcsinh v der\n    output: u v der u 2 ^ 1 + sqrt /\n    explain: "The derivative of {u}, over the square root of {u} squared plus 1."\n    tags: [chain, hyperbolic]\n\n  - id: d-arccosh\n    name: Derivative of inverse hyperbolic cosine\n    mode: differentiate\n    pattern: u arccosh v der\n    output: u v der u 2 ^ 1 - sqrt /\n    explain: "The derivative of {u}, over the square root of {u} squared minus 1."\n    tags: [chain, hyperbolic]\n\n  - id: d-arctanh\n    name: Derivative of inverse hyperbolic tangent\n    mode: differentiate\n    pattern: u arctanh v der\n    output: u v der 1 u 2 ^ - /\n    explain: "The derivative of {u}, over 1 minus {u} squared."\n    tags: [chain, hyperbolic]\n\n  # log base a: `a u log` in MathZetta\'s postfix\n  - id: d-log\n    name: Derivative of a logarithm\n    mode: differentiate\n    pattern: a u log v der\n    require:\n      - [freeOf, a, v]\n    output: u v der u a ln * /\n    explain: "The derivative of {u}, over {u} times the natural log of {a}."\n    tags: [chain]\n\n  - id: d-abs\n    name: Derivative of an absolute value\n    mode: differentiate\n    pattern: u abs v der\n    output: u u abs / u v der *\n    explain: "The absolute value keeps the sign of {u}: {u} over its absolute value, times the derivative of {u}."\n    tags: [chain]\n\n  # n-th root: `n u root`\n  - id: d-root\n    name: Derivative of a root\n    mode: differentiate\n    pattern: n u root v der\n    require:\n      - [freeOf, n, v]\n    output: u v der n n u root n 1 - ^ * /\n    fold: true\n    explain: "The derivative of {u}, over {n} times the {n}th root of {u} to the power {n} minus 1."\n    tags: [chain]\n','name: More derivatives\n\n# Derivative rules the CAS adds to rules/calculus.yaml (which calcvidnew also uses).\n\nrules:\n  - id: d-general-power\n    name: A power with the variable in both places\n    mode: differentiate\n    pattern: u w ^ v der\n    require:\n      - [dependsOn, u, v]\n      - [dependsOn, w, v]\n    output: u w ^ w v der u ln * w u v der * u / + *\n    explain: "Write {u}^{w} as e^({w} ln {u}) and use the chain rule: the derivative is {u}^{w} times the derivative of {w} ln {u}."\n    tags: [chain]\n','name: Integrals\n\n# mode: integrate. `f v 2 int` is ∫ f dv (MathZetta\'s postfix; the 2 counts f and v).\n# Each rule turns one integral into smaller ones, or into an answer. The CAS also\n# tries, when no rule fits: multiplying out, a linear substitution (∫ f(ax + b) dx\n# from the rule for ∫ f(u) du), and integration by parts for polynomial × e^, sin,\n# cos, ln, arctan. Definite integrals are the antiderivative between the limits.\n#\n# Patterns see the normal form: 1/x is "1 x /", x^-2 is "1 x 2 ^ /", √x is "x sqrt",\n# 1/√(1 - x²) is "1 x 2 ^ neg 1 + sqrt /".\n\nrules:\n  # ---- splitting -----------------------------------------------------------\n  - id: i-sum\n    name: Sum rule\n    mode: integrate\n    pattern: u w + v 2 int\n    output: u v 2 int w v 2 int +\n    explain: "Integrate {u} and {w} separately, and add."\n    tags: [split]\n\n  - id: i-difference\n    name: Difference rule\n    mode: integrate\n    pattern: u w - v 2 int\n    output: u v 2 int w v 2 int -\n    explain: "Integrate {u} and {w} separately, and subtract."\n    tags: [split]\n\n  - id: i-negate\n    name: A minus sign\n    mode: integrate\n    pattern: u neg v 2 int\n    output: u v 2 int neg\n    explain: "The minus sign comes out of the integral."\n    tags: [split]\n\n  - id: i-constant-multiple\n    name: Constant multiple\n    mode: integrate\n    pattern: c u * v 2 int\n    commutative: true\n    require:\n      - [freeOf, c, v]\n      - [dependsOn, u, v]\n    output: c u v 2 int *\n    priority: 5\n    explain: "{c} doesn\'t change with {v}, so it comes out of the integral."\n    tags: [split]\n\n  - id: i-constant-divisor\n    name: Dividing by a constant\n    mode: integrate\n    pattern: u c / v 2 int\n    require:\n      - [freeOf, c, v]\n      - [dependsOn, u, v]\n    output: u v 2 int c /\n    priority: 5\n    explain: "Dividing by {c} comes out of the integral."\n    tags: [split]\n\n  - id: i-constant-over\n    name: A constant over something\n    mode: integrate\n    pattern: c u / v 2 int\n    require:\n      - [freeOf, c, v]\n      - [dependsOn, u, v]\n    when:\n      c: [notOne]\n    output: c 1 u / v 2 int *\n    priority: 4\n    explain: "{c} comes out of the integral, leaving ∫ 1/{u}."\n    tags: [split]\n\n  # ---- the table -----------------------------------------------------------\n  - id: i-constant\n    name: Integral of a constant\n    mode: integrate\n    pattern: c v 2 int\n    require:\n      - [freeOf, c, v]\n    output: c v *\n    priority: 10\n    explain: "A constant {c} integrates to {c}·{v}."\n    tags: [base]\n\n  - id: i-variable\n    name: Integral of the variable\n    mode: integrate\n    pattern: u v 2 int\n    require:\n      - [sameAs, u, v]\n    output: v 2 ^ 2 /\n    priority: 10\n    explain: "∫ {v} d{v} = {v}²/2."\n    tags: [base]\n\n  - id: i-power\n    name: Power rule\n    mode: integrate\n    pattern: u n ^ v 2 int\n    require:\n      - [sameAs, u, v]\n      - [freeOf, n, v]\n    when:\n      n: [notMinusOne]\n    output: v n 1 + ^ n 1 + /\n    explain: "Add 1 to the power and divide by the new power: {result}."\n    tags: [power]\n\n  - id: i-reciprocal\n    name: Integral of 1/x\n    mode: integrate\n    pattern: 1 u / v 2 int\n    require:\n      - [sameAs, u, v]\n    output: v abs ln\n    explain: "∫ 1/{v} d{v} = ln|{v}| (the power rule can\'t divide by 0 here)."\n    tags: [log]\n\n  - id: i-reciprocal-power\n    name: Power rule (a power below)\n    mode: integrate\n    pattern: 1 u n ^ / v 2 int\n    require:\n      - [sameAs, u, v]\n      - [freeOf, n, v]\n    when:\n      n: [notOne]\n    output: v 1 n - ^ 1 n - /\n    explain: "1/{v}^{n} is {v}^(-{n}): add 1 to the power and divide by it."\n    tags: [power]\n\n  - id: i-sqrt\n    name: Integral of a square root\n    mode: integrate\n    pattern: u sqrt v 2 int\n    require:\n      - [sameAs, u, v]\n    output: 2 v 3 2 / ^ * 3 /\n    explain: "√{v} is {v}^(1/2): add 1 to the power and divide by 3/2."\n    tags: [power]\n\n  - id: i-reciprocal-sqrt\n    name: Integral of 1/√x\n    mode: integrate\n    pattern: 1 u sqrt / v 2 int\n    require:\n      - [sameAs, u, v]\n    output: 2 v sqrt *\n    explain: "1/√{v} is {v}^(-1/2): add 1 to the power and divide by 1/2."\n    tags: [power]\n\n  - id: i-root\n    name: Integral of a root\n    mode: integrate\n    pattern: n u root v 2 int\n    require:\n      - [sameAs, u, v]\n      - [freeOf, n, v]\n    output: v 1 n / 1 + ^ 1 n / 1 + /\n    explain: "The {n}th root of {v} is {v}^(1/{n}): add 1 to the power and divide by it."\n    tags: [power]\n\n  - id: i-exp\n    name: Integral of e^x\n    mode: integrate\n    pattern: =e u ^ v 2 int\n    require:\n      - [sameAs, u, v]\n    output: =e v ^\n    explain: "e^{v} is its own integral."\n    tags: [exp]\n\n  - id: i-exponential\n    name: Integral of a^x\n    mode: integrate\n    pattern: a u ^ v 2 int\n    require:\n      - [sameAs, u, v]\n      - [freeOf, a, v]\n    output: a v ^ a ln /\n    explain: "∫ {a}^{v} d{v} = {a}^{v} / ln {a}."\n    tags: [exp]\n\n  - id: i-ln\n    name: Integral of ln x\n    mode: integrate\n    pattern: u ln v 2 int\n    require:\n      - [sameAs, u, v]\n    output: v v ln * v -\n    explain: "By parts: ∫ ln {v} d{v} = {v} ln {v} - {v}."\n    tags: [log]\n\n  - { id: i-sin, name: Integral of sin, mode: integrate, pattern: u sin v 2 int, require: [[sameAs, u, v]], output: v cos neg, tags: [trig], explain: "The derivative of -cos {v} is sin {v}." }\n  - { id: i-cos, name: Integral of cos, mode: integrate, pattern: u cos v 2 int, require: [[sameAs, u, v]], output: v sin, tags: [trig], explain: "The derivative of sin {v} is cos {v}." }\n  - { id: i-tan, name: Integral of tan, mode: integrate, pattern: u tan v 2 int, require: [[sameAs, u, v]], output: v cos abs ln neg, tags: [trig], explain: "tan = sin/cos, and -ln|cos {v}| differentiates to it." }\n  - { id: i-sec-squared, name: Integral of sec², mode: integrate, pattern: 1 u cos 2 ^ / v 2 int, require: [[sameAs, u, v]], output: v tan, tags: [trig], explain: "The derivative of tan {v} is 1/cos²{v}." }\n  - { id: i-sec-squared-written, name: Integral of sec², mode: integrate, pattern: u sec 2 ^ v 2 int, require: [[sameAs, u, v]], output: v tan, tags: [trig], explain: "The derivative of tan {v} is sec²{v}." }\n  - { id: i-csc-squared, name: Integral of csc², mode: integrate, pattern: 1 u sin 2 ^ / v 2 int, require: [[sameAs, u, v]], output: v cot neg, tags: [trig], explain: "The derivative of -cot {v} is 1/sin²{v}." }\n  - { id: i-sin-squared, name: Integral of sin², mode: integrate, pattern: u sin 2 ^ v 2 int, require: [[sameAs, u, v]], output: v 2 / 2 v * sin 4 / -, tags: [trig], explain: "sin²{v} = (1 - cos 2{v})/2, which integrates term by term." }\n  - { id: i-cos-squared, name: Integral of cos², mode: integrate, pattern: u cos 2 ^ v 2 int, require: [[sameAs, u, v]], output: v 2 / 2 v * sin 4 / +, tags: [trig], explain: "cos²{v} = (1 + cos 2{v})/2, which integrates term by term." }\n  - { id: i-sinh, name: Integral of sinh, mode: integrate, pattern: u sinh v 2 int, require: [[sameAs, u, v]], output: v cosh, tags: [trig], explain: "The derivative of cosh {v} is sinh {v}." }\n  - { id: i-cosh, name: Integral of cosh, mode: integrate, pattern: u cosh v 2 int, require: [[sameAs, u, v]], output: v sinh, tags: [trig], explain: "The derivative of sinh {v} is cosh {v}." }\n  - { id: i-arctan, name: Integral of 1/(1 + x²), mode: integrate, pattern: 1 u 2 ^ 1 + / v 2 int, require: [[sameAs, u, v]], output: v arctan, tags: [trig], explain: "The derivative of arctan {v} is 1/(1 + {v}²)." }\n  - { id: i-arcsin, name: Integral of 1/√(1 - x²), mode: integrate, pattern: 1 u 2 ^ neg 1 + sqrt / v 2 int, require: [[sameAs, u, v]], output: v arcsin, tags: [trig], explain: "The derivative of arcsin {v} is 1/√(1 - {v}²)." }\n','name: Values and identities\n\n# The CAS\'s simplify rules (mode: simplify). They see expressions in normal form\n# (src/cas/form.js), as drawn: x/2, 2x - 3, -x, √x is "x sqrt",\n# ∛x is "3 x root", π/6 is "=pi 6 /", 2π/3 is "2 =pi * 3 /", e^x is "=e x ^".\n#\n#   - a word is a slot (u, v, n …); a fixed symbol is written =pi, =e\n#   - a number is itself; "1 2 /" is ½\n#   - arithmetic (2 + 3, x + x, √8 = 2√2) is done for you before and after every rule\n#   - a rule for a sum or product of two parts also fires on any two parts of a longer one\n#   - tags pick what a step list shows (steps: \'brief\' hides value, identity, tidy …)\n#\n# To change one: copy it into a file of your own with the same id and edit it;\n# to turn one off: give its id with enabled: false. createCAS({ rules: [yours] }).\n\nrules:\n  # ---- logs and exponentials -------------------------------------------\n  - id: ln-one\n    name: Log of 1\n    pattern: 1 ln\n    output: "0"\n    explain: "e^0 = 1, so ln 1 = 0."\n    tags: [value]\n\n  - id: ln-e\n    name: Log of e\n    pattern: =e ln\n    output: "1"\n    explain: "e^1 = e, so ln e = 1."\n    tags: [value]\n\n  - id: ln-exp\n    name: Log undoes e to a power\n    pattern: =e u ^ ln\n    output: u\n    explain: "ln and e^ undo each other: ln(e^{u}) = {u}."\n    tags: [inverse]\n\n  - id: exp-ln\n    name: e to a power undoes the log\n    pattern: =e u ln ^\n    output: u\n    explain: "e^ and ln undo each other: e^(ln {u}) = {u} (for {u} > 0)."\n    tags: [inverse]\n\n  - id: exp-multiple-ln\n    name: A multiple of a log in a power\n    pattern: =e n u ln * ^\n    when:\n      n: [isRational]\n    output: u n ^\n    explain: "e^({n} ln {u}) = (e^(ln {u}))^{n} = {result}."\n    tags: [inverse]\n\n  - id: log-same-base\n    name: Log of its own base\n    pattern: b b log\n    output: "1"\n    explain: "{b}^1 = {b}, so the log of {b} to base {b} is 1."\n    tags: [value]\n\n  - id: log-of-power\n    name: Log undoes a power of its base\n    pattern: b b u ^ log\n    output: u\n    explain: "log base {b} and {b}^ undo each other."\n    tags: [inverse]\n\n  - id: power-of-log\n    name: A power undoes the log to its base\n    pattern: b b u log ^\n    output: u\n    explain: "{b}^(log base {b} of {u}) = {u}."\n    tags: [inverse]\n\n  # ---- roots and absolute value --------------------------------------------\n  - id: sqrt-of-square\n    name: Square root of a square\n    pattern: u 2 ^ sqrt\n    output: u abs\n    explain: "√({u}²) is the size of {u}, whatever its sign: |{u}|."\n    tags: [identity]\n\n  - id: abs-squared\n    name: Absolute value squared\n    pattern: u abs 2 ^\n    output: u 2 ^\n    explain: "Squaring removes the sign anyway: |{u}|² = {u}²."\n    tags: [identity]\n\n  - id: abs-of-abs\n    name: Absolute value twice\n    pattern: u abs abs\n    output: u abs\n    explain: "|{u}| is never negative, so its absolute value is itself."\n    tags: [identity]\n\n  - id: abs-positive\n    name: Absolute value of a positive number\n    pattern: u abs\n    when:\n      u: [isPositiveConstant]\n    output: u\n    explain: "{u} is positive, so |{u}| = {u}."\n    tags: [value]\n\n  - id: abs-negative\n    name: Absolute value of a negative number\n    pattern: u abs\n    when:\n      u: [isNegativeConstant]\n    output: u neg\n    explain: "{u} is negative, so |{u}| = -({u})."\n    tags: [value]\n\n  - id: abs-of-exp\n    name: e to a power is positive\n    pattern: =e u ^ abs\n    output: =e u ^\n    explain: "e to any power is positive, so the absolute value changes nothing."\n    tags: [identity]\n\n  # ---- trigonometry: signs and turns ---------------------------------------\n  - id: sin-odd\n    name: Sine of a negative angle\n    pattern: u sin\n    when:\n      u: [hasMinusSign]\n    output: u neg sin neg\n    explain: "sin(-θ) = -sin θ."\n    tags: [symmetry]\n\n  - id: cos-even\n    name: Cosine of a negative angle\n    pattern: u cos\n    when:\n      u: [hasMinusSign]\n    output: u neg cos\n    explain: "cos(-θ) = cos θ."\n    tags: [symmetry]\n\n  - id: tan-odd\n    name: Tangent of a negative angle\n    pattern: u tan\n    when:\n      u: [hasMinusSign]\n    output: u neg tan neg\n    explain: "tan(-θ) = -tan θ."\n    tags: [symmetry]\n\n  - id: sin-full-turn\n    name: A full turn\n    pattern: u sin\n    when:\n      u: [[piMultipleAtLeast, 2]]\n    output: u 2 =pi * - sin\n    explain: "Turning a full 2π brings you back: sin({u}) = sin({u} - 2π)."\n    tags: [symmetry]\n\n  - id: cos-full-turn\n    name: A full turn\n    pattern: u cos\n    when:\n      u: [[piMultipleAtLeast, 2]]\n    output: u 2 =pi * - cos\n    explain: "Turning a full 2π brings you back: cos({u}) = cos({u} - 2π)."\n    tags: [symmetry]\n\n  - id: tan-half-turn\n    name: Tangent repeats every π\n    pattern: u tan\n    when:\n      u: [[piMultipleAtLeast, 1]]\n    output: u =pi - tan\n    explain: "tan repeats every π: tan({u}) = tan({u} - π)."\n    tags: [symmetry]\n\n  - id: sin-half-turn\n    name: A half turn\n    pattern: u sin\n    when:\n      u: [[piMultipleAbove, 1], [piMultipleBelow, 2]]\n    output: u =pi - sin neg\n    explain: "Half a turn flips the sign: sin({u}) = -sin({u} - π)."\n    tags: [symmetry]\n\n  - id: cos-half-turn\n    name: A half turn\n    pattern: u cos\n    when:\n      u: [[piMultipleAbove, 1], [piMultipleBelow, 2]]\n    output: u =pi - cos neg\n    explain: "Half a turn flips the sign: cos({u}) = -cos({u} - π)."\n    tags: [symmetry]\n\n  - id: sin-second-quarter\n    name: The second quarter turn\n    pattern: u sin\n    when:\n      u: [[piMultipleAbove, 1/2], [piMultipleAtMost, 1]]\n    output: =pi u - sin\n    explain: "sin({u}) = sin(π - {u})."\n    tags: [symmetry]\n\n  - id: cos-second-quarter\n    name: The second quarter turn\n    pattern: u cos\n    when:\n      u: [[piMultipleAbove, 1/2], [piMultipleAtMost, 1]]\n    output: =pi u - cos neg\n    explain: "cos({u}) = -cos(π - {u})."\n    tags: [symmetry]\n\n  - id: tan-second-quarter\n    name: The second quarter turn\n    pattern: u tan\n    when:\n      u: [[piMultipleAbove, 1/2], [piMultipleBelow, 1]]\n    output: =pi u - tan neg\n    explain: "tan({u}) = -tan(π - {u})."\n    tags: [symmetry]\n\n  # ---- trigonometry: the standard angles -------------------------------------\n  - { id: sin-0, name: Sine of 0, pattern: 0 sin, output: "0", tags: [value], explain: "sin 0 = 0." }\n  - { id: sin-pi-6, name: Sine of π/6, pattern: =pi 6 / sin, output: 1 2 /, tags: [value], explain: "sin(π/6) = 1/2." }\n  - { id: sin-pi-4, name: Sine of π/4, pattern: =pi 4 / sin, output: 2 sqrt 2 /, tags: [value], explain: "sin(π/4) = √2/2." }\n  - { id: sin-pi-3, name: Sine of π/3, pattern: =pi 3 / sin, output: 3 sqrt 2 /, tags: [value], explain: "sin(π/3) = √3/2." }\n  - { id: sin-pi-2, name: Sine of π/2, pattern: =pi 2 / sin, output: "1", tags: [value], explain: "sin(π/2) = 1." }\n  - { id: sin-pi, name: Sine of π, pattern: =pi sin, output: "0", tags: [value], explain: "sin π = 0." }\n  - { id: cos-0, name: Cosine of 0, pattern: 0 cos, output: "1", tags: [value], explain: "cos 0 = 1." }\n  - { id: cos-pi-6, name: Cosine of π/6, pattern: =pi 6 / cos, output: 3 sqrt 2 /, tags: [value], explain: "cos(π/6) = √3/2." }\n  - { id: cos-pi-4, name: Cosine of π/4, pattern: =pi 4 / cos, output: 2 sqrt 2 /, tags: [value], explain: "cos(π/4) = √2/2." }\n  - { id: cos-pi-3, name: Cosine of π/3, pattern: =pi 3 / cos, output: 1 2 /, tags: [value], explain: "cos(π/3) = 1/2." }\n  - { id: cos-pi-2, name: Cosine of π/2, pattern: =pi 2 / cos, output: "0", tags: [value], explain: "cos(π/2) = 0." }\n  - { id: cos-pi, name: Cosine of π, pattern: =pi cos, output: "-1", tags: [value], explain: "cos π = -1." }\n  - { id: tan-0, name: Tangent of 0, pattern: 0 tan, output: "0", tags: [value], explain: "tan 0 = 0." }\n  - { id: tan-pi-6, name: Tangent of π/6, pattern: =pi 6 / tan, output: 3 sqrt 3 /, tags: [value], explain: "tan(π/6) = √3/3." }\n  - { id: tan-pi-4, name: Tangent of π/4, pattern: =pi 4 / tan, output: "1", tags: [value], explain: "tan(π/4) = 1." }\n  - { id: tan-pi-3, name: Tangent of π/3, pattern: =pi 3 / tan, output: 3 sqrt, tags: [value], explain: "tan(π/3) = √3." }\n  - { id: arcsin-0, name: Arcsine of 0, pattern: 0 arcsin, output: "0", tags: [value], explain: "sin 0 = 0, so arcsin 0 = 0." }\n  - { id: arcsin-half, name: Arcsine of 1/2, pattern: 1 2 / arcsin, output: =pi 6 /, tags: [value], explain: "sin(π/6) = 1/2, so arcsin(1/2) = π/6." }\n  - { id: arcsin-1, name: Arcsine of 1, pattern: 1 arcsin, output: =pi 2 /, tags: [value], explain: "sin(π/2) = 1, so arcsin 1 = π/2." }\n  - { id: arccos-0, name: Arccosine of 0, pattern: 0 arccos, output: =pi 2 /, tags: [value], explain: "cos(π/2) = 0, so arccos 0 = π/2." }\n  - { id: arccos-half, name: Arccosine of 1/2, pattern: 1 2 / arccos, output: =pi 3 /, tags: [value], explain: "cos(π/3) = 1/2, so arccos(1/2) = π/3." }\n  - { id: arccos-1, name: Arccosine of 1, pattern: 1 arccos, output: "0", tags: [value], explain: "cos 0 = 1, so arccos 1 = 0." }\n  - { id: arccos-minus-1, name: Arccosine of -1, pattern: -1 arccos, output: =pi, tags: [value], explain: "cos π = -1, so arccos(-1) = π." }\n  - { id: arctan-0, name: Arctangent of 0, pattern: 0 arctan, output: "0", tags: [value], explain: "tan 0 = 0, so arctan 0 = 0." }\n  - { id: arctan-1, name: Arctangent of 1, pattern: 1 arctan, output: =pi 4 /, tags: [value], explain: "tan(π/4) = 1, so arctan 1 = π/4." }\n  - { id: sinh-0, name: sinh of 0, pattern: 0 sinh, output: "0", tags: [value], explain: "sinh 0 = 0." }\n  - { id: cosh-0, name: cosh of 0, pattern: 0 cosh, output: "1", tags: [value], explain: "cosh 0 = 1." }\n  - { id: tanh-0, name: tanh of 0, pattern: 0 tanh, output: "0", tags: [value], explain: "tanh 0 = 0." }\n\n  - { id: arctan-sqrt3, name: Arctangent of √3, pattern: 3 sqrt arctan, output: =pi 3 /, tags: [value], explain: "tan(π/3) = √3, so arctan √3 = π/3." }\n  - { id: arctan-sqrt3-over-3, name: Arctangent of √3/3, pattern: 3 sqrt 3 / arctan, output: =pi 6 /, tags: [value], explain: "tan(π/6) = √3/3, so arctan(√3/3) = π/6." }\n\n  # atan2(y, x): the angle of the point (x, y), in (-π, π]\n  - id: atan2-right\n    name: The angle of a point to the right\n    pattern: y x atan2\n    when:\n      x: [isPositiveConstant]\n    output: y x / arctan\n    explain: "For x > 0, atan2(y, x) = arctan(y/x)."\n    tags: [value]\n  - id: atan2-left-up\n    name: The angle of a point to the left\n    pattern: y x atan2\n    when:\n      x: [isNegativeConstant]\n      y: [isPositiveConstant]\n    output: y x / arctan =pi +\n    explain: "For x < 0 and y > 0, atan2(y, x) = arctan(y/x) + π."\n    tags: [value]\n  - id: atan2-left-down\n    name: The angle of a point to the left\n    pattern: y x atan2\n    when:\n      x: [isNegativeConstant]\n      y: [isNegativeConstant]\n    output: y x / arctan =pi -\n    explain: "For x < 0 and y < 0, atan2(y, x) = arctan(y/x) - π."\n    tags: [value]\n  - { id: atan2-negative-axis, name: The angle of a point on the negative axis, pattern: 0 x atan2, when: { x: [isNegativeConstant] }, output: =pi, tags: [value], explain: "A point on the negative x-axis is at angle π." }\n  - { id: atan2-up, name: Straight up, pattern: y 0 atan2, when: { y: [isPositiveConstant] }, output: =pi 2 /, tags: [value], explain: "A point straight up is at angle π/2." }\n  - { id: atan2-down, name: Straight down, pattern: y 0 atan2, when: { y: [isNegativeConstant] }, output: =pi 2 / neg, tags: [value], explain: "A point straight down is at angle -π/2." }\n\n  # the odd inverses take the sign out; arccos(-u) = π - arccos u\n  - id: arcsin-odd\n    name: Arcsine of a negative\n    pattern: u arcsin\n    when:\n      u: [hasMinusSign]\n    output: u neg arcsin neg\n    explain: "arcsin(-x) = -arcsin x."\n    tags: [symmetry]\n\n  - id: arctan-odd\n    name: Arctangent of a negative\n    pattern: u arctan\n    when:\n      u: [hasMinusSign]\n    output: u neg arctan neg\n    explain: "arctan(-x) = -arctan x."\n    tags: [symmetry]\n\n  - id: arccos-negative\n    name: Arccosine of a negative\n    pattern: u arccos\n    when:\n      u: [hasMinusSign]\n    output: =pi u neg arccos -\n    explain: "arccos(-x) = π - arccos x."\n    tags: [symmetry]\n\n  # ---- trigonometry: identities --------------------------------------------\n  - id: pythagorean\n    name: sin² + cos² = 1\n    pattern: u sin 2 ^ u cos 2 ^ +\n    commutative: true\n    output: "1"\n    explain: "sin²{u} + cos²{u} = 1 for every angle."\n    tags: [identity]\n\n  - id: pythagorean-multiple\n    name: sin² + cos² = 1\n    pattern: a u sin 2 ^ * a u cos 2 ^ * +\n    commutative: true\n    output: a\n    explain: "{a}·sin²{u} + {a}·cos²{u} = {a}(sin²{u} + cos²{u}) = {a}."\n    tags: [identity]\n\n  - id: one-minus-sin-squared\n    name: 1 - sin² = cos²\n    pattern: u sin 2 ^ neg 1 +\n    output: u cos 2 ^\n    explain: "1 - sin²{u} = cos²{u}."\n    tags: [identity]\n\n  - id: one-minus-cos-squared\n    name: 1 - cos² = sin²\n    pattern: u cos 2 ^ neg 1 +\n    output: u sin 2 ^\n    explain: "1 - cos²{u} = sin²{u}."\n    tags: [identity]\n\n  - id: sin-over-cos\n    name: sin/cos is tan\n    pattern: u sin u cos /\n    output: u tan\n    explain: "sin {u} / cos {u} = tan {u}."\n    tags: [identity]\n\n  - id: cos-over-sin\n    name: cos/sin is cot\n    pattern: u cos u sin /\n    output: u cot\n    explain: "cos {u} / sin {u} = cot {u}."\n    tags: [identity]\n\n  - id: sin-arcsin\n    name: sin undoes arcsin\n    pattern: u arcsin sin\n    output: u\n    explain: "sin(arcsin {u}) = {u}."\n    tags: [inverse]\n\n  - id: cos-arccos\n    name: cos undoes arccos\n    pattern: u arccos cos\n    output: u\n    explain: "cos(arccos {u}) = {u}."\n    tags: [inverse]\n\n  - id: tan-arctan\n    name: tan undoes arctan\n    pattern: u arctan tan\n    output: u\n    explain: "tan(arctan {u}) = {u}."\n    tags: [inverse]\n\n  - id: cosh-squared-minus-sinh-squared\n    name: cosh² - sinh² = 1\n    pattern: u cosh 2 ^ u sinh 2 ^ -\n    output: "1"\n    explain: "cosh²{u} - sinh²{u} = 1."\n    tags: [identity]\n\n  # ---- logs: opt-in rewrites (run with cas.rewrite(x, \'expand-log\') or \'combine-log\') ----\n  - id: ln-product\n    name: Log of a product\n    mode: expand-log\n    pattern: a b * ln\n    output: a ln b ln +\n    explain: "ln({a}·{b}) = ln {a} + ln {b} (for positive {a}, {b})."\n    tags: [log]\n\n  - id: ln-quotient\n    name: Log of a quotient\n    mode: expand-log\n    pattern: a b / ln\n    output: a ln b ln -\n    explain: "ln({a}/{b}) = ln {a} - ln {b} (for positive {a}, {b})."\n    tags: [log]\n\n  - id: ln-power\n    name: Log of a power\n    mode: expand-log\n    pattern: a n ^ ln\n    output: n a ln *\n    explain: "ln({a}^{n}) = {n}·ln {a} (for positive {a})."\n    tags: [log]\n\n  - id: ln-root\n    name: Log of a square root\n    mode: expand-log\n    pattern: a sqrt ln\n    output: a ln 2 /\n    explain: "ln √{a} = ½ ln {a}."\n    tags: [log]\n\n  - id: ln-sum-combine\n    name: Sum of logs\n    mode: combine-log\n    pattern: a ln b ln +\n    commutative: true\n    output: a b * ln\n    explain: "ln {a} + ln {b} = ln({a}·{b})."\n    tags: [log]\n\n  - id: ln-difference-combine\n    name: Difference of logs\n    mode: combine-log\n    pattern: a ln b ln -\n    output: a b / ln\n    explain: "ln {a} - ln {b} = ln({a}/{b})."\n    tags: [log]\n\n  - id: ln-multiple-combine\n    name: A multiple of a log\n    mode: combine-log\n    pattern: n a ln *\n    when:\n      n: [isRational]\n    output: a n ^ ln\n    explain: "{n}·ln {a} = ln({a}^{n})."\n    tags: [log]\n','name: Solving by undoing\n\n# mode: isolate. Each rule undoes the outermost function on the side with the\n# unknown: ln u = r → u = e^r. The predicates hasVariable and freeOfVariable mean\n# "mentions / doesn\'t mention the unknown being solved for".\n#\n# The CAS does the rest itself: it moves terms without the unknown to the other\n# side, divides by a coefficient, and when the unknown appears in several places\n# it uses the polynomial methods (factoring, the quadratic formula) or splits a\n# product equal to 0. Every answer is checked in the original equation, so a\n# rule may give extra candidates (squaring both sides does): the wrong ones are\n# dropped, with a step saying so.\n#\n# An output can give several answers: `pm` (±) splits into two, and so does\n# `|` (or): "u r arcsin = u =pi r arcsin - = |".\n\nrules:\n  - id: iso-swap\n    name: Unknown on the left\n    mode: isolate\n    pattern: l r =\n    when:\n      l: [freeOfVariable]\n      r: [hasVariable]\n    output: r l =\n    priority: 20\n    explain: "Swap the sides so the unknown is on the left."\n    tags: [tidy]\n\n  - id: iso-ln\n    name: Undo the log\n    mode: isolate\n    pattern: u ln r =\n    when:\n      u: [hasVariable]\n    output: u =e r ^ =\n    explain: "Raise e to each side: e^(ln {u}) = {u}."\n    tags: [inverse]\n\n  - id: iso-log\n    name: Undo the log\n    mode: isolate\n    pattern: b u log r =\n    when:\n      u: [hasVariable]\n      b: [freeOfVariable]\n    output: u b r ^ =\n    explain: "Raise {b} to each side: {b}^(log base {b} of {u}) = {u}."\n    tags: [inverse]\n\n  - id: iso-exp\n    name: Undo e to a power\n    mode: isolate\n    pattern: =e u ^ r =\n    when:\n      u: [hasVariable]\n    output: u r ln =\n    explain: "Take ln of each side: ln(e^{u}) = {u}."\n    tags: [inverse]\n\n  - id: iso-exponential\n    name: Undo a power of a number\n    mode: isolate\n    pattern: b u ^ r =\n    when:\n      u: [hasVariable]\n      b: [freeOfVariable]\n    output: u b r log =\n    explain: "Take the log to base {b} of each side."\n    tags: [inverse]\n\n  - id: iso-square\n    name: Undo the square\n    mode: isolate\n    pattern: u 2 ^ r =\n    when:\n      u: [hasVariable]\n    output: u r sqrt pm =\n    priority: 1\n    explain: "Take the square root of each side; {u} can be positive or negative."\n    tags: [inverse]\n\n  - id: iso-even-power\n    name: Undo an even power\n    mode: isolate\n    pattern: u n ^ r =\n    when:\n      u: [hasVariable]\n      n: [isEvenWhole]\n    output: u r 1 n / ^ pm =\n    explain: "Take the root of each side (the {n}th root); {u} can be positive or negative."\n    tags: [inverse]\n\n  - id: iso-odd-power\n    name: Undo an odd power\n    mode: isolate\n    pattern: u n ^ r =\n    when:\n      u: [hasVariable]\n      n: [isOddWhole]\n    output: u r 1 n / ^ =\n    explain: "Take the root of each side (the {n}th root)."\n    tags: [inverse]\n\n  - id: iso-fractional-power\n    name: Undo a power\n    mode: isolate\n    pattern: u n ^ r =\n    when:\n      u: [hasVariable]\n      n: [isRational]\n    output: u r 1 n / ^ =\n    priority: -1\n    explain: "Raise each side to the power 1/{n}."\n    tags: [inverse]\n\n  - id: iso-sqrt\n    name: Undo the square root\n    mode: isolate\n    pattern: u sqrt r =\n    when:\n      u: [hasVariable]\n    output: u r 2 ^ =\n    explain: "Square each side (this can add answers, so each one is checked at the end)."\n    tags: [inverse]\n\n  - id: iso-root\n    name: Undo the root\n    mode: isolate\n    pattern: n u root r =\n    when:\n      u: [hasVariable]\n      n: [freeOfVariable]\n    output: u r n ^ =\n    explain: "Raise each side to the power {n}."\n    tags: [inverse]\n\n  - id: iso-reciprocal\n    name: Undo the fraction\n    mode: isolate\n    pattern: c u / r =\n    when:\n      u: [hasVariable]\n      c: [freeOfVariable]\n    output: u c r / =\n    explain: "{c}/{u} = {r}, so {u} = {c}/{r}."\n    tags: [inverse]\n\n  - id: iso-abs\n    name: Undo the absolute value\n    mode: isolate\n    pattern: u abs r =\n    when:\n      u: [hasVariable]\n    output: u r pm =\n    explain: "|{u}| = {r} means {u} is {r} or -{r}."\n    tags: [inverse]\n\n  - id: iso-sin\n    name: Undo the sine\n    mode: isolate\n    pattern: u sin r =\n    when:\n      u: [hasVariable]\n    output: u r arcsin = u =pi r arcsin - = |\n    explain: "Two angles in a turn have sine {r}: arcsin {r} and π - arcsin {r} (add any whole number of turns, 2πk, for the rest)."\n    tags: [inverse, trig]\n\n  - id: iso-cos\n    name: Undo the cosine\n    mode: isolate\n    pattern: u cos r =\n    when:\n      u: [hasVariable]\n    output: u r arccos pm =\n    explain: "Two angles in a turn have cosine {r}: ±arccos {r} (add 2πk for the rest)."\n    tags: [inverse, trig]\n\n  - id: iso-tan\n    name: Undo the tangent\n    mode: isolate\n    pattern: u tan r =\n    when:\n      u: [hasVariable]\n    output: u r arctan =\n    explain: "Take arctan of each side (add πk for the rest)."\n    tags: [inverse, trig]\n\n  - { id: iso-arcsin, name: Undo arcsin, mode: isolate, pattern: u arcsin r =, when: { u: [hasVariable] }, output: u r sin =, tags: [inverse], explain: "Take sin of each side." }\n  - { id: iso-arccos, name: Undo arccos, mode: isolate, pattern: u arccos r =, when: { u: [hasVariable] }, output: u r cos =, tags: [inverse], explain: "Take cos of each side." }\n  - { id: iso-arctan, name: Undo arctan, mode: isolate, pattern: u arctan r =, when: { u: [hasVariable] }, output: u r tan =, tags: [inverse], explain: "Take tan of each side." }\n  - { id: iso-sinh, name: Undo sinh, mode: isolate, pattern: u sinh r =, when: { u: [hasVariable] }, output: u r arcsinh =, tags: [inverse], explain: "Take arcsinh of each side." }\n  - { id: iso-tanh, name: Undo tanh, mode: isolate, pattern: u tanh r =, when: { u: [hasVariable] }, output: u r arctanh =, tags: [inverse], explain: "Take arctanh of each side." }\n  - { id: iso-cosh, name: Undo cosh, mode: isolate, pattern: u cosh r =, when: { u: [hasVariable] }, output: u r arccosh pm =, tags: [inverse], explain: "Two numbers have cosh {r}: ±arccosh {r}." }\n  - { id: iso-arcsinh, name: Undo arcsinh, mode: isolate, pattern: u arcsinh r =, when: { u: [hasVariable] }, output: u r sinh =, tags: [inverse], explain: "Take sinh of each side." }\n  - { id: iso-arccosh, name: Undo arccosh, mode: isolate, pattern: u arccosh r =, when: { u: [hasVariable] }, output: u r cosh =, tags: [inverse], explain: "Take cosh of each side." }\n  - { id: iso-arctanh, name: Undo arctanh, mode: isolate, pattern: u arctanh r =, when: { u: [hasVariable] }, output: u r tanh =, tags: [inverse], explain: "Take tanh of each side." }\n'],Lr=['name: Complex numbers\n\n# Loaded only by createCAS({ complex: true }), where i is the imaginary unit.\n# The arithmetic (i² = -1, √-4 = 2i, (2 + 3i)(1 - i) = 5 + i, 1/(1 + i) = (1 - i)/2)\n# and Re, Im, conj, arg, |z| are the normal form\'s (src/cas/form.js); unknowns\n# count as real, so Re(x + 3i) = x. Write a complex unknown as a + bi.\n\noperators:\n  - { symbol: re, arity: 1, label: real part }\n  - { symbol: im, arity: 1, label: imaginary part }\n  - { symbol: conj, arity: 1, label: complex conjugate }\n  - { symbol: arg, arity: 1, label: argument }\n\nrules:\n  - id: euler\n    name: Euler\'s formula\n    pattern: =e u ^\n    when:\n      u: [hasImaginaryPart, isConstant]\n    output: =e u re ^ u im cos =i u im sin * + *\n    explain: "e^(a + bi) = e^a (cos b + i sin b)."\n    tags: [complex]\n\n  - id: ln-negative\n    name: Log of a negative number\n    pattern: u ln\n    when:\n      u: [isNegativeConstant]\n    output: u neg ln =i =pi * +\n    explain: "{u} = |{u}|·e^(iπ), so ln {u} = ln |{u}| + iπ (the principal value)."\n    tags: [complex]\n\n  - id: ln-complex\n    name: Log of a complex number\n    pattern: u ln\n    when:\n      u: [hasImaginaryPart, isConstant]\n    output: u abs ln u arg =i * +\n    explain: "ln z = ln |z| + i arg z (the principal value)."\n    tags: [complex]\n\n  - id: i-to-the-i\n    name: i to the power i\n    pattern: =i =i ^\n    output: =e =pi 2 / neg ^\n    explain: "i = e^(iπ/2), so i^i = e^(i·iπ/2) = e^(-π/2) (the principal value)."\n    tags: [complex]\n'],_r={normal:"Work out the numbers and collect like terms.","simplify.cancel":"Cancel the factors the top and bottom share.","simplify.expand":"Multiplying out gives something shorter.","simplify.factor":"Factoring gives something shorter.",expand:"Multiply out the brackets.",factor:"Write it as a product of factors.","factor.common":"Take out the common factor {common}.","factor.polynomial":"Factor what is left: {roots}.",together:"Write it as one fraction over a common denominator.",cancel:"Cancel the factors the top and bottom share: {cancelled}.","solve.move":"Keep the terms with {v} on the left and move the others to the right.","solve.divide":"Divide both sides by {by}.","solve.multiply":"Multiply both sides by {by}.","solve.one-side":"Move everything to one side, so the equation reads … = 0.","solve.factor":"Factor the left side.","solve.zero-product":"A product is 0 exactly when one of its factors is 0.","solve.split":"There are two cases; solve each.","solve.quadratic":"The quadratic formula with a = {a}, b = {b}, c = {c}: {v} = (-b ± √(b² - 4ac)) / 2a, and b² - 4ac = {D}.","solve.no-real":"b² - 4ac = {D} is negative, so there are no real solutions.","solve.quadratic-complex":"The quadratic formula with a = {a}, b = {b}, c = {c}: {v} = (-b ± √(b² - 4ac)) / 2a. b² - 4ac = {D} is negative, so its square root is imaginary: √({D}) = {root}.","solve.roots":"{v}^{n} = {r} has {n} solutions, spaced evenly round a circle: |{r}|^(1/{n})·(cos θ + i sin θ), with θ = (arg {r} + 2πk)/{n} for k = 0 … {last}.","solve.clear-denominator":"Multiply both sides by {den} (an answer that makes it 0 is not allowed).","solve.substitute":"Write {u} for {kernel}: this is a polynomial in {u}.","solve.reject":"{value} is not a solution: {reason}.","solve.always":"This is true for every {v}.","solve.never":"This is never true, so there is no solution.","solve.unsolved":"There is no exact method here for solving for {v}; try { approximate: true }.","solve.numeric":"Found {count} solution(s) numerically.","system.eliminate":"Use one equation to remove {v} from the others.","integrate.expand":"Multiply out, then integrate term by term.","integrate.substitute":"Let {u} = {inside}, so d{u} = {a} d{v}.","integrate.back":"Put {inside} back for {u}.","integrate.parts":"By parts with u = {u} and dv = {dv} d{v}: ∫ u dv = u·v - ∫ v du, where v = {V} and du = {du} d{v}.","integrate.constant":"The number {c} comes out of the integral.","integrate.divide":"The top has at least the degree of the bottom: divide first.","integrate.partial-fractions":"Split the fraction into simple fractions, one for each factor below.","integrate.evaluate":"Evaluate {F} at {b} and at {a}, and subtract.",substitute:"Put in the values."},Hr=["normal","rewrite","simplify","expand","factor","cancel","diff","differentiate","integrate","solve","solveSystem","substitute","approximate","equivalent"];class Jr{constructor(e,n,t,r,s){this.set=e,this.messages=n,this.format=t,this.read=r,this.complex=s;for(const e of Hr){const n=Jr.prototype[e];this[e]=(...e)=>ae(s,()=>n.apply(this,e))}}node(e){if(e&&"object"==typeof e&&!Array.isArray(e)&&e.type)return e;if(Array.isArray(e))return U(e);if("string"!=typeof e)throw new Q("expected an expression");if(this.read){const n=this.read(e);return Array.isArray(n)?U(n):"string"==typeof n?R(n):n}return R(e)}step(e,n,t,r,s,o={}){return function(e,n,t,r,s,o,i,a,u=[]){const l=(o??"").replace(/\{(\w+)\}/g,(e,n)=>{const t=i?.[n];return void 0===t?e:"string"==typeof t?t:a(t)});return{id:e,name:n,tags:t,from:r.s,to:s.s,path:u,before:r.s,after:s.s,explanation:l}}(e,n,t,r,s,this.messages[e]??"",o,this.format)}result(e,n,t,r={}){let s=e;if(t.approximate&&s){s=$t(s,{digits:"number"==typeof t.approximate?t.approximate:15})}return{result:s,text:s?this.format(s):"",steps:Rt(n,Vt(t.steps)),...r}}variable(e,n){if(n)return n;const t=[...Ke(e)];return 1===t.length?t[0]:t.includes("x")?"x":t[0]??"x"}tidyStep(e,n){return e.s===n.s?[]:[this.step("normal","Work out and collect",["tidy","arithmetic"],e,n)]}normal(e){return Le(this.node(e))}rewrite(e,n,t={}){const r=this.node(e),s=this.set.view(n);if(!s.size)throw new Q(`no rules in mode "${n}" (modes: ${this.set.modes().join(", ")})`);const o=Bt(r,s,{format:this.format,order:t.order,variable:t.variable});return this.result(o.node,[...this.tidyStep(r,Le(r)),...o.steps],t)}simplify(e,n={}){const t=this.node(e),r=[...this.tidyStep(t,Le(t))],s=Bt(t,this.set.view("simplify"),{format:this.format});r.push(...s.steps);let o=s.node,i=null;const a=We(o),u=[],l=0===Ke(o).size,c=l?null:cr(a);c&&u.push({node:ar(c),id:"simplify.cancel",name:"Cancel"});const f=Bt(Qe(Wt(a)),this.set.view("simplify"),{format:this.format});u.push({node:f.node,id:"simplify.expand",name:"Multiply out"}),"add"!==a.k||l||u.push({node:ar(nr(a)),id:"simplify.factor",name:"Factor"});for(const e of u)W(e.node)<W(o)&&(i=this.step(e.id,e.name,["shorter"],o,e.node),o=e.node);return i&&r.push(i),this.result(o,r,n)}expand(e,n={}){const t=this.node(e),r=Qe(Wt(We(t)));return this.result(r,t.s===r.s?[]:[this.step("expand","Multiply out",["expand"],t,r)],n)}factor(e,n={}){const t=this.node(e),r=nr(We(t)),s=ar(r),o=[],i=r.notes.find(e=>"factor.common"===e.id);return s.s!==Le(t).s&&o.push(this.step(i?"factor.common":"factor","Factor",["factor"],t,s,i?{common:Qe(i.values.common)}:{})),this.result(s,o,n)}cancel(e,n={}){const t=this.node(e),r=We(t),s=cr(r);if(!s)return this.result(Le(t),this.tidyStep(t,Le(t)),n);const o=ar(s),i=lr(r),a=C("/",[Qe(i.num),Qe(i.den)]),u=[];return"/"===t.op&&a.s===Le(t).s||u.push(this.step("together","One fraction",["together"],t,a)),o.s!==a.s&&u.push(this.step("cancel","Cancel",["cancel"],a,o,{cancelled:s.cancelled.length?s.cancelled.map(e=>this.format(Qe(e))).join(", "):"numbers"})),this.result(o,u,n)}diff(e,n,t={}){const r=this.node(e),s=this.variable(r,n);let o=r;const i=[];for(let e=0;e<(t.order??1);e++){const e=Bt(C("der",[o,Z(s)]),this.set.view("differentiate"),{format:this.format,order:"pre"}),n=Bt(e.node,this.set.view("simplify"),{format:this.format});i.push(...e.steps,...n.steps),o=n.node}return this.result(o,i,t,{variable:s})}differentiate(e,n){const t=Bt(C("der",[e,Z(n)]),this.set.view("differentiate"),{format:this.format,order:"pre"});return Bt(t.node,this.set.view("simplify"),{format:this.format}).node}integrate(e,n,t={}){let r=this.node(e);const s=e=>void 0===e?null:this.node("number"==typeof e?String(e):e);let o=s(t.from),i=s(t.to);if(V(r)&&"int"===r.op){const[e,t,s,a]=r.args;n=t.token,s&&a&&(o=s,i=a),r=e}const a=this.variable(r,n),u=Er(this,r,a),l=Bt(u.node,this.set.view("simplify"),{format:this.format}),c=[...u.steps,...l.steps];if(o&&i){const e=zr(this,l.node,a,o,i,r);c.push(e.step);const n=Bt(e.value,this.set.view("simplify"),{format:this.format});return c.push(...n.steps),this.result(n.node,c,t,{antiderivative:l.node,variable:a})}let f=l.node;return t.constant&&(f=C("+",[f,Z(String(t.constant))])),this.result(f,c,t,{variable:a})}solve(e,n,t={}){const r=this.node(e),s=this.variable(r,n),o=hr(this,r,s,t);let i=o.solutions;if(t.approximate){const e="number"==typeof t.approximate?t.approximate:15;i=i.map(n=>$t(n,{digits:e}))}const a=i.map(e=>C("=",[Z(s),e])),u=a.length?a.reduce((e,n)=>e?C("|",[e,n]):n,null):null;return{solutions:i,result:u,variable:s,text:o.always?`true for every ${s}`:u?this.format(u):o.unsolved.length?`no exact solution found for ${s}`:"no solution",unsolved:o.unsolved.filter(e=>!e.always),always:o.always,steps:Rt(o.steps,Vt(t.steps))}}solveSystem(e,n,t={}){const r=e.map(e=>this.node(e)),s=n??[...new Set(r.flatMap(e=>[...Ke(e)]))].sort(),o=function(e,n,t){const r=[],s=n.map(n=>{const r=dr(n)?Se([We(n.args[0]),Ce(We(n.args[1]))]):We(n),s=t.map(n=>{const t=_t(r,n);if(!t||Math.max(0,...t.keys())>1)throw new Q(`not linear in ${n}`);const s=t.get(1)??ce(ee.ZERO);if(!he(s))throw new Q(`the coefficient of ${n} must be a number (it is ${e.format(Qe(s))})`);return s.q});let o=r;for(let e=0;e<t.length;e++)o=Se([o,Ne([ce(s[e].neg()),We(Z(t[e]))])]);for(const e of t)if(Ht(o,e))throw new Q(`not linear in ${e}`);return{coeffs:s,rhs:Ce(o)}}),o=()=>s.map(e=>fr(Qe(Se(e.coeffs.map((e,n)=>Ne([ce(e),We(Z(t[n]))])))),Qe(e.rhs))).reduce((e,n)=>e?C("&",[e,n]):n,null);let i=0;const a=[];for(let n=0;n<t.length&&i<s.length;n++){const u=s.findIndex((e,t)=>t>=i&&!e.coeffs[n].isZero);if(u<0)continue;[s[i],s[u]]=[s[u],s[i]];const l=o(),c=s[i],f=c.coeffs[n];c.coeffs=c.coeffs.map(e=>e.div(f)),c.rhs=Ne([ce(f.inv()),c.rhs]);for(let e=0;e<s.length;e++){if(e===i||s[e].coeffs[n].isZero)continue;const t=s[e].coeffs[n];s[e].coeffs=s[e].coeffs.map((e,n)=>e.sub(t.mul(c.coeffs[n]))),s[e].rhs=Se([s[e].rhs,Ne([ce(t.neg()),c.rhs])])}r.push(e.step("system.eliminate","Eliminate",["eliminate"],l,o(),{v:Z(t[n])})),a.push(n),i++}for(const e of s.slice(i))if(!he(e.rhs)||!e.rhs.q.isZero)return{solution:null,steps:r,status:"none"};if(a.length<t.length)return{solution:null,steps:r,status:"many"};const u=new Map;return a.forEach((e,n)=>u.set(t[e],Qe(s[n].rhs))),{solution:u,steps:r,status:"one"}}(this,r,s),i=o.solution?Object.fromEntries([...o.solution].map(([e,n])=>[e,n])):null,a="one"===o.status?s.map(e=>`${e} = ${this.format(i[e])}`).join(", "):"none"===o.status?"no solution":"infinitely many solutions";return{solution:i,status:o.status,text:a,steps:Rt(o.steps,Vt(t.steps))}}substitute(e,n,t={}){const r=this.node(e);let s=r;for(const[e,t]of Object.entries(n))s=br(s,e,"number"==typeof t?this.node(String(t)):this.node(t));const o=this.simplify(s,t),i=[this.step("substitute","Substitute",["substitute"],r,s)];return{...o,steps:[...Rt(i,Vt(t.steps)),...o.steps]}}defined(e){const n=this.set.view("define");return n.size?Bt(e,n,{format:this.format}).node:e}approximate(e,n={}){const t=this.defined(this.node(e));if(Ke(t).size&&!n.scope){const e=$t(Le(t),n);return{result:e,text:this.format(e),number:null,exact:!1}}const r=xt(t,n);return{result:r.node,text:r.text,number:r.number,re:r.re,im:r.im,exact:r.exact}}equivalent(e,n){const t=this.node(e),r=this.node(n),s=this.defined(C("-",[t,r])),o=this.simplify(s,{steps:!1}).result;if("num"===o.type&&"0"===o.token)return{equal:!0,how:"exact"};const i=[...Ke(s)];if(!i.length)try{return{equal:Math.abs(xt(s,{digits:30}).number)<1e-25,how:"numeric"}}catch{return{equal:!1,how:"numeric"}}const a=[.31,1.7,2.9,-.83,.57];let u=0;for(let e=0;e<a.length;e++){const n=Object.fromEntries(i.map((n,t)=>[n,String(a[(e+t)%a.length]+t/7)]));let t;try{t=xt(s,{digits:30,scope:n}).number}catch{continue}if(u++,Math.abs(t)>1e-20)return{equal:!1,how:"numeric"}}return{equal:u>0,how:"numeric"}}modes(){return this.set.modes()}}function Qr(e={}){const n=!1===e.defaults?[]:[...Wr,...e.complex?Lr:[]],{specs:t,messages:r}=function(e){Et();const n=[],t={};for(const r of e){const e="string"==typeof r?Tn(r):r,s=Array.isArray(e)?e:e?.rules??[];for(const n of(Array.isArray(e)?[]:e?.operators)??[]){if(!n?.symbol)throw new Error("an operator needs a symbol");S(n.symbol)||q({symbol:String(n.symbol),arity:"any"===n.arity?y:Number(n.arity),label:n.label})}Object.assign(t,(Array.isArray(e)?null:e?.messages)??{});for(const e of s)n.push(e)}return{specs:n,messages:t}}([...n,...e.rules??[]]),s=ae(e.complex,()=>Tt(t));return new Jr(s,{..._r,...r},e.format??Cr,e.read??null,Boolean(e.complex))}function Yr(p){const d=p=>e("op"===p.type?function(e,o){const{op:p,args:d}=e,[h,m]=d,g=(e,n)=>i(p,a(e),n)?`(${o(e)})`:o(e),v=()=>d.map(o).join(", "),b=e=>f(p,e)?`(${o(e)})`:o(e);if("*"===p&&n(e))return"-"+o(m);if("+"===p&&n(m))return`${g(h,"left")} ${o(m)}`;if("-"===p&&n(m))return`${g(h,"left")} + ${o(m).slice(1)}`;if("^"===p){const e=g(h,"left");return`${e.startsWith("-")?`(${e})`:e}^${g(m,"right")}`}if(Gr[p])return`${g(h,"left")} ${Gr[p]} ${g(m,"right")}`;if(Kr[p])return Kr[p][0]+v()+Kr[p][1];if(t(p))return`${r(p)}(${v()})`;switch(p){case"neg":return"-"+b(h);case"pmod":return`${o(h)} (mod ${o(m)})`;case"pm":return"±"+b(h);case"!":return b(h)+"!";case"\\not":return"¬"+b(h);case"sqrt":return`√(${o(h)})`;case"root":{const e="number"===h.type?{3:"∛",4:"∜"}[h.value]:null;return e?`${e}(${o(m)})`:`root(${o(h)}, ${o(m)})`}case"abs":return`|${o(h)}|`;case"floor":return`⌊${o(h)}⌋`;case"ceil":return`⌈${o(h)}⌉`;case"log":return`log_${"op"===h.type?`(${o(h)})`:o(h)}(${o(m)})`;case"loglog":return`ln(ln(${o(h)}))`;case"logloglog":return`ln(ln(ln(${o(h)})))`;case"perm":return`perm(${o(h)}, ${o(m)})`;case"comb":return`comb(${o(h)}, ${o(m)})`;case"_":return`${o(h)}_${o(m)}`;case"?:":return`${o(h)} ? ${o(m)} : ${o(d[2])}`;case"sum":return d.map(o).join(" + ");case"prod":return d.map(o).join(" * ");case"der":{const{n:n,body:t,by:r,variable:i}=s(e),a=n>1?`^${n}`:"";return i?`d${a}${o(t)}/d${o(r)}${a}`:`d${a}/d${o(r)}${a}(${o(t)})`}case"int":return function(e,n){const[t,r,s,o]=e,i=e.length>=2&&l(r)?n(r):c(t)||"x",a=`${n(t)} d${i}`;switch(e.length){case 1:return`∫ ${a}`;case 2:return l(r)?`∫ ${a}`:`∫[${n(r)}] ${a}`;case 3:return l(r)?`∫[${n(s)}] ${a}`:`∫[${n(r)}, ${n(s)}] ${a}`;default:return`∫[${n(s)}, ${n(o)}] ${a}`}}(d,o);case"lim":return function(e,n){const[t,r,s,o]=e;switch(e.length){case 1:return`lim ${n(t)}`;case 2:return`lim[${n(r)}] ${n(t)}`;case 3:return`lim[${n(r)}→${n(s)}] ${n(t)}`;default:{const e=u(o).toLowerCase().startsWith("r")?"⁺":"⁻";return`lim[${n(r)}→${n(s)}${e}] ${n(t)}`}}}(d,o);case"\\sumidx":return Xr("∑",d,o);case"\\prodidx":return Xr("∏",d,o)}return`${p}(${v()})`}(p,d):function(e){const n=e.value;return"symbol"!==e.type?n:o[n]?o[n]:"\\+"===n||"\\-"===n?n[1]:n}(p));return d(p)}const Gr={"+":"+","-":"-","*":"*","/":"/","^":"^",mod:"mod","±":"±","=":"=","===":"===","≠":"≠","≠=":"≠","<":"<",">":">","≤":"≤","≥":"≥","∈":"∈","∉":"∉","→":"→","≈":"≈","∪":"∪","∩":"∩","≡":"≡","&":"∧","|":"∨","\\xor":"⊕","\\nand":"↑","\\nor":"↓"},Kr={"()":["(",")"],"[]":["[","]"],"{}":["{","}"],"(]":["(","]"],"[)":["[",")"]};function Xr(e,n,t){const[r,s,o]=n,i="∑"===e?"sum":"prod";return 3===n.length&&"op"===s.type&&"="===s.op?`${i}(${t(r)}, ${t(s)}, ${t(o)})`:2===n.length&&"op"===s.type&&"∈"===s.op?`${i}(${t(r)}, ${t(s)})`:1===n.length?`${e} ${t(r)}`:2===n.length?`${e}[${t(s)}] ${t(r)}`:`${e}[${t(s)}, ${t(o)}] ${t(r)}`}for(const e of[Gr,Kr])Object.setPrototypeOf(e,null);const es=new Map,ns={re:"Re",im:"Im",conj:"overline"},ts=["Re","Im","re","im","conj","arg"];function rs(e){const n=[],t=e=>{"op"===e.type?(e.args.forEach(t),m(e.op)?(-1===m(e.op).arity&&n.push(String(e.args.length)),n.push(e.op)):n.push(String(e.args.length),`\\:${ns[e.op]??e.op}`)):n.push(e.token)};return t(e),n}function ss(e={}){const n=e.rules||e.names?null:e.complex?"complex":"default";if(n&&es.has(n))return es.get(n);const t=[...e.names??[],...e.complex?ts:[]],r=Qr({...e,read:e=>h(e,t.length?{names:t}:{}),format:Cr});return n&&es.set(n,r),r}const os=new Set(["simplify","expand","factor","cancel","diff","integrate","solve","approx"]);function is(e={}){const n=e.cas??ss(),t=(e,r)=>{let s;try{s=n.read(e)}catch{return g(e,r)}const o=s[s.length-1],i=o?.startsWith("\\:")?o.slice(2):null,a=i&&os.has(i);if(!a&&!("der"===o||"int"===o))return g(e,r);const u=n.node(s),l=a?u.args:[u],c=e=>e?e.s:void 0;let f;try{switch(a?i:"der"===o?"derivative":"antiderivative"){case"simplify":f=n.simplify(l[0]);break;case"expand":f=n.expand(l[0]);break;case"factor":f=n.factor(l[0]);break;case"cancel":f=n.cancel(l[0]);break;case"diff":f=n.diff(l[0],c(l[1]),{order:l[2]?Number(l[2].s):1});break;case"integrate":f=n.integrate(l[0],c(l[1]),l[3]?{from:l[2].s,to:l[3].s}:{});break;case"solve":f=n.solve(l[0],c(l[1]));break;case"approx":f=n.approximate(l[0],{digits:l[1]?Number(l[1].s):15});break;case"derivative":f=n.diff(u.args[0],u.args[1].s);break;default:f=n.integrate(u)}}catch(n){return`(${e})`}if(t.last=f,!f.result)return`"${f.text}"`;const h=(m=rs(f.result),p(m)?"":Yr(d(m)));var m;return/^-?\d+(\.\d+)?$/.test(h)?h:`(${h})`};return t}const as=e=>String(e).replace(/&/g,"&amp;").replace(/</g,"&lt;").replace(/>/g,"&gt;").replace(/"/g,"&quot;");function us(e,n={}){return`<ol class="mv-steps">${e.map(e=>{const t=rs(R(e.from)),r=v(t,{...n,highlight:e.path&&e.path.length?{path:e.path}:void 0}),s=e.to===e.from?"":`<span class="mv-step-arrow" aria-label="becomes">⟶</span>${v(rs(R(e.to)),n)}`,o=e.explanation?`<p class="mv-step-why">${as(e.explanation)}</p>`:"";return`<li class="mv-step" data-rule="${as(e.id)}"><span class="mv-step-name">${as(e.name)}</span><div class="mv-step-maths">${r}${s}</div>${o}</li>`}).join("")}</ol>`}export{is as casBlocks,ss as mathCAS,us as stepsHTML,rs as tokensOf};

package/dist/chunks/calculate-BX9lJ7Ff.js

@@ -0,0 +1,2 @@
+/* Inline Math 0.1.0 — MIT licence; includes mathzetta (and @mathzetta/cas, for the CAS) */
+import{f as t,g as n,h as e,i as r,j as o,k as i}from"./lite-Bj_L6w2g.js";const a={pi:Math.PI,"π":Math.PI,e:Math.E,tau:2*Math.PI,"τ":2*Math.PI,"∞":1/0,"-∞":-1/0,true:!0,false:!1},s=(t,n)=>n?s(n,t%n):Math.abs(t),u={neg:t=>-t,pm:t=>[t,-t],"!":t=>{if(!Number.isInteger(t)||t<0)return NaN;if(t>170)return 1/0;let n=1;for(let e=2;e<=t;e++)n*=e;return n},"\\not":t=>!t,sqrt:Math.sqrt,abs:Math.abs,floor:Math.floor,ceil:Math.ceil,ln:Math.log,loglog:t=>Math.log(Math.log(t)),logloglog:t=>Math.log(Math.log(Math.log(t))),expm1:Math.expm1,sin:Math.sin,cos:Math.cos,tan:Math.tan,csc:t=>1/Math.sin(t),sec:t=>1/Math.cos(t),cot:t=>1/Math.tan(t),sinh:Math.sinh,cosh:Math.cosh,tanh:Math.tanh,arcsin:Math.asin,arccos:Math.acos,arctan:Math.atan,arcsinh:Math.asinh,arccosh:Math.acosh,arctanh:Math.atanh},c={"+":(t,n)=>t+n,"-":(t,n)=>t-n,"*":(t,n)=>t*n,"/":(t,n)=>t/n,"^":(t,n)=>t**n,mod:(t,n)=>t%n,log:(t,n)=>Math.log(n)/Math.log(t),root:(t,n)=>n**(1/t),atan2:Math.atan2,round:(t,n)=>Math.round(t*10**n)/10**n,"±":(t,n)=>[t+n,t-n],fixed:(t,n)=>Number(t).toFixed(n),perm:(t,n)=>{if(!Number.isInteger(t)||!Number.isInteger(n)||n<0||n>t)return NaN;let e=1;for(let r=t-n+1;r<=t&&e!==1/0;r++)e*=r;return e},comb:(t,n)=>{if(!Number.isInteger(t)||!Number.isInteger(n)||n<0||n>t)return NaN;n=Math.min(n,t-n);let e=1;for(let r=1;r<=n&&e!==1/0;r++)e=e*(t-n+r)/r;return Math.round(e)},"=":(t,n)=>t==n,"===":(t,n)=>t===n,"≠":(t,n)=>t!=n,"≠=":(t,n)=>t!==n,"≈":(t,n)=>Math.abs(t-n)<1e-9,"≡":(t,n)=>t===n,"∪":(t,n)=>[...new Set([...[].concat(t),...[].concat(n)])],"∩":(t,n)=>[].concat(t).filter(t=>[].concat(n).includes(t)),"<":(t,n)=>t<n,">":(t,n)=>t>n,"≤":(t,n)=>t<=n,"≥":(t,n)=>t>=n,"&":(t,n)=>t&&n,"|":(t,n)=>t||n,"\\xor":(t,n)=>Boolean(t)!==Boolean(n),"\\nor":(t,n)=>!(t||n),"\\nand":(t,n)=>!(t&&n)},l={max:t=>Math.max(...t),min:t=>Math.min(...t),sum:t=>t.reduce((t,n)=>t+n,0),prod:t=>t.reduce((t,n)=>t*n,1),mean:t=>t.reduce((t,n)=>t+n,0)/t.length,gcd:t=>t.reduce(s),"()":t=>t,"[]":t=>t,"{}":t=>t};const f=5e6;let h=0,m=0;function g(t,e){if(++h>f)throw new Error("Too much to work out");switch(t.type){case"number":return Number(t.value);case"string":return t.value.slice(1,-1);case"symbol":if(Object.prototype.hasOwnProperty.call(e,t.value))return e[t.value];if(Object.prototype.hasOwnProperty.call(a,t.value))return a[t.value];throw new Error(`No value for '${t.value}'`)}const r=n(t);if(r)return g(r,e);const{op:o,args:i}=t;if("?:"===o)return g(i[0],e)?g(i[1],e):g(i[2],e);if("∈"===o||"∉"===o)return function(t,n){const[e,r]=t.args,o=g(e,n);let i;if("symbol"===r.type&&b[r.value])return i=b[r.value](o),"∈"===t.op?i:!i;if("op"===r.type&&["()","[]","(]","[)"].includes(r.op)&&2===r.args.length){const[t,e]=r.args.map(t=>g(t,n)),a="["===r.op[0]?o>=t:o>t,s="]"===r.op[1]?o<=e:o<e;i=a&&s}else{const t=g(r,n);i=Array.isArray(t)&&t.includes(o)}return"∈"===t.op?i:!i}(t,e);if("random"===o)return Math.random();if(("\\sumidx"===o||"\\prodidx"===o)&&3===i.length&&"op"===i[1].type&&"="===i[1].op&&"symbol"===i[1].args[0].type){const[t,n,r]=[i[0],i[1].args[1],i[2]],a=i[1].args[0].value,s=g(n,e),u=g(r,e);if(!Number.isFinite(s)||!Number.isFinite(u)||u-s>1e6)throw new Error(`Can't evaluate '${o}' over that range`);let c="\\sumidx"===o?0:1;for(let n=s;n<=u;n++){const r=g(t,{...e,[a]:n});c="\\sumidx"===o?c+r:c*r}return c}if("pmod"===o){const[t,n]=i,r=g(n,e);if("op"===t.type&&["≡","="].includes(t.op)){const[n,o]=t.args.map(t=>g(t,e));return(n-o)%r===0}return g(t,e)%r}const s=i.map(t=>g(t,e));if(u[o]&&1===s.length)return u[o](s[0]);if(c[o]&&2===s.length)return c[o](s[0],s[1]);if(l[o])return l[o](s);throw new Error(`Can't evaluate '${o}'`)}const b={"ℝ":Number.isFinite,"ℤ":Number.isInteger,"ℕ":t=>Number.isInteger(t)&&t>0,"ℚ":Number.isFinite,"ℂ":t=>"number"==typeof t&&!Number.isNaN(t)};for(const t of[a,u,c,l,b])Object.setPrototypeOf(t,null);class d extends Error{}class p extends Error{}const M={bits:2e5,factorial:2e4,terms:1e5,rootIndex:12,steps:5e6};let w=0,y=0;const N=t=>t<0n?-t:t,v=t=>0n===t?0:N(t).toString(2).length;function x(t,n){for(t=N(t),n=N(n);n;)[t,n]=[n,t%n];return t}function I(t,n=1n){if(0n===n)throw new d("division by zero");n<0n&&(t=-t,n=-n);const e=x(t,n)||1n;return{n:t/e,d:n/e}}const $=I(0n),B=I(1n),F=(t,n)=>I(t.n*n.d+n.n*t.d,t.d*n.d),A=(t,n)=>I(t.n*n.n,t.d*n.d),O=(t,n)=>I(t.n*n.d,t.d*n.n),S=t=>({n:-t.n,d:t.d}),P=(t,n)=>{const e=t.n*n.d-n.n*t.d;return e<0n?-1:e>0n?1:0},j=t=>{const n=t.n/t.d;return t.n<0n&&n*t.d!==t.n?n-1n:n};function k(t){const n=v(t.n),e=v(t.d);if(n<1e3&&e<1e3)return Number(t.n)/Number(t.d);const r=BigInt(Math.max(n,e)-1e3),o=t.d>>r;return 0n===o?t.n<0n?-1/0:1/0:Number(t.n>>r)/Number(o)}function E(t){let n=String(t).trim(),e=1n;n.startsWith("-")&&(e=-1n,n=n.slice(1));const r=n.endsWith("..");r&&(n=n.slice(0,-2));let o=0;const i=n.search(/e/i);i>=0&&(o=Number(n.slice(i+1)),n=n.slice(0,i));const[a="",s=""]=n.split(".");if(!(/^\d*$/.test(a)&&/^\d*$/.test(s)&&a+s))throw new p(`'${t}' is not a number`);let u;if(r&&s){let t=1;for(let n=1;2*n<=s.length;n++)if(s.slice(-n)===s.slice(-2*n,-n)){t=n;break}const n=s.slice(-t);let e=s;for(;e.endsWith(n);)e=e.slice(0,-t);const r=BigInt(e||"0"),o=BigInt(e+n);u=F(I(BigInt(a||"0")),I(o-r,(10n**BigInt(t)-1n)*10n**BigInt(e.length)))}else u=I(BigInt(a+s||"0"),10n**BigInt(s.length));return o&&(u=o>0?A(u,I(10n**BigInt(o))):O(u,I(10n**BigInt(-o)))),e<0n?S(u):u}function z(t,n=0,e=0,r=null){return 0n===t.n?{kind:"x",c:$,pi:0,e:0,root:null}:{kind:"x",c:t,pi:n,e:e,root:r}}const W=(t,n=1n)=>z(I(BigInt(t),BigInt(n))),q=z(B,1),D=z(B,0,1),_=t=>null!==t&&"object"==typeof t&&"x"===t.kind,C=t=>_(t)&&!t.pi&&!t.e&&!t.root,T=t=>_(t)?0n===t.c.n:0===t,G=t=>C(t)&&1n===t.c.d?t.c.n:null,J=(t,n)=>Object.assign(t,{group:n}),L=t=>({kind:"s",text:t});function R(t,n){return t&&n?t.n===n.n&&t.r===n.r:t===n}const X=(t,n)=>0===P(t.c,n.c)&&t.pi===n.pi&&t.e===n.e&&R(t.root,n.root);function Z(t){if("number"==typeof t)return t;if(_(t)){let n=k(t.c);return t.pi&&(n*=Math.PI**t.pi),t.e&&(n*=Math.E**t.e),t.root&&(n*=k(I(t.root.r))**(1/t.root.n)),n}if("boolean"==typeof t)return t?1:0;throw new p("not a number")}function H(t,...n){if(Number.isNaN(t))throw new d("undefined");if(!Number.isFinite(t)&&n.every(Number.isFinite))throw new p("too large to calculate");return t}const K=t=>"boolean"==typeof t?t:_(t)?0n!==t.c.n:"number"==typeof t?0!==t:Boolean(t);function Q(t){if("bigint"==typeof t)return z(I(t));if("number"==typeof t)return Number.isFinite(t)?z(Number.isInteger(t)&&Math.abs(t)<=Number.MAX_SAFE_INTEGER?I(BigInt(t)):E(String(t))):t;if("boolean"==typeof t)return t;if("string"==typeof t)return/^-?(\d+\.?\d*|\.\d+)(\.\.)?$/.test(t.trim())?z(E(t)):kt(e(t),{});if(Array.isArray(t))return J(t.map(Q),t.group||"()");if(_(t)||"function"==typeof t)return t;throw new p(`can't use ${typeof t} as a value`)}function U(t){return Array.isArray(t)?t.map(U):t&&"s"===t.kind?t.text:"boolean"==typeof t?t:Z(t)}function V(t,n){if(t<2n)return t;const e=BigInt(n);let r=1n<<BigInt(Math.ceil(v(t)/n));for(;;){const n=((e-1n)*r+t/r**(e-1n))/e;if(n>=r)return r;r=n}}function Y(t,n,e,r=0,o=0){if(0n===e)return z($);const{coef:i,root:a}=function(t,n){let e=1n,r=n;const o=BigInt(t),i=r<10n**24n?1000000n:10000n;for(let t=2n;t<=i&&t*t<=r;t+=2n===t?1n:2n){const n=t**o;for(;r%n===0n;)r/=n,e*=t}const a=V(r,t);if(a**o===r&&(e*=a,r=1n),1n===r)return{coef:e,root:null};for(let n=t;n>=2;n--){if(t%n!==0)continue;const o=V(r,n);if(o**BigInt(n)===r)return{coef:e,root:{n:t/n,r:o}}}return{coef:e,root:{n:t,r:r}}}(n,e);return a&&1===a.n?z(A(t,I(i*a.r)),r,o):z(A(t,I(i)),r,o,a)}function tt(t,n){const e=A(t.c,n.c);if(v(e.n)+v(e.d)>M.bits)return null;const r=t.pi+n.pi,o=t.e+n.e;if(!t.root||!n.root)return z(e,r,o,t.root||n.root);const i=Bt(t.root.n,n.root.n);if(i>M.rootIndex)return null;const a=t.root.r**BigInt(i/t.root.n)*n.root.r**BigInt(i/n.root.n);return v(a)>4e3?null:Y(e,i,a,r,o)}function nt(t){if(0n===t.c.n)throw new d("division by zero");const n=z(O(B,t.c),-t.pi,-t.e);if(!t.root)return n;const{n:e,r:r}=t.root;return tt(n,Y(I(1n,r),e,r**BigInt(e-1)))}function et(t,n){if(0n===n){if(0n===t.c.n)throw new d("0^0 is undefined");return z(B)}if(n<0n){const e=nt(t);return e&&et(e,-n)}if((v(t.c.n)+v(t.c.d)+(t.root?v(t.root.r):0))*Number(n)>M.bits||n>100000n)return null;const e=I(t.c.n**n,t.c.d**n),r=Number(n);if(!t.root)return z(e,t.pi*r,t.e*r);const{n:o,r:i}=t.root,a=i**(n/BigInt(o));return Y(A(e,I(a)),o,i**(n%BigInt(o))||1n,t.pi*r,t.e*r)}function rt(t,n,e){if(!C(t)||e>1000n)return null;const r=et(t,n);if(!r)return null;let{n:o,d:i}=r.c,a=!1;if(o<0n){if(e%2n==0n)throw new d("not a real number");a=!0,o=-o}const s=function(t,n){const e=t.n*t.d**BigInt(n-1);return v(e)>4e3?null:Y(I(1n,t.d),n,e)}({n:o,d:i},Number(e));return s?a?z(S(s.c),s.pi,s.e,s.root):s:null}function ot(t){return n=>Array.isArray(n)?J(n.map(n=>ot(t)(n)),n.group):t(n)}function it(t){const n=(e,r)=>{const o=Array.isArray(e),i=Array.isArray(r);if(o&&i){if(e.length!==r.length)throw new d("lists of different lengths");return J(e.map((t,e)=>n(t,r[e])),e.group)}return o?J(e.map(t=>n(t,r)),e.group):i?J(r.map(t=>n(e,t)),r.group):t(e,r)};return n}const at=it((t,n)=>_(t)&&_(n)&&function(t,n){if(0n===t.c.n)return n;if(0n===n.c.n)return t;if(t.pi!==n.pi||t.e!==n.e||!R(t.root,n.root))return null;const e=F(t.c,n.c);return v(e.n)+v(e.d)>M.bits?null:z(e,t.pi,t.e,t.root)}(t,n)||H(Z(t)+Z(n),Z(t),Z(n))),st=ot(t=>_(t)?z(S(t.c),t.pi,t.e,t.root):-Z(t)),ut=(t,n)=>at(t,st(n)),ct=it((t,n)=>_(t)&&_(n)&&tt(t,n)||H(Z(t)*Z(n),Z(t),Z(n))),lt=it((t,n)=>{if(T(n))throw new d("division by zero");if(_(t)&&_(n)){const e=nt(n),r=e&&tt(t,e);if(r)return r}return H(Z(t)/Z(n),Z(t),Z(n))}),ft=it((t,n)=>{if(_(t)&&C(n)){const{n:e,d:r}=n.c,o=1n===r?et(t,e):rt(t,e,r);if(o)return o}const e=Z(t),r=Z(n);if(0===e&&0===r)throw new d("0^0 is undefined");if(0===e&&r<0)throw new d("division by zero");if(e<0&&!Number.isInteger(r))throw new d("not a real number");return H(e**r,e,r)});function ht(t,n){if(C(t)&&C(n))return P(t.c,n.c);if(_(t)&&_(n)&&X(t,n))return 0;const e=Z(t),r=Z(n);return e===r||Number.isFinite(e)&&Number.isFinite(r)&&Math.abs(e-r)<=1e-12*Math.max(1,Math.abs(e),Math.abs(r))?0:e<r?-1:1}function mt(t,n){return Array.isArray(t)||Array.isArray(n)?Array.isArray(t)&&Array.isArray(n)&&t.length===n.length&&t.every((t,e)=>mt(t,n[e])):t&&"s"===t.kind?Boolean(n&&"s"===n.kind&&t.text===n.text):"boolean"==typeof t||"boolean"==typeof n?t===n:0===ht(t,n)}const gt={0:()=>W(0),2:()=>W(1,2),3:()=>Y(I(1n,2n),2,2n),4:()=>Y(I(1n,2n),2,3n),6:()=>W(1)};function bt(t){t=(t%24+24)%24;const[n,e]=t<=6?[t,1]:t<=12?[12-t,1]:t<=18?[t-12,-1]:[24-t,-1],r=gt[n];return r?e<0?st(r()):r():null}const dt=t=>bt(t+6),pt=(t,n)=>Math.abs(t)<1e-15&&Math.abs(n)>1e-6?0:t;function Mt(t,n){const e=function(t){if(T(t))return 0;if(!_(t)||1!==t.pi||t.e||t.root)return null;const n=A(t.c,I(12n));return 1n!==n.d?null:Number((n.n%24n+24n)%24n)}(n);if(null!==e){const n=bt(e),r=dt(e);if(n&&r)switch(t){case"sin":return n;case"cos":return r;case"tan":return lt(n,r);case"cot":return lt(r,n);case"sec":return lt(W(1),r);case"csc":return lt(W(1),n)}}const r=Z(n),o=pt(Math.sin(r),r),i=pt(Math.cos(r),r),a={sin:o,cos:i,tan:o/i,cot:i/o,sec:1/i,csc:1/o}[t];if(!Number.isFinite(a)&&Number.isFinite(r))throw new d("undefined");return H(a)}function wt(t,n){const e={arcsin:[-6,6],arccos:[0,12],arctan:[-5,5]}[t],r=n=>"arcsin"===t?bt(n):"arccos"===t?dt(n):dt(n)&&!T(dt(n))?lt(bt(n),dt(n)):null;if(_(n))for(let t=e[0];t<=e[1];t++){const e=r(t);if(e&&_(e)&&X(e,n))return z(I(BigInt(t),12n),1)}const o=Z(n);if("arctan"!==t&&Math.abs(o)>1)throw new d("not a real number");return{arcsin:Math.asin,arccos:Math.acos,arctan:Math.atan}[t](o)}function yt(t){if(_(t)){if(0n===t.c.n)throw new d("ln 0 is undefined");if(t.c.n<0n)throw new d("not a real number");if(1n===t.c.n&&1n===t.c.d&&!t.pi&&!t.root)return W(t.e)}const n=Z(t);if(n<=0)throw new d(0===n?"ln 0 is undefined":"not a real number");return Math.log(n)}function Nt(t){if(t<.5)return Math.PI/(Math.sin(Math.PI*t)*Nt(1-t));const n=[.9999999999998099,676.5203681218851,-1259.1392167224028,771.3234287776531,-176.6150291621406,12.507343278686905,-.13857109526572012,9984369578019572e-21,1.5056327351493116e-7];let e=n[0];const r=(t-=1)+7+.5;for(let r=1;r<9;r++)e+=n[r]/(t+r);return Math.sqrt(2*Math.PI)*r**(t+.5)*Math.exp(-r)*e}function vt(t,n){const e=G(t)??("number"==typeof t&&Number.isInteger(t)?BigInt(t):null);if(null===e)throw new d(`${n} needs whole numbers`);return e}const xt=ot(t=>C(t)?z(I(j(t.c))):Q(Math.floor(Z(t)))),It=ot(t=>C(t)?z(I(-j(S(t.c)))):Q(Math.ceil(Z(t))));function $t(t,n){const e=Number(vt(n,"round"));if(C(t)){const n=I(10n**BigInt(Math.abs(e))),r=e>=0?A(t.c,n):O(t.c,n),o=j(F(A(I(2n),I(N(r.n),r.d)),B))/2n,i=I(r.n<0n?-o:o);return z(e>=0?O(i,n):A(i,n))}const r=10**e;return Math.sign(Z(t))*Math.round(Math.abs(Z(t))*r)/r}const Bt=(t,n)=>t/Ft(t,n)*n,Ft=(t,n)=>n?Ft(n,t%n):t,At={pi:()=>q,"π":()=>q,tau:()=>z(I(2n),1),"τ":()=>z(I(2n),1),e:()=>D,"∞":()=>1/0,"-∞":()=>-1/0,true:()=>!0,false:()=>!1},Ot={abs:t=>_(t)?z(I(N(t.c.n),t.c.d),t.pi,t.e,t.root):Math.abs(Z(t)),floor:xt,ceil:It,sqrt:t=>ft(t,W(1,2)),ln:ot(yt),loglog:ot(t=>yt(yt(t))),logloglog:ot(t=>yt(yt(yt(t)))),expm1:ot(t=>T(t)?W(0):Math.expm1(Z(t))),sinh:ot(t=>T(t)?W(0):Math.sinh(Z(t))),cosh:ot(t=>T(t)?W(1):Math.cosh(Z(t))),tanh:ot(t=>T(t)?W(0):Math.tanh(Z(t))),arcsinh:ot(t=>T(t)?W(0):Math.asinh(Z(t))),arccosh:ot(t=>{const n=Z(t);if(n<1)throw new d("not a real number");return 1===n?W(0):Math.acosh(n)}),arctanh:ot(t=>{const n=Z(t);if(Math.abs(n)>=1)throw new d("undefined");return T(t)?W(0):Math.atanh(n)}),"!":ot(function(t){const n=G(t)??("number"==typeof t&&Number.isInteger(t)?BigInt(t):null);if(null!==n){if(n<0n)throw new d("factorial of a negative number");if(n>BigInt(M.factorial))throw new p("too large to calculate");let t=1n;for(let e=2n;e<=n;e++)t*=e;return z(I(t))}const e=Z(t);if(e<0&&Number.isInteger(e))throw new d("factorial of a negative number");return H(Nt(e+1),e)}),neg:st,pm:t=>J([t,st(t)],"pm"),"\\not":t=>!K(t)};for(const t of["sin","cos","tan","csc","sec","cot"])Ot[t]=ot(n=>Mt(t,n));for(const t of["arcsin","arccos","arctan"])Ot[t]=ot(n=>wt(t,n));const St={"+":at,"-":ut,"*":ct,"/":lt,"^":ft,"±":(t,n)=>J([at(t,n),ut(t,n)],"pm"),mod:it(function(t,n){if(T(n))throw new d("division by zero");if(C(t)&&C(n))return z(F(t.c,S(A(n.c,I(j(O(t.c,n.c)))))));const e=Z(t),r=Z(n);return H((e%r+r)%r)}),log:it(function(t,n){if(_(t)&&X(t,D))return yt(n);const e=Z(t),r=Z(n);if(e<=0||1===e)throw new d("not a valid base");if(r<=0)throw new d(0===r?"log 0 is undefined":"not a real number");const o=Math.log(r)/Math.log(e);if(C(t)&&C(n))for(let e=1n;e<=12n;e++){const r=Math.round(o*Number(e));if(Math.abs(o*Number(e)-r)>1e-9)continue;const i=et(t,BigInt(r)),a=et(n,e);if(i&&a&&X(i,a))return z(I(BigInt(r),e))}return o}),root:it((t,n)=>ft(n,lt(W(1),t))),atan2:(t,n)=>Math.atan2(Z(t),Z(n)),round:$t,fixed:function(t,n){const e=Number(vt(n,"fixed"));if(e<0)throw new d("fixed needs 0 or more decimal places");const r=$t(t,n);if(C(r)){const t=r.c.n*10n**BigInt(e)/r.c.d,n=(t<0n?-t:t).toString().padStart(e+1,"0");return L(`${t<0n?"-":""}${n.slice(0,n.length-e)}${e?"."+n.slice(-e):""}`)}return L(Number(Z(r)).toFixed(e))},perm:function(t,n){const e=vt(t,"permutations"),r=vt(n,"permutations");if(r<0n||r>e)throw new d("undefined");if(r>BigInt(M.terms))throw new p("too large to calculate");let o=1n;for(let t=e-r+1n;t<=e;t++)o*=t;return z(I(o))},comb:function(t,n){const e=vt(t,"combinations");let r=vt(n,"combinations");if(r<0n||r>e)throw new d("undefined");if(e-r<r&&(r=e-r),r>BigInt(M.terms))throw new p("too large to calculate");let o=1n;for(let t=1n;t<=r;t++)o=o*(e-r+t)/t;return z(I(o))},"=":(t,n)=>mt(t,n),"===":(t,n)=>mt(t,n),"≠":(t,n)=>!mt(t,n),"≠=":(t,n)=>!mt(t,n),"≈":(t,n)=>Math.abs(Z(t)-Z(n))<=1e-9*Math.max(1,Math.abs(Z(t)),Math.abs(Z(n))),"<":(t,n)=>ht(t,n)<0,">":(t,n)=>ht(t,n)>0,"≤":(t,n)=>ht(t,n)<=0,"≥":(t,n)=>ht(t,n)>=0,"&":(t,n)=>K(t)&&K(n),"|":(t,n)=>K(t)||K(n),"\\xor":(t,n)=>K(t)!==K(n),"\\nor":(t,n)=>!(K(t)||K(n)),"\\nand":(t,n)=>!(K(t)&&K(n))},Pt={sum:t=>t.reduce((t,n)=>at(t,n),W(0)),prod:t=>t.reduce((t,n)=>ct(t,n),W(1)),max:t=>t.reduce((t,n)=>ht(n,t)>0?n:t),min:t=>t.reduce((t,n)=>ht(n,t)<0?n:t),mean:t=>lt(Pt.sum(t),W(t.length)),gcd:t=>z(I(t.map(t=>vt(t,"gcd")).reduce(x)))},jt=["∈","∉","∪","∩","pmod","random","(]","[)"];function kt(e,o){if(++w>M.steps)throw new p("too much to work out");switch(e.type){case"number":return z(E(e.value));case"string":return L(e.value.slice(1,-1));case"symbol":return Et(e.value,o)}const i=n(e);if(i)return kt(i,o);const{op:a,args:s}=e,u=t=>kt(t,o);if("?:"===a)return K(u(s[0]))?u(s[1]):u(s[2]);if("_"===a){if(s.every(t=>"op"!==t.type))return Et(`${s[0].value}_${s[1].value}`,o);throw new p("subscript")}if("der"===a)return function(t,n){const[e,r]=t.args,o=Wt(r);if(!o||!Object.prototype.hasOwnProperty.call(n,o))throw new p(`a derivative needs a value for ${o||"its variable"}`);const i=Z(n[o]),a=zt(e,o,n);if(!Number.isFinite(a(i)))throw new d("undefined");const s=1e-6*Math.max(1,Math.abs(i)),u=(a(i+s)-a(i))/s,c=(a(i)-a(i-s))/s;if(!Number.isFinite(u)||!Number.isFinite(c)||Math.abs(u-c)>.001*Math.max(1,Math.abs(u)))throw new d("not differentiable here");let l=.1*Math.max(1,Math.abs(i));const f=1.4,h=[[(a(i+l)-a(i-l))/(2*l)]];let m=h[0][0],g=1/0;for(let t=1;t<12;t++){l/=f,h[t]=[(a(i+l)-a(i-l))/(2*l)];let n=f*f;for(let e=1;e<=t;e++){h[t][e]=(h[t][e-1]*n-h[t-1][e-1])/(n-1),n*=f*f;const r=Math.max(Math.abs(h[t][e]-h[t][e-1]),Math.abs(h[t][e]-h[t-1][e-1]));r<=g&&(g=r,m=h[t][e])}if(Math.abs(h[t][t]-h[t-1][t-1])>=2*g)break}if(!Number.isFinite(m))throw new d("not differentiable here");return Ct(m,g)}(e,o);if("int"===a)return function(t,n){const[e,o,i,a]=t,s=t=>t&&"symbol"===t.type&&/^[a-zA-Zʹ-Ͽ]$/.test(t.value);let u,c,l;if(4===t.length)[u,c,l]=[Wt(o),i,a];else{if(3!==t.length||s(o))throw new p("an integral needs limits to be calculated");[u,c,l]=[r(e)||"x",o,i]}let f=Z(kt(c,n)),h=Z(kt(l,n));if(f===h)return W(0);let m=1;f>h&&([f,h]=[h,f],m=-1);const g=zt(e,u,n),{value:b,error:d}=qt(g,f,h);return Ct(m*b,d)}(s,o);if("lim"===a)return function(t,n){const[e,o,i,a]=t;let s,u,c=0;if(2===t.length&&"op"===o.type&&"→"===o.op)[s,u]=[Wt(o.args[0]),o.args[1]];else if(2===t.length)[s,u]=[r(e),o];else{if(!(t.length>=3))throw new p("a limit needs a point");[s,u]=[Wt(o),i]}if(4===t.length){const t="op"===a.type?"":String(a.value).toLowerCase();c=t.startsWith("r")||"+"===t?1:-1}if(!s)throw new p("a limit needs a variable");const l=kt(u,n),f=Z(l),h=zt(e,s,n);let m,g;if(Number.isFinite(f))if(c)({value:m,error:g}=_t(t=>h(f+c*t*Math.max(1,Math.abs(f)))));else{const t=_t(t=>h(f+t*Math.max(1,Math.abs(f)))),r=_t(t=>h(f-t*Math.max(1,Math.abs(f))));if(t.value===r.value)({value:m,error:g}=t);else{if(!(Number.isFinite(t.value)&&Number.isFinite(r.value)&&Math.abs(t.value-r.value)<=1e-6*Math.max(1,Math.abs(t.value))))throw new d("the limit does not exist");m=(t.value+r.value)/2,g=Math.max(t.error,r.error,Math.abs(t.value-r.value))}try{const t=kt(e,{...n,[s]:l});if(Math.abs(Z(t)-m)<=1e-6*Math.max(1,Math.abs(m)))return t}catch(t){if(!(t instanceof d))throw t}}else({value:m,error:g}=_t(t=>h(Math.sign(f)/t)));return Number.isFinite(m)?Ct(m,g):m}(s,o);if("\\sumidx"===a||"\\prodidx"===a)return function(t,n,e){const[o,i,a]=n;if(n.length<3)throw new p("a sum needs its limits");let s,u=i;"op"===i.type&&"="===i.op&&"symbol"===i.args[0].type?[s,u]=[i.args[0].value,i.args[1]]:s=r(o)||"n";const c=vt(kt(u,e),"a sum"),l=kt(a,e),f=t?at:ct;if(l===1/0){if(!t)throw new p("an infinite product");return function(t,n){const e=[],r=[];let o=0,i=n;for(let a=4;a<=17;a++){const s=2**a;for(;i<n+s;i++){const n=t(i);if(!Number.isFinite(n))throw new d("undefined");o+=n}if(e.push(s),r.push(o),a>=6&&r[r.length-1]===r[r.length-2])return Ct(o,0)}const a=r[r.length-1]-r[r.length-2],s=r[r.length-2]-r[r.length-3];if(Math.abs(a)>=.9*Math.abs(s)&&Math.abs(a)>1e-12)throw new d("the series diverges");const u=e.slice(-7).map(t=>1/t),c=r.slice(-7),l=r.slice(-8,-1),f=(t,n)=>{const e=n.slice();for(let n=1;n<t.length;n++)for(let r=0;r<t.length-n;r++)e[r]=(t[r+n]*e[r]-t[r]*e[r+1])/(t[r+n]-t[r]);return e[0]},h=f(u,c),m=f(e.slice(-8,-1).map(t=>1/t),l);return Ct(h,Math.abs(h-m))}(zt(o,s,e),Number(c))}const h=vt(l,"a sum");if(h-c+1n>BigInt(M.terms))throw new p("too many terms to add one by one");let m=W(t?0:1);for(let t=c;t<=h;t++)m=f(m,kt(o,{...e,[s]:z(I(t))}));return m}("\\sumidx"===a,s,o);if(["()","[]","{}"].includes(a))return J(s.map(u),a);if(jt.includes(a))return function(n,e){const r={};for(const[t,n]of Object.entries(e))"function"!=typeof n&&(r[t]=U(n));try{return Q(function(n,e={}){0===m&&(h=0),m++;try{return g(t(n),e)}finally{m--}}(n,r))}catch(t){throw new p(t.message)}}(e,o);if(a.startsWith("\\:")){const t=o[a.slice(2)];if("function"!=typeof t)throw new p(`no function ${a.slice(2)}`);return Q(t(...s.map(t=>U(u(t)))))}if(Ot[a]&&1===s.length)return Ot[a](u(s[0]));if(St[a]&&2===s.length)return St[a](u(s[0]),u(s[1]));if(Pt[a])return Pt[a](s.map(u));throw new p(`can't calculate '${a}'`)}function Et(t,n){if(Object.prototype.hasOwnProperty.call(n,t))return n[t];if(Object.prototype.hasOwnProperty.call(At,t))return At[t]();throw new p(`no value for '${t}'`)}function zt(t,n,e){return r=>{try{return Z(kt(t,{...e,[n]:r}))}catch(t){if(t instanceof d)return NaN;throw t}}}const Wt=t=>"symbol"===t.type?t.value:null;function qt(t,n,e){if(n===-1/0&&e===1/0){const n=qt(t,-1/0,0),e=qt(t,0,1/0);return{value:n.value+e.value,error:n.error+e.error}}return e===1/0?Dt(e=>t(n+e/(1-e))/((1-e)*(1-e)),0,1):n===-1/0?Dt(n=>t(e-n/(1-n))/((1-n)*(1-n)),0,1):Dt(t,n,e)}function Dt(t,n,e){const r=(e-n)/2;let o=1,i=0,a=NaN,s=NaN,u=0;const c=o=>{const i=Math.PI/2*Math.sinh(Math.abs(o)),a=Math.exp(-2*i),s=2*r*a/(1+a),c=r*(Math.PI/2)*Math.cosh(o)*(4*a)/((1+a)*(1+a));if(0===s||0===c)return 0;const l=o>0?e-s:n+s;if(0!==o&&(l<=n||l>=e))return 0;const f=t(0===o?(n+e)/2:l);if(!Number.isFinite(f)){if(s<1e-6*r)return 0;throw new d("the integral diverges")}return Math.abs(o)>3.5&&(u=Math.max(u,Math.abs(c*f))),c*f};i=c(0);for(let t=o;t<=4;t+=o)i+=c(t)+c(-t);a=o*i;for(let t=1;t<=8;t++){o/=2;for(let t=o;t<=4;t+=2*o)i+=c(t)+c(-t);s=a,a=o*i;const n=Math.max(1,Math.abs(a));if(u>1e-8*n)throw new d("the integral diverges");if(t>=3&&Math.abs(a-s)<=1e-13*n)break}const l=Math.abs(a-s);if(!Number.isFinite(a)||l>1e-6*Math.max(1,Math.abs(a)))throw new d("the integral does not converge");return{value:a,error:l}}function _t(t){const n=[],e=[];for(let r=0;r<30;r++){const o=.1/2**r,i=t(o);Number.isFinite(i)&&(n.push(o),e.push(i))}if(e.length<6)throw new d("undefined");const r=e.slice(1).map((t,n)=>t-e[n]).slice(-6);if(r.every(t=>Math.sign(t)===Math.sign(r[0])&&0!==t)&&Math.abs(r[5])>=.95*Math.abs(r[0])&&Math.abs(e[e.length-1])>10)return{value:Math.sign(r[0])*(1/0),error:0};const o=t=>{const r=n.slice(t,t+8),o=e.slice(t,t+8);for(let t=1;t<r.length;t++)for(let n=0;n<r.length-t;n++)o[n]=(r[n+t]*o[n]-r[n]*o[n+1])/(r[n+t]-r[n]);return o[0]};let i=null;for(let t=0;t+9<=e.length;t++){const n=o(t),e=o(t+1),r=Math.abs(n-e);Number.isFinite(n)&&(!i||r<i.error)&&(i={value:n,error:r})}const a=e[e.length-1],s=Math.abs(a-e[e.length-2]);if(!i||i.error>Math.max(s,1e-6*Math.max(1,Math.abs(a)))){if(s<=1e-6*Math.max(1,Math.abs(a)))return{value:a,error:10*s};throw new d("the limit could not be found")}return i}function Ct(t,n=0){if(!Number.isFinite(t))return t;const e=Math.max(1e-10,4*n)*Math.max(1,Math.abs(t));if(e>1e-6*Math.max(1,Math.abs(t)))return t;if(Math.abs(t)<e)return W(0);const r=[[W(1),1e3],[q,100],[z(B,2),100],[D,100],...[2n,3n,5n,6n,7n].map(t=>[z(B,0,0,{n:2,r:t}),100])];for(const[n,o]of r){const r=t/Z(n),i=Tt(r,o);if(i&&Math.abs(Number(i.n)/Number(i.d)-r)*Z(n)<=e)return tt(z(i),n)}return t}function Tt(t,n){if(!Number.isFinite(t)||Math.abs(t)>1e12)return null;let[e,r,o,i]=[0,1,1,0],a=t;for(let t=0;t<40;t++){const t=Math.floor(a),s=t*r+e,u=t*i+o;if(u>n)break;if([e,r,o,i]=[r,s,i,u],Math.abs(a-t)<1e-15)break;a=1/(a-t)}return i?I(BigInt(r),BigInt(i)):null}function Gt(t,n){if("boolean"==typeof t)return String(t);if(t&&"s"===t.kind)return JSON.stringify(t.text);if(Array.isArray(t)){if("pm"===t.group&&2===t.length&&mt(t[1],st(t[0])))return`±${Jt(t[0],n)}`;const[e,r]={"()":"()","[]":"[]","{}":"{}"}[t.group]||"{}";return e+t.map(t=>Gt(t,n)).join(", ")+r}return"number"==typeof t?Lt(t,n):n.decimal||null!==n.places?function(t,n){if(C(t)&&null===n.places){let e=t.c.d,r=0,o=0;for(;e%2n==0n;)e/=2n,r++;for(;e%5n==0n;)e/=5n,o++;const i=Math.max(r,o);if(1n===e&&i<=n.digits){const e=t.c.n*10n**BigInt(i)/t.c.d,r=e<0n,o=N(e).toString().padStart(i+1,"0"),a=i?`${o.slice(0,-i)}.${o.slice(-i)}`:o;if(o.length<=n.maxDigits)return(r?"-":"")+Rt(a)}}return C(t)&&null!==n.places?Gt($t(t,W(n.places)),{...n,places:null,decimal:!0}):Lt(Z(t),n)}(t,n):function(t,n){const{c:e}=t,r=N(e.n);if(r.toString().length>n.maxDigits||e.d.toString().length>n.maxDigits)return function(t,n){if(!C(t)||1n!==t.c.d)return Lt(Z(t),n);const e=N(t.c.n).toString(),r=BigInt(e.slice(0,n.digits+1)),o=((r+5n)/10n).toString(),i=e.length-1+(o.length>n.digits?1:0),a=Rt(`${o[0]}.${o.slice(1,n.digits)}`);return`${t.c.n<0n?"-":""}${a}*10^(${i})`}(t,n);const o=(t,n)=>1===n?t:`${t}^${n}`,i=[],a=[];t.pi>0&&i.push(o("pi",t.pi));t.pi<0&&a.push(o("pi",-t.pi));t.e>0&&i.push(o("e",t.e));t.e<0&&a.push(o("e",-t.e));t.root&&i.push(2===t.root.n?`sqrt(${t.root.r})`:`root(${t.root.n}, ${t.root.r})`);1n===r&&i.length||i.unshift(r.toString());1n!==e.d&&a.unshift(e.d.toString());const s=t=>t.length>1?`(${t.join("*")})`:t[0],u=a.length?`${s(i)}/${s(a)}`:i.join("*");return(e.n<0n?"-":"")+u}(t,n)}const Jt=(t,n)=>{const e=Gt(t,n);return/^[\d.]+$|^[a-zπ]+$/i.test(e)?e:`(${e})`};function Lt(t,n){if(t===1/0)return"∞";if(t===-1/0)return"-∞";if(Object.is(t,-0)||0===t)return"0";if(null!==n.places&&Math.abs(t)<1e21)return t.toFixed(n.places);const e=Math.abs(t);if(e>=1e15||e<1e-6){const[e,r]=t.toExponential(n.digits-1).split("e");return`${Rt(e)}*10^(${Number(r)})`}return String(Number(t.toPrecision(n.digits)))}const Rt=t=>t.includes(".")?t.replace(/\.?0+$/,""):t;function Xt(t,n={},r={}){const i="string"==typeof t?e(t):o(t),a={};for(const[t,e]of Object.entries(n))a[t]=Q(e);let s;0===y&&(w=0),y++;try{s=kt(i,a)}finally{y--}const u={decimal:!1,places:null,digits:12,maxDigits:50,...r};void 0===u.places&&(u.places=null);const c=t=>Array.isArray(t)?t.every(c):"number"!=typeof t;return{value:U(s),text:Gt(s,u),exact:c(s)}}function Zt(t,n={}){let e=t.trim();const r={},o=/^(~|≈)\s*(\d+)?/.exec(e);o&&(r.decimal=!0,void 0!==o[2]&&(r.places=Number(o[2])),e=e.slice(o[0].length)),/\d\.\d|\.\d/.test(e)&&(r.decimal=!0);try{const{text:t}=Xt(e,n,r);return t.startsWith('"')||/^\d+(\.\d+)?$/.test(t)?t:`(${t})`}catch(t){if(t instanceof d)return'"undefined"';if(t instanceof p)return`(${e})`;throw t}}i(Zt);for(const t of[M,gt,At,Ot,St,Pt])Object.setPrototypeOf(t,null);export{Xt as a,Zt as c};

package/dist/chunks/draw-Cm4exbg7.js

@@ -0,0 +1,2 @@
+/* Inline Math 0.1.0 — MIT licence; MathZetta not included: map mathzetta, mathzetta/speech and mathzetta/calculate with an import map */
+import{GREEK as t,isCustomFunction as e,unaryNeedsBrackets as r,isMinusOneTimes as s,bracketedSides as i,isImplicitProduct as n,customName as l,needsBrackets as a,derivativeOrder as h,bracketOp as o,precedenceOf as c,isVariable as u,firstLeaf as d,firstVariable as p,walk as f,toPostfix as m,operator as g,VARIADIC as b,parseTree as $,fromPostfix as v}from"mathzetta";import{toSpeech as x}from"mathzetta/speech";const y={"(":{fill:"M0.92,0 C0.1,0.2 0.1,0.8 0.92,1 L1,0.975 C0.48,0.75 0.48,0.25 1,0.025 Z"},")":{fill:"M0.08,0 C0.9,0.2 0.9,0.8 0.08,1 L0,0.975 C0.52,0.75 0.52,0.25 0,0.025 Z"},"{":{fill:"M1,0 C0.45,0 0.5,0.05 0.5,0.2 L0.5,0.38 C0.5,0.46 0.35,0.49 0.05,0.5 C0.35,0.51 0.5,0.54 0.5,0.62 L0.5,0.8 C0.5,0.95 0.45,1 1,1 L1,0.985 C0.68,0.985 0.64,0.95 0.64,0.8 L0.64,0.62 C0.64,0.53 0.55,0.505 0.3,0.5 C0.55,0.495 0.64,0.47 0.64,0.38 L0.64,0.2 C0.64,0.05 0.68,0.015 1,0.015 Z"},"}":{fill:"M0,0 C0.55,0 0.5,0.05 0.5,0.2 L0.5,0.38 C0.5,0.46 0.65,0.49 0.95,0.5 C0.65,0.51 0.5,0.54 0.5,0.62 L0.5,0.8 C0.5,0.95 0.55,1 0,1 L0,0.985 C0.32,0.985 0.36,0.95 0.36,0.8 L0.36,0.62 C0.36,0.53 0.45,0.505 0.7,0.5 C0.45,0.495 0.36,0.47 0.36,0.38 L0.36,0.2 C0.36,0.05 0.32,0.015 0,0.015 Z"},sqrt:{stroke:"M0,0.62 L0.18,0.55 L0.45,0.97 L0.98,0.015"},int:{stroke:"M0.95,0.09 C0.9,0.01 0.72,0 0.66,0.12 C0.6,0.25 0.42,0.75 0.34,0.88 C0.28,0.99 0.1,0.99 0.05,0.91"},sum:{fill:"M0.04,0 L0.94,0 L0.98,0.24 L0.94,0.24 C0.9,0.1 0.8,0.07 0.62,0.07 L0.24,0.07 L0.58,0.5 L0.2,0.93 L0.64,0.93 C0.82,0.93 0.92,0.9 0.95,0.76 L0.99,0.76 L0.94,1 L0.02,1 L0.47,0.47 Z"},prod:{fill:"M0,0 L1,0 L1,0.04 L0.88,0.05 L0.88,0.95 L1,0.96 L1,1 L0.62,1 L0.62,0.96 L0.74,0.95 L0.74,0.06 L0.26,0.06 L0.26,0.95 L0.38,0.96 L0.38,1 L0,1 L0,0.96 L0.12,0.95 L0.12,0.05 L0,0.04 Z"}};function w(t,e,r=""){const s=y[t];return`<svg class="${e}" viewBox="0 0 1 1" preserveAspectRatio="none" aria-hidden="true"${r?` style="${r}"`:""}>${s.fill?`<path d="${s.fill}" fill="currentColor"/>`:`<path d="${s.stroke}" fill="none" stroke="currentColor" vector-effect="non-scaling-stroke"/>`}</svg>`}const L=.8,M=.2,C=.25,k=.7,z="−",P={"+":"add","-":"sub","*":"mul","/":"frac","^":"pow",_:"subscript",neg:"neg",pm:"pm","!":"factorial","=":"eq","===":"eq","≠":"neq","≠=":"neq","<":"lt",">":"gt","≤":"le","≥":"ge","≈":"approx","≡":"equiv","∈":"in","∉":"notin","∪":"union","∩":"intersect","→":"to","±":"pm","&":"and","|":"or","\\not":"not","\\xor":"xor","\\nand":"nand","\\nor":"nor","?:":"cases","()":"tuple","[]":"list","{}":"set","(]":"interval","[)":"interval","\\sumidx":"sum","\\prodidx":"prod",mod:"mod",pmod:"pmod"},j={"=":"=","===":"=","≠":"≠","≠=":"≠","<":"<",">":">","≤":"≤","≥":"≥","∈":"∈","∉":"∉","→":"→","≈":"≈","≡":"≡","&":"∧","|":"∨","\\xor":"⊕","\\nand":"↑","\\nor":"↓"},A={"±":"±","∪":"∪","∩":"∩",mod:"mod"},Z=["sin","cos","tan","csc","sec","cot","sinh","cosh","tanh"],q={arcsin:"sin",arccos:"cos",arctan:"tan",arcsinh:"sinh",arccosh:"cosh",arctanh:"tanh"},B=[...Z,...Object.keys(q),"ln"],F={atan2:"atan2",max:"max",min:"min",gcd:"gcd",mean:"mean",round:"round",fixed:"fixed",expm1:"expm1",random:"random"},O={"()":["(",")"],"[]":["[","]"],"{}":["{","}"],"(]":["(","]"],"[)":["[",")"]},W={"\\partial":"∂","\\nabla":"∇","\\ell":"ℓ","\\infty":"∞","\\emptyset":"∅"},S={vec:"→",hat:"ˆ",widehat:"ˆ",tilde:"˜",dot:"˙",ddot:"¨"},E=["boxed","overline","bar","underline","cancel",...Object.keys(S)],_={prefix:"mv",display:!0,fnParens:"auto",trigPowers:!0,inverseTrig:"arc",derivativeBrackets:"[",scriptFractions:"slash",cases:!0,label:!0,attributes:null,style:null},I=t=>String(t).replace(/&/g,"&amp;").replace(/</g,"&lt;").replace(/>/g,"&gt;").replace(/"/g,"&quot;"),N=t=>+t.toFixed(3)+"em",T=(t,e=L,r=M)=>({html:t,h:e,d:r});function D(...t){const e=t.flat().filter(Boolean);return{html:e.map(t=>t.html).join(""),h:Math.max(L,...e.map(t=>t.h)),d:Math.max(M,...e.map(t=>t.d))}}function R(t,e={}){const r=new G({..._,...e}),s=r.render(t,{level:0}),i=`${r.p}${r.o.display?` ${r.p}-display`:""}`,n=function(t,e){if(!1===e||null===e)return null;if("string"==typeof e)return e;try{return"function"==typeof e?e(t):x(t)}catch{return null}}(t,r.o.label),l=null===n?s.html:s.html.replace(/^<span /,'<span aria-hidden="true" ');return{html:`<span class="${i}"${null===n?"":` role="math" aria-label="${I(n)}"`}>${l}</span>`,h:s.h,d:s.d}}class G{constructor(t){this.o=t,this.p=t.prefix}el(t,e,r="",s=""){return`<span class="${t.split(" ").map(t=>`${this.p}-${t}`).join(" ")}"${r?` style="${r}"`:""}${s}>${e}</span>`}text(t,e){const r=T(this.el(t,I(e)));return r.leaf={kind:t,content:I(e)},r}bin(t){return this.text("bin",t)}rel(t){return this.text("rel",t)}punct(t){return this.text("punct",t)}word(t){return this.text("word",t)}wrap(t,e,r){let s="";void 0!==t.id&&(s+=` data-node="${t.id}"`),"op"===t.type&&(s+=` data-op="${I(t.op)}"`);const i=this.o.attributes&&this.o.attributes(t)||{};for(const[t,e]of Object.entries(i))"class"!==t&&"style"!==t&&(s+=` ${t}="${I(e)}"`);const n=[i.style,this.o.style&&this.o.style(t)].filter(Boolean).join(";"),l=(r.leaf?["n",...e,r.leaf.kind]:["n",...e]).map(t=>`${this.p}-${t}`).concat(i.class?[i.class]:[]).join(" "),a=r.leaf?r.leaf.content:r.html;return{html:`<span class="${l}"${n?` style="${I(n)}"`:""}${s}>${a}</span>`,h:r.h,d:r.d}}render(t,e){return"number"===t.type?this.number(t):"string"===t.type?this.wrap(t,["text"],this.text("text",t.value.replace(/^["'`]|["'`]$/g,""))):"symbol"===t.type?this.wrap(t,["symbol"],this.symbol(t.value)):this.operation(t,e)}number(t){return this.wrap(t,["number"],this.text("num",t.value.replace(/^-/,z)))}symbol(e){return t[e]&&(e=t[e]),W[e]?this.text("ord",W[e]):"\\+"===e||"\\-"===e?this.text("ord",e[1]):"-∞"===e?this.text("ord","−∞"):/^[A-Za-z]$/.test(e)||/^[α-ω]$/.test(e)?this.text("var",e):/^[A-Za-z]′+$/.test(e)?T(this.el("var",I(e[0]))+this.el("prime",I(e.slice(1)))):/^Δ[A-Za-z]$/.test(e)?T(this.el("ord","Δ")+this.el("var",e[1])):e.endsWith("⁻¹")?this.script(this.symbol(e.slice(0,-2)),this.text("num","−1"),{level:0}):/^\\[A-Za-z]+$/.test(e)?this.text("name",e.slice(1)):/^[A-Za-z]{2,}$/.test(e)?this.text("name",e):this.text("ord",e)}operation(t,i){const{op:n,args:l}=t,[h,o]=l,c=(t,e={})=>this.render(t,{...i,...e,forceParens:e.forceParens||!1,before:e.before||!1}),u=P[n]??(/^[a-z0-9]+$/i.test(n)?n:"other"),d=e=>this.wrap(t,["op",`op-${u}`],e),p=(e,r)=>a(t,e,r)&&"/"!==e.op?this.fence(c(e)):c(e);if(B.includes(n)||"log"===n)return d(this.named(t,i));if(F[n])return d(D(this.word(F[n]),this.fence(this.list(l,i))));if(e(n))return d(this.custom(t,i));if(O[n])return d(this.fence(this.list(l,i),...O[n]));if(j[n]&&2===l.length)return d(D(p(h,"left"),this.rel(j[n]),p(o,"right")));if(A[n]&&2===l.length)return d(D(p(h,"left"),this.bin(A[n]),p(o,"right")));switch(n){case"+":return s(o)?d(D(p(h,"left"),this.bin(z),this.wrap(o,["op","op-mul"],this.render(o.args[1],i)))):"number"===o.type&&o.value.startsWith("-")?d(D(p(h,"left"),this.bin(z),this.wrap(o,["number"],this.text("num",o.value.slice(1))))):d(D(p(h,"left"),this.bin("+"),this.signed(o,p(o,"right"))));case"-":return d(D(p(h,"left"),this.bin(z),this.signed(o,p(o,"right"))));case"*":return d(this.product(t,i));case"/":return d(this.fraction(t,i));case"^":return d(this.power(t,i));case"_":return d(this.script(c(h,{forceParens:!0}),null,i,c(o,{level:i.level+1})));case"neg":return d(D(this.text("sign",z),r(n,h)?this.fence(c(h)):c(h,{before:i.before})));case"pm":return d(D(this.text("sign","±"),r(n,h)?this.fence(c(h)):c(h)));case"!":return d(D(r(n,h)?this.fence(c(h)):c(h,{forceParens:!0}),this.text("ord","!")));case"\\not":return d(D(this.text("sign","¬"),r(n,h)?this.fence(c(h)):c(h)));case"pmod":return d(D(c(h),this.fence(D(this.word("mod"),this.text("space",""),c(o)))));case"sqrt":return d(this.radical(c(h)));case"root":return d(this.radical(c(o),this.render(h,{level:i.level+2})));case"abs":return d(this.fence(c(h),"|","|"));case"floor":return d(this.fence(c(h),"⌊","⌋"));case"ceil":return d(this.fence(c(h),"⌈","⌉"));case"loglog":return d(D(this.word("ln"),this.fence(D(this.word("ln"),this.fence(c(h))))));case"logloglog":return d(D(this.word("ln"),this.fence(D(this.word("ln"),this.fence(D(this.word("ln"),this.fence(c(h))))))));case"perm":return d(this.script(this.script(T(""),this.render(h,{level:i.level+1}),i),null,i,this.render(o,{level:i.level+1}),this.word("P")));case"comb":return d(this.fence(this.stack(c(h),c(o),!1)));case"?:":return d(this.o.cases?this.cases(t,i):D(c(h),this.bin("?"),c(o),this.bin(":"),c(l[2])));case"sum":return d(D(l.map((t,e)=>e?[this.bin("+"),c(t)]:[c(t)]).flat()));case"prod":return d(D(l.map((t,e)=>e?[this.bin("·"),c(t)]:[c(t)]).flat()));case"der":return d(this.derivative(t,i));case"int":return d(this.integral(l,i));case"lim":return d(this.limit(l,i));case"\\sumidx":return d(this.bigOperator("sum",l,i));case"\\prodidx":return d(this.bigOperator("prod",l,i))}throw new Error(`Inline Math cannot draw '${n}'`)}signed(t,e){return("number"===t.type?t.value.startsWith("-"):"neg"===t.op||"pm"===t.op)&&!e.fenced?this.fence(e):e}list(t,e){return D(t.map((t,r)=>r?[this.punct(","),this.render(t,e)]:[this.render(t,e)]).flat())}product(t,r){const[l,a]=t.args;if(s(t))return D(this.text("sign",z),this.render(a,r));const h=i(t),o=n(l,a,h),c=function(t){for(;"op"===t.type;){if(B.includes(t.op)||"log"===t.op||F[t.op]||e(t.op))return!0;if("^"!==t.op&&"*"!==t.op)return!1;t=t.args[0]}return!1}(a),u=h.left?this.fence(this.render(l,r)):this.render(l,{...r,before:o&&!c}),d=h.right?this.fence(this.render(a,r)):this.render(a,{...r,before:r.before});return D(u,o?T(this.el(c?"fnspace":"times","")):this.bin("·"),d)}named(t,e,r=null){const{op:s,args:i}=t,n="log"===s?i[1]:i[0];let l;l=q[s]&&"power"===this.o.inverseTrig?this.script(this.word(q[s]),this.text("num","−1"),e):this.word(s),r&&(l=this.script(l,r,e)),"log"===s&&(l=this.script(l,null,e,this.render(i[0],{level:e.level+1})));const a="auto"===this.o.fnParens&&!e.forceParens&&!e.before&&H(n),h=this.render(n,{...e,forceParens:!1,before:!1});return D(l,a?D(T(this.el("fnspace","")),h):this.fence(h))}custom(t,e){const r=l(t.op),s=t.args;if(E.includes(r)&&1===s.length)return this.decorate(r,this.render(s[0],e));return D(/^[A-Za-z]′*(⁻¹)?$/.test(r)?this.symbol(r):this.word(r),T(this.el("callgap","")),this.fence(this.list(s,e)))}decorate(t,e){switch(t){case"boxed":return T(this.el("boxed",e.html),e.h+.22,e.d+.22);case"cancel":return T(this.el("cancel",e.html),e.h,e.d);case"underline":return T(this.el("underline",e.html),e.h,e.d+.12);case"overline":case"bar":return T(this.el("over",e.html),e.h+.14,e.d);default:return T(this.el("accent",this.el("accent-mark",S[t])+e.html),e.h+.3,e.d)}}fraction(t,e){const[r,s]=t.args;if(e.level>0&&"slash"===this.o.scriptFractions&&H(r,!0)&&H(s,!0)){const t=t=>"op"!==t.type||H(t)?this.render(t,e):this.fence(this.render(t,e));return D(t(r),this.text("ord","/"),t(s))}const i=e.level>0||e.inFraction||!this.o.display?.8:1,n={...e,inFraction:!0,forceParens:!1,before:!1};return this.stack(this.render(r,n),this.render(s,n),!0,i)}stack(t,e,r=!0,s=1){const i=t=>1===s?t.html:this.el("scaled",t.html,`font-size:${N(s)}`),n=r?.06:0,l=C+n/2+.1+Math.max(M,s*t.d),a=this.el("frac"+(r?"":" frac-nobar"),this.el("numer",i(t))+this.el("denom",i(e)),`vertical-align:${N(l)}`);return T(a,l+Math.max(L,s*t.h),Math.max(L,s*e.h)+Math.max(M,s*e.d)+.1+n/2-C)}power(t,e){const[r,s]=t.args,i={...e,level:e.level+1,forceParens:!1,before:!1,inFraction:!1},n=()=>this.render(s,i);if(this.o.trigPowers&&Z.includes(r.op)&&K(s))return this.wrap(r,["op",`op-${r.op}`],this.named(r,e,n()));if("_"===r.op){const t=this.render(r.args[1],i);return this.wrap(r,["op","op-subscript"],this.script(this.render(r.args[0],{...e,forceParens:!0}),n(),e,t))}const l=a(t,r,"left")||"number"===r.type&&r.value.startsWith("-")||"/"===r.op||B.includes(r.op)?this.fence(this.render(r,e)):this.render(r,{...e,forceParens:!0});return this.script(l,n(),e)}script(t,e,r,s=null,i=null){const n=r.level>0?.8:k,l=i?D(t,i):t;let a=e?Math.max(.42,t.h-.38,n*e.d+.12):0,h=s?Math.max(.2,t.d-.05,n*s.h-.45):0;if(e&&s){const t=.1-(a-n*e.d-(n*s.h-h));t>0&&(h+=t)}const o=t=>this.el("script",t.html,`font-size:${N(n)}`);if(e&&!s)return T(t.html+this.el("sup",e.html,`font-size:${N(n)};vertical-align:${N(a/n)}`),Math.max(l.h,a+n*e.h),Math.max(l.d,n*e.d-a));if(s&&!e)return T(l.html+this.el("sub",s.html,`font-size:${N(n)};vertical-align:${N(-h/n)}`),Math.max(l.h,n*s.h-h),Math.max(l.d,h+n*s.d));const c=a+n*e.h,u=h+n*s.d,d=this.el("subsup",o(e)+o(s),`vertical-align:${N(a)};height:${N(c+u)}`);return T(l.html+d,Math.max(l.h,c),Math.max(l.d,u))}radical(t,e=null){const r=Math.max(t.h,L)+.1+.06,s=Math.max(t.d,M),i=r+s,n=w("sqrt",`${this.p}-hook`,`height:${N(i)};width:${N(.5+.06*(i-1))};vertical-align:${N(-s)}`);let l="";if(e){const t=.8*k;l=this.el("rootindex",e.html,`font-size:${N(t)};vertical-align:${N((.55*i-s)/t)}`)}return T(this.el("sqrt",l+n+this.el("radicand",t.html)),r,s)}fence(t,e="(",r=")"){const s=Math.max(t.h,L)+Math.max(t.d,M),i=s>1.25,n=i?.1:0,l=(t,e)=>{if(!t)return{html:"",w:0};const r={"(":.26,")":.26,"{":.36,"}":.36}[t];if(r){const n=r+(i?Math.min(.2,.05*(s-1)):0);return{html:w(t,`${this.p}-delim ${this.p}-delim-${e}`,`width:${N(n)}`),w:n+.1}}const n={"[":"lbrack","]":"rbrack","|":"vert","⌊":"lfloor","⌋":"rfloor","⌈":"lceil","⌉":"rceil"}[t];return{html:this.el(`delim ${n} delim-${e}`,""),w:"|"===t?.18:.3}},a=l(e,"l"),h=l(r,"r"),o=`padding:${N(n)} ${N(h.w)} ${N(n)} ${N(a.w)}`,c=T(this.el("fence",a.html+t.html+h.html,o),Math.max(t.h,L)+n,Math.max(t.d,M)+n);return c.fenced=!0,c}derivative(t,e){const{n:r,body:s,by:i,variable:n}=h(t),l={...e,level:e.level+1},a=()=>this.text("d","d"),u=()=>r>1?this.render({type:"number",value:String(r)},l):null,d={...e,inFraction:!0},p=D(a(),u()?this.script(this.render(i,d),u(),d):this.render(i,d)),f=t=>D(u()?this.script(a(),u(),d):a(),t||[]),m=e.level>0||e.inFraction||!this.o.display?.8:1;if(n)return this.stack(f(this.render(s,d)),p,!0,m);const g=o(s)&&c(o(s))<100,[b,$]="("===this.o.derivativeBrackets?["(",")"]:["[","]"],v=g?this.fence(this.render(s,e),b,$):this.render(s,e);return D(this.stack(f(),p,!0,m),v)}integral(t,e){const[r,s,i,n]=t,l=t=>this.render(t,e),a=t=>t?this.render(t,{...e,level:e.level+1}):null,h=()=>this.text("var",p(r)||"x");let o,c=null,d=null;switch(t.length){case 1:o=h();break;case 2:u(s)?o=l(s):(o=h(),c=a(s));break;case 3:u(s)?(o=l(s),c=a(i)):(o=h(),c=a(s),d=a(i));break;default:o=l(s),c=a(i),d=a(n)}const f=l(r),m=T(this.el("dx",this.el("d","d")+o.html),o.h,o.d);return D(this.bigSign("int",f,c,d,e),f,m)}bigSign(t,e,r,s,i){const n=!this.o.display||i.level>0,l=Math.max(n?.6:.95,Math.max(e.h-C,e.d+C)+.12),a=w(t,`${this.p}-bigsign ${this.p}-${t}sign`,`height:${N(2*l)};width:${N(n?.38:.5)};vertical-align:${N(C-l)}`);if(!r&&!s)return T(a,C+l,l-C);const h=k,o=s||T(""),c=r||T(""),u=C+l-h*o.h;c.d;const d=this.el("limits",this.el("script",o.html,`font-size:${N(h)}`)+this.el("script lower",c.html,`font-size:${N(h)}`),`vertical-align:${N(u)};height:${N(2*l)}`);return T(a+d,C+l,l-C)}limit(t,e){const[r,s,i,n]=t,l={...e,level:e.level+1},a=t=>this.render(t,l),h=()=>this.rel("→");let o;switch(t.length){case 1:return D(this.word("lim"),T(this.el("fnspace","")),this.render(r,e));case 2:{const t=p(r),e="number"===s.type||"symbol"===s.type&&s.value.includes("∞");o=t&&e?D(this.text("var",t),h(),a(s)):a(s);break}case 3:o=D(a(s),h(),a(i));break;default:{const t=d(n).toLowerCase().startsWith("r")?"+":z;o=D(a(s),h(),this.script(a(i),this.wrap(n,["symbol"],this.text("ord",t)),l))}}return D(this.under(this.word("lim"),o,e),this.render(r,e))}under(t,e,r){return!this.o.display||r.level>0?this.script(t,null,r,e):T(this.el("under",t.html+this.el("script",e.html,`font-size:${N(.7)}`)),t.h,Math.max(t.d,M)+.05+.7*(e.h+e.d))}bigOperator(t,e,r){const[s,i,n]=e,l={...r,level:r.level+1},a=i?this.render(i,l):null,h=n?this.render(n,l):null,o=this.render(s,r);if(!this.o.display||r.level>0){const e=.55,s=w(t,`${this.p}-bigsign`,`height:${N(2*e)};width:${N(.8)};vertical-align:${N(C-e)}`);return D(this.script(T(s,C+e,e-C),h,r,a),o)}const c=.68,u=k,d=w(t,`${this.p}-bigsign`,`height:${N(1.36)};width:${N("sum"===t?1:1.05)}`),p=h?1.01+u*h.d:C-c,f=t=>this.el("script",t.html,`font-size:${N(u)}`),m=this.el("bigop",(h?f(h):"")+d+(a?f(a):""),`vertical-align:${N(p)}`);return D(T(m,.93+(h?.08+u*(h.h+h.d):0),c-C+(a?.08+u*(a.h+a.d):0)),o)}cases(t,e){const r=[];let s=t;for(;"op"===s.type&&"?:"===s.op;)r.push([this.render(s.args[1],e),D(this.word("if"),this.text("space",""),this.render(s.args[0],e))]),s=s.args[2];r.push([this.render(s,e),this.word("otherwise")]);const i=r.map(([t,e])=>({h:Math.max(t.h,e.h),d:Math.max(t.d,e.d)})),n=i.reduce((t,e)=>t+e.h+e.d,0)+.25*(r.length-1),l=C+n/2-i[0].h,a=r.map(([t,e])=>this.el("case-value",t.html)+this.el("case-cond",e.html)).join(""),h=T(this.el("cases",a,`vertical-align:${N(l)}`),C+n/2,n/2-C);return this.fence(h,"{","")}}function H(t,e=!1){return"number"===t.type?!t.value.startsWith("-")||e:"symbol"===t.type||"op"===t.type&&("abs"===t.op?!e:!e||"neg"!==t.op&&"^"!==t.op?"*"===t.op&&!s(t)&&t.args.every(t=>"number"===t.type&&!t.value.startsWith("-")||J(t)||"*"===t.op&&H(t,e)):H(t.args[0],!0)&&("neg"===t.op||H(t.args[1],!0)))}const J=e=>"symbol"===e.type&&(1===e.value.length||/^[α-ωΑ-Ω]$/.test(e.value)||Boolean(t[e.value])),K=t=>"number"===t.type&&/^[1-9]\d*$/.test(t.value)&&Number(t.value)<10;const Q=t=>m(t).join(" ");function U(t,e){const r=new Set;V(t,e,r);const s=[];return f(t,t=>{r.has(t)&&s.push(t)}),s}function V(t,e,r){if(null!=e)if(Array.isArray(e))e.forEach(e=>V(t,e,r));else if("number"!=typeof e)if("function"!=typeof e){if("string"==typeof e){const s=Q($(e));return void f(t,t=>{Q(t)===s&&r.add(t)})}if(e.path){const s=X(t,e.path);return void(s&&r.add(s))}if(e.type){const s=Q(e);f(t,t=>{t!==e&&Q(t)!==s||r.add(t)})}}else f(t,(t,s)=>{e(t,s)&&r.add(t)});else f(t,t=>{t.id===e&&r.add(t)})}function X(t,e){let r=t;for(const t of e){if("op"!==r.type||!r.args[t])return null;r=r.args[t]}return r}function Y(t,e){const r=[],s=t=>{if(t.id===e)return!0;if("op"!==t.type)return!1;for(let e=0;e<t.args.length;e++){if(r.push(e),s(t.args[e]))return!0;r.pop()}return!1};return s(t)?r:null}function tt(t){let e=0;const r=t=>{"op"===t.type&&(t.args.forEach(r),g(t.op)?.arity===b&&e++),t.id=e++};return r(t),t}function et(t,e={}){const{blocks:r}=e;if("string"==typeof t)return $(t,r?{blocks:t=>r(t,e.scope)}:{});if(Array.isArray(t))return v(t);if(t&&"object"==typeof t&&t.type)return void 0===t.id?tt(structuredClone(t)):t;throw new Error("Inline Math needs an expression, standard postfix or a tree")}function rt(t,e={}){const r=et(t,e);return R(r,st(r,e)).html}function st(t,e){const{colors:r,highlight:s,scope:i,blocks:n,...l}=e,a=new Map;for(const[e,s]of function(t){return t?Array.isArray(t)?t:Object.entries(t).map(([t,e])=>[/^\d+$/.test(t)?Number(t):t,e]):[]}(r))if(it(s))for(const r of U(t,e))a.set(r,`color:${s}`);const h=new Set(void 0===s?[]:U(t,s));if(!a.size&&!h.size)return l;const o=l.prefix||"mv";return{...l,style:t=>[a.get(t),l.style&&l.style(t)].filter(Boolean).join(";"),attributes:t=>{const e={...l.attributes?l.attributes(t):{}};return h.has(t)&&(e.class=[e.class,`${o}-hl`].filter(Boolean).join(" ")),e}}}function it(t){return"string"==typeof t&&/^[#\w\s.,%()/+-]+$/.test(t)&&!/url\s*\(|expression/i.test(t)}export{_ as D,y as P,et as a,st as b,Q as c,X as d,R as l,tt as n,Y as p,U as s,rt as t};

package/dist/chunks/html-CvuSyhzk.js

@@ -0,0 +1,2 @@
+/* Inline Math 0.1.0 — MIT licence; MathZetta not included: map mathzetta, mathzetta/speech and mathzetta/calculate with an import map */
+import{calculateBlock as t}from"mathzetta/calculate";import{t as a,a as o}from"./draw-Cm4exbg7.js";function r(a,r={}){return o(a,{...r,blocks:r.blocks??t})}function s(o,r={}){return a(o,{...r,blocks:r.blocks??t})}export{s as a,r as t};

package/dist/chunks/html-DWtfytjA.js

@@ -0,0 +1,2 @@
+/* Inline Math 0.1.0 — MIT licence; includes mathzetta (and @mathzetta/cas, for the CAS) */
+import{c as o}from"./calculate-BX9lJ7Ff.js";import{d as s,e as t}from"./lite-Bj_L6w2g.js";function c(s,c={}){return t(s,{...c,blocks:c.blocks??o})}function r(t,c={}){return s(t,{...c,blocks:c.blocks??o})}export{r as a,c as t};