>
|
use plots in
module()
local g,a,b,c,P1,P2,p1,p2,p3,p4,p5,p6,p7;
g := x-> 2+sin(x);
a := .5; b := -1.;c:=.1;
P1:=.2347601860*x^2+1.076231171*x+2.;
P2:=-.2347601860*x^2+1.076231171*x+2.;
p1:=plot([[x,P1,x =-1..0],[x,P2,x=0..1],[x,g(x),x=-1..1],[x,0,x=-1..1]],color=[red,green,black,black], thickness=[8,8,1,1],transparency=[.7,.7,0,0],scaling=constrained,axes=none,view=[-1..1,-.2..3]);
p2:=textplot({[-1,-.1,x[k]],[0,-.1,x[k+2]],[1,-.1,x[k+4]],[-.5,-.1,x[k+1]],[.5,-.1,x[k+3]]});
p3:=plot([[[-1,g(-1)],[-1,0]],[[-.5,g(-.5)],[-.5,0]],[[0,0],[0,g(0)]]],color=[red$3],thickness=[1,1,4]);
p4:=plot([[[.01,0],[.01,g(.01)]],[[.5,0],[.5,g(.5)]],[[1,0],[1,g(1)]]],color=[green$3],thickness=[2,1,1]);
p5:=plot([[-1,g(-1)],[-.5,g(-.5)],[0,g(0)],[.5,g(.5)],[1,g(1)]],style=point,symbol=solidcircle,symbolsize=15,color=black):
p6:=textplot({[-1,g(-1)+c,typeset(``(x[k],f[k]))],[-.5,g(-.5)+c,typeset(``(x[k+1],f[k+1]))],[0,g(0)+c,typeset(``(x[k+2],f[k+2]))],[.5,g(.5)+c,typeset(``(x[k+3],f[k+3]))],[1,g(1)+c,typeset(``(x[k+4],f[k+4]))]},align=left);
p7:=display(p1,p2,p3,p4,p5,p6);
print(p7);
end module:
end use:
|