ceptreee

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

def ellipse(ax,x,y,width,height,angle=0,color="k"):
    e = patches.Ellipse(xy = (x,y), width = width, height = height,
                    angle = angle, color = color)
    ax.add_patch(e)

def plot_rw():

ceptreee

hai1 = "(｣・ω・)｣ﾊｲ"
hai2 = "(｣・ω・)｢ﾊｲ"
hai3 = "(｢・ω・)｢ﾊｲ"
hai4 = "(｢・ω・)｣ﾊｲ"

do for [i=1:100]{ 
    if (i % 4 == 0) {
        set title hai1 font "Arial, 24"
        plot x 
    }

ceptreee

hai1 = "(｣・ω・)｣ﾊｲ"
hai2 = "(｣・ω・)｢ﾊｲ"
hai3 = "(｢・ω・)｢ﾊｲ"
hai4 = "(｢・ω・)｣ﾊｲﾎﾞｰﾙ"

do for [i=1:100]{ 
    if (i % 4 == 0) {
        set title hai1 font "Arial, 24"
        plot x 
    }

ceptreee

a = "(｣・ω・)｢"
b = "(｢・ω・)｣"

do for [i=1:100]{ 
    if (i % 2 == 0) {
        set title a font "Arial, 24"
        plot x 
    } else {
        set title b font "Arial, 24"
        plot x**2

ceptreee

ceptreee

ceptreee

include("LinearVolterraIntegralEquation.jl")

using PyPlot 

problem = "1"
Lt,K,g,y_exc = Problems_LVIE(problem)

N_NI = 3

Nm = 3

ceptreee

ceptreee

ceptreee

function hoge(n)
    if n == 1
        Lt    = 1.0
        K     = (t,s) -> 2.0 * cos( t - s )
        g     = (t) -> exp(t)
        y_exc = (t) -> exp(t) * ( 1.0 + t )^2 
    elseif n == 2
        Lt    = 3.0
        K     = (t,s) -> (t-s)
        g     = (t) -> 1.0 - t - t^2/2.0