AndreyAkinshin

################################################################################
# A simulation of the following system:
#
# dm1/dt = -k1 * m1 + f1(p3)   dp1/dt = mu1 * (m1 - p1)
# dm2/dt = -k2 * m2 + f2(p1)   dp2/dt = mu2 * (m2 - p2)
# dm3/dt = -k3 * m3 + f3(p2)   dp3/dt = mu3 * (m3 - p3)
#
# Parameters: 
#   k1  = 1.5   k2  = 1.7   k3  = 2.0
#   mu1 = 3.0   mu2 = 3.5   mu3 = 4.0

AndreyAkinshin

library(deSolve)
library(dplyr)
library(ggplot2)

# Model
f <- function(x) 258 / (1 + x ^ 7)
model <- function(t, state, parameters) {
  with(as.list(c(state)), {
    dx1 <- f(x5) - x1
    dx2 <- f(x1) - x2

AndreyAkinshin

################################################################################
# A simulation of the following system:
#
# dm1/dt = -k1 * m1 + f1(p5)   dp1/dt = mu1 * m1 - nu1 * p1
# dm2/dt = -k2 * m2 + f2(p1)   dp2/dt = mu2 * m2 - nu2 * p2
# dm3/dt = -k3 * m3 + f3(p2)   dp3/dt = mu3 * m3 - nu3 * p3
# dm4/dt = -k4 * m4 + f4(p3)   dp4/dt = mu4 * m4 - nu4 * p4
# dm5/dt = -k5 * m5 + f5(p4)   dp5/dt = mu5 * m5 - nu5 * p5
#
# Parameters: 

AndreyAkinshin

################################################################################
# A simulation of the following system:
#
# dm1/dt = -k1 * m1 + f1(p9)   dp1/dt = mu1 * m1 - nu1 * p1
# dm2/dt = -k2 * m2 + f2(p1)   dp2/dt = mu2 * m2 - nu2 * p2
# dm3/dt = -k3 * m3 + f3(p2)   dp3/dt = mu3 * m3 - nu3 * p3
# dm4/dt = -k1 * m1 + f1(p3)   dp4/dt = mu1 * m1 - nu1 * p4
# dm5/dt = -k2 * m2 + f2(p4)   dp5/dt = mu2 * m2 - nu2 * p5
# dm6/dt = -k3 * m3 + f3(p5)   dp6/dt = mu3 * m3 - nu3 * p6
# dm7/dt = -k1 * m1 + f1(p6)   dp7/dt = mu1 * m1 - nu1 * p7

AndreyAkinshin

################################################################################
# A simulation of the following system:
#
# dx1/dt =          - k1 * x1 + f1(x9)
# dx2/dt = mu2 * x1 - k2 * x2
# dx3/dt = mu3 * x2 - k3 * x3
# dx4/dt = mu4 * x3 - k4 * x4
# dx5/dt =          - k5 * x5 + f5(x4)
# dx6/dt = mu6 * x5 - k6 * x6
# dx7/dt = mu7 * x6 - k7 * x7

AndreyAkinshin

################################################################################
# A simulation of the following system:
#
# dm1/dt = L1(p3) - k1 * m1
# dm2/dt = L2(p1) - k2 * m2
# dm3/dt = L3(p2) - k3 * m3
# dp1/dt = G1(m1) - l1 * p1
# dp2/dt = G2(m2) - l2 * p2
# dp3/dt = G3(m3) - l3 * p3
#

AndreyAkinshin

################################################################################
# A simulation of the following system:
#
# dx/dt = L1(z) - k1 * x
# dy/dt = L2(x) - k2 * y
# dz/dt = L3(y) - k3 * z
#
# Parameters:
#   k1 = 1.2  k2 = 1.4  k3 = 1.5
#   a1 = 2    a2 = 3    a3 = 4

AndreyAkinshin

################################################################################
# A simulation of the following system:
#
# dx/dt = L1(z) - k1 * x
# dy/dt = L2(x) - k2 * y
# dz/dt = L3(y) - k3 * z
#
# Parameters:
#   k1 = 1.1  k2 = 2.2  k3 = 3.3
#   a1 = 2    a2 = 3    a3 = 4

AndreyAkinshin

// Copyright (c) 2019 Andrey Akinshin
// Licensed under The MIT License https://opensource.org/licenses/MIT
using System;
using System.Collections.Generic;
using System.Linq;

/// <summary>
/// The ED-PELT algorithm for changepoint detection.
///
/// <remarks>

AndreyAkinshin

using System.Reflection;
using Xunit;

namespace AssemblyNameTests
{
    public class ParsingTests
    {
        [Theory]
        [InlineData("  Por ta bility")]
        [InlineData("'Por, tability!'")]