doraneko94

// This code was inspired by https://github.com/rp-rs/rp-hal/blob/main/rp2040-hal/examples/watchdog.rs

#![no_std]
#![no_main]

use panic_halt as _;
use rp2040_hal as hal;
use hal::pac;

use embedded_hal::digital::v2::OutputPin;

doraneko94

from matplotlib import pyplot as plt
from matplotlib.patches import Ellipse
import numpy as np

fig, ax = plt.subplots(1)
ax.set_aspect('equal')

m = np.array([8, 12])
S = np.array([[16, np.sqrt(78)], [np.sqrt(78), 9]])
r = np.random.multivariate_normal(m, S, size=100) # random points to plot

doraneko94

import numpy as np
import pandas as pd
from scipy import stats

class TwoWayAnovaRM:
    
    def __init__(self, data: pd.DataFrame, subject, samples: list):
        
        self.data = data
        self.subject = subject

doraneko94

doraneko94

import numpy as np

class PolarDiff:

    def __init__(self, n, r2=1, x=None, clip=None, epsilon=1e-7):

        if x is not None:
            if n != len(x):
                raise ValueError
            self.theta = np.zeros(n-1)

doraneko94

# ２種類の気体分子が円形の容器内で動く様子のシミュレーションを行い、
# その様子をアニメーションとして保存する。
# 詳細は https://ushitora.net/archives/2197 を参照

import numpy as np
from matplotlib import pyplot as plt
from matplotlib import animation as animation

""" データの準備 """
# 乱数で２種類の気体の極座標生成

doraneko94

import numpy as np
from matplotlib import pyplot as plt
from scipy.stats import norm

def bic_raw(x, k):
    n = len(x)
    return -2 * np.sum(np.log(norm.pdf(x))) + k * np.log(n)

def bic_rss(x, k):
    n = len(x)

doraneko94

import numpy as np

def tau_num(X: np.ndarray, Y: np.ndarray):
    XpY = X + Y
    XmY = X - Y
    A = (XpY == 2).sum()
    S = (XpY == 0).sum()
    dX = (XmY == 1).sum()
    dY = (XmY == -1).sum()
    return A * S - dX * dY

doraneko94

lst = [11, 10, 22, 10, 25]
ans = 1
tolerance = 1e-4
size = len(lst)

def n_merge(N):
    ans = N
    for i in range(1, N):
        ans *= i**2
    return ans

doraneko94

#include <iostream>
using namespace std;

int N, M, Q;
int a[2010];
long long dp[2010][2010] = {0};
long long rdp[2010][2010] = {0};
long long accum[2010];
const long long mod = 1000000007;