2119.sagews

Test-berekening voor het project van Bart en Marcel

(a,b)=(2,1)

x = var('x')

f0(x) = 1/2; f0

x |--> 1/2

def F(x,f): return a - b*(x-a)/(f(x)-b)

def plotF():
    f1(x) = F(x,f0)
    f2(x) = F(x,f1)
    f3(x) = F(x,f2)
    f4(x) = F(x,f3)
    f5(x) = F(x,f4)
    f6(x) = F(x,f5)
    d1(x) = 1/2*(a+b+sqrt((a+b)^2-4*x*b))
    d2(x) = 1/2*(a+b-sqrt((a+b)^2-4*x*b))
    xaxis = (x,-10,10)
    p  = plot(f0(x), xaxis, detect_poles='show', rgbcolor='#000000') #zwart
    p += plot(f1(x), xaxis, detect_poles='show', rgbcolor='#0000ff') #blauw
    p += plot(f2(x), xaxis, detect_poles='show', rgbcolor='#00ff00') #groen
    p += plot(f3(x), xaxis, detect_poles='show', rgbcolor='#ff0000') #rood
    p += plot(f4(x), xaxis, detect_poles='show', rgbcolor='#00ffff') #cyaan
    p += plot(f5(x), xaxis, detect_poles='show', rgbcolor='#ff00ff') #magenta
    p += plot(f6(x), xaxis, detect_poles='show', rgbcolor='#ffff00') #geel
    p += plot(d1(x), xaxis, detect_poles='show', rgbcolor='#000000') #zwart
    p += plot(d2(x), xaxis, detect_poles='show', rgbcolor='#000000') #zwart
    p.set_aspect_ratio(1)
    show(p, xmin = -10, xmax = 10, ymin=-10, ymax = 10)

(a,b)=(1,2); f0(x) = 1; plotF()

verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 88 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: '' verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 89 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: ''

(a,b)=(2,1); f0(x) = 1/2; plotF()

verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 77 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: '' verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 77 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: ''

(a,b)=(1,2); f0(x) = x; plotF()

verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 88 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: '' verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 89 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: ''

(a,b)=(2,1); f0(x) = x; plotF()

verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 78 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: '' verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 77 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: ''

(a,b)=(2,1); f0(x) = x^2; plotF()

verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 78 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: '' verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 78 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: ''

(a,b)=(2,1); f0(x) = gamma(x); plotF()

verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 78 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: '' verbose 0 (3495: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 77 points. verbose 0 (3495: plot.py, generate_plot_points) Last error message: ''