<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN">
<html lang="en">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<link rel="shortcut icon" href="/favicon.ico">
<link rel="stylesheet" type="text/css" href="https://pari.math.u-bordeaux.fr/pari.css">
<title>PARI/GP Development Headquarters</title>
</head>
<body>
<h1 class="head">Run PARI/GP in your browser
<br> <a href="gp.html">Stable Version</a> | <a href="gpexp.html">Development version</a> | WebAssembly version
</h1>
<div id="gpshell">
<pre class="gpoutput" id="print"></pre>
<pre class="gpoutput" id="output"></pre>
<textarea class="code" id="input" onkeydown="return shiftenter(this,event)"></textarea>
<button type="button" onclick="return eval(this)">Evaluate with PARI</button>
<button type="button" onclick="return clearcell(this)">Clear messages</button>
</div>
<button type="button" onclick="return newcell(this)">New Cell</button>
<br>
Type a GP expression in the third (blue) box then hit
<code>&lt;Shift&gt;-&lt;Enter&gt;</code> or click on the
Evaluate button below it. A trailing semi-colon <code>;</code> prevents output.
You can assign to variables and define user functions. (Speed is approximately one
eighth of native speed.)
<h2>Examples</h2>
<ul>
<li><code>?eulerphi</code> (short help), <code>??factor</code> (long help in
new tab)
<li><code>2*3</code>, <code>10!</code>, <code>(1+x)^10</code>,
<code>exp(Pi*sqrt(163))</code>
  </li>
<li><code>factor(2^128+1)</code>, <code>factor(x^4+4)</code>, <code>isprime(10^50+151)</code>
  </li>
<li><code>
ploth(x=1,100,[sqrt(2/(Pi*x)),besselj(0,x),-sqrt(2/(Pi*x))])</code>
  </li>
<li><code>
ploth(x=0,2*Pi,[sin(2*x),cos(3*x)],1)</code>
  </li>
<li><code>intnum(x=-oo,oo,1/(1+x^2))</code>, <code>sumpos(n=1, n^(-3))</code>
  </li>
<li><code>sin(x)</code>, <code>ellj(x + O(x^20))</code>
  </li>
<li><code>Mod(x,x^2-x-1)^100000==fibonacci(100000)*x+fibonacci(99999)</code>
  </li>
<li><code>g=gamma(1/4+O(5^20))</code>, <code> algdep(g,4)</code>, <code>g^4+4*g^2+5</code>
  </li>
<li><code>\p100</code>: change default accuracy to 100 decimal digits
  </li>
<li><code>K=bnfinit(a^2+23);</code> (huge: don't forget the <code>;</code>),
  <code>K.clgp</code>, <code>bnrclassfield(K)</code>
  </li>
<li><code>E=ellinit([1,3]*Mod(1,10^50+151)); ellcard(E)</code>
  </li>
<li><code>E=ellinit([-157^2,0]); P=ellheegner(E)</code>, <code>ellisoncurve(E,P)</code>
  </li>
<li><code>polgalois(x^8-2)</code>, <code>ellsearch(11)</code> (loads packages
on demand)</li>
<li><code>lfun(x^3-2, 2)</code> (Dedekind zeta function)
<li><code>E=ellinit("11a1");L=lfuninit(E,[0,100]);
ploth(x=0,100,lfunhardy(L,x))</code>
<li>
<pre>
ecm(N, B = 1000!, nb = 100)=
{
  for(a = 1, nb,
    iferr(ellmul(ellinit([a,1], N), [0,1], B),
      E, return(gcd(lift(component(E,2)),N))))
}
ecm(2^101-1)
  </pre>
  </li>
</ul>
<br>
<p>
PARI is compiled using <a href="http://emscripten.org">emscripten</a>.
Enter <code>\v</code> for version information.
</p>
<script type="text/javascript">
var currPrint = document.getElementById("print");
var currOutput = document.getElementById("output");
function cvtohtml(s) {
  s = s.replace(/&/g, "&amp;");
  s = s.replace(/</g, "&lt;");
  s=s.replace( />/g, "&gt;");
  s = s.replace(/\n/g, "<br>");
  return s;
}
var rawPrint = false;
function myprint(text)
{
  if (rawPrint)
  {
    currOutput.innerHTML += (text)+"\n";
  }
  else
  {
    currOutput.innerHTML += cvtohtml(text)+"\n";
  }
}
function myprinterr(text)
{
  currPrint.innerHTML += cvtohtml(text)+"\n";
}
var Module =
{
  noInitialRun: true,
  print: myprint,
  printErr: myprinterr,
  wasmBinaryFile: './gp-sta.wasm'
};
function myevalm(s) { return Module.ccall("gp_embedded","string",["string"],[s]); }
function myinit() { Module.ccall("gp_embedded_init",null,["number","number"],[100000000,100000000]); }
function myplotinit() { Module.ccall("pari_emscripten_plot_init",null,["number","number"],[640,480]); }
var initonce = true;
function searchid(node, id)
{
  var children = node.childNodes;
  for (i in children)
  {
    if (children[i].id == id) { return children[i]; }
  }
}
function evalgp(where, x) {
  var s = myevalm(x);
  if (s.slice(0,6) == "Help: ")
  {
    var url = s.slice(6,-1);
    window.open(url);
  }
  else if (s == "retry.")
  {
    setTimeout(function (){ evalgp(where, x); }, 500);
  }
  else if (s.length> 1)
  {
    var out = searchid(where.parentNode,"output");
    var t = cvtohtml(s);
    var c = s.charAt(0);
    if (c != "%" && c != "\n")
    {
      t = "<em style=\"color: red\">*** " + t + "</em>\n";
    }
    out.innerHTML += t;
  }
}
function eval(where) {
  if (initonce) { myinit(); myplotinit(); initonce = false; }
  currPrint = searchid(where.parentNode,"print");
  var x = searchid(where.parentNode,"input").value;
  currOutput = searchid(where.parentNode,"output");
  currOutput.innerHTML += "? " + cvtohtml(x) + "\n";
  evalgp(where, x);
}
function clearcell(where) {
  searchid(where.parentNode,"print").innerHTML = "";
}
function shiftenter(where,e) {
   if (e.shiftKey && e.keyCode == 13)
   {
     eval(where);
     return false;
   }
}
function newcell(where) {
  var main = document.getElementById("gpshell");
  var clone = main.cloneNode(true);
  searchid(clone,"output").innerHTML = "";
  searchid(clone,"print").innerHTML = "";
  where.parentNode.insertBefore(clone, where);
  searchid(clone,"input").focus();
}
</script>
<script async type="text/javascript" src="./gp-sta.js"></script>
</body>


Collaborative Calculation and Data Science

colby

sagemath

cornellcollege

facebook

github

google

twitter

depaul

Testing latest pari + WASM + node.js... and it works?! Wow.