mskashi

import std/sequtils
import std/strutils
import std/tables

type board = array[1..19, array[1..19, char]]

var table = initTable[string, int]()
var count:int = 0

proc print(a:board) =

mskashi

# calculate coefficients of Newton-Cotes integral formula
#  for \int_a^b f(x)dx
#  using sampling points x0(=a), x1, ..., xn-1, xn(=b)

import sympy as sym

n = 7;

x = sym.Symbol('x');
a = sym.Symbol('a');

mskashi

#include <stdio.h>
#include <float.h>
#include <stdint.h>
#include <inttypes.h>

int64_t d2i(double x)
{
	union {
		int64_t i;
		double d;

mskashi

#include <iostream>
#include <cfenv>
#include <cmath>

int check_rounding()
{
	volatile double x, y, z;

	x = 1;
	y = pow(2., -55);

mskashi

#include <kv/kraw-approx.hpp>
#include <kv/defint.hpp>

namespace ub = boost::numeric::ublas;

typedef kv::interval<double> itv;

struct Func {
	template <class T> ub::vector<T> operator() (const ub::vector<T>& x) {
		static T p = kv::constants<T>::pi();