1040ForShowB10.mws

>    restart;

1040ForShowB10.mws
Thu Nov 10 10:39:51 MST 2005

>    with(plots,display):  with(plottools, disk):

########################################

>    f := sqrt;

f := sqrt

>    #f(4);

>    XRG := -1..3:  YRG := XRG:

>    Pf := plot(f, XRG, YRG, colour=BLACK, thickness=2):

>    display(Pf);

[Maple Plot]

>    A := 1:

>    #Taylor( f(x), x=A, 5);

>    series( f(x), x=A, 5);

series(1+1/2*(x-1)-1/8*(x-1)^2+1/16*(x-1)^3-5/128*(x-1)^4+O((x-1)^5),x=-(-1),5)
series(1+1/2*(x-1)-1/8*(x-1)^2+1/16*(x-1)^3-5/128*(x-1)^4+O((x-1)^5),x=-(-1),5)

>    PD := disk( [A, f(A)], 1/30, colour=BLACK ):

>    #display([Pf,PD]);

Several Taylor Polynomials and how they relate to f

>    N := 10:

>    SEQ := NULL:

>    for k from 1 to N
  do
    series( f(x), x=A, k):
    convert( %, polynom);
    plot(%, x=XRG, y=YRG):
    SEQ := SEQ, %:
  end do:

>    SEQR := display( [SEQ], insequence=true):

>    Enchilada := display([Pf, PD, SEQR]):

>    display(Enchilada, insequence=true);

[Maple Plot]

######################################################

>    binomial(4,2);

6

>    N := 15:

>    for n from 0 to N
  do
     for k from 0 to n
       do
         binomial(n,k);
         printf(" %d ",%);
       end do;
       printf("\n");
  end do;

 1 

 1  1 

 1  2  1 

 1  3  3  1 

 1  4  6  4  1 

 1  5  10  10  5  1 

 1  6  15  20  15  6  1 

 1  7  21  35  35  21  7  1 

 1  8  28  56  70  56  28  8  1 

 1  9  36  84  126  126  84  36  9  1 

 1  10  45  120  210  252  210  120  45  10  1 

 1  11  55  165  330  462  462  330  165  55  11  1 

 1  12  66  220  495  792  924  792  495  220  66  12  1 

 1  13  78  286  715  1287  1716  1716  1287  715  286  78  13  1 

 1  14  91  364  1001  2002  3003  3432  3003  2002  1001  364  91  14  1 

 1  15  105  455  1365  3003  5005  6435  6435  5005  3003  1365  455  105  15  1 

#########################################################

>    G := (x,y,z) -> (x+y+z)/3;

G := proc (x, y, z) options operator, arrow; 1/3*x+1/3*y+1/3*z end proc

>    G(3,2,6);

11/3

>    AVE := proc(a,b,c)
  description "Compute and print the average of a,b,c.";
  local q;

  q := a + b + c;
  q := q/3;
  print(q);
  return(q);
    end:

>    q := 1728;

q := 1728

>    GG := AVE(5,6,7); %;

6

GG := 6

6

>    q;

1728

>    #print(AVE);

>   

>    ##############################################

>   

>   

>   

>