please i want to write the maple programme  for solving y^/ =y, y(0)=1, 0≤x≤1 with h=0.1 using the following three schemes simultaneously:

1. (7/12)y_n+1 - y_n+1/2 +5/12y_n =(h/12)f_n+2

2. y_n+2 - (4/3)y_n+1 +1/3y_n =(2h/3)f_n+2

3. y_n+1 -y_n =hf_n+1/2

ifplease how can i factorise in powers of (x-x_n), the expression  y(t):=(1/3h²)[(-x+x_n+h)(-x+x_n+3h)]y_n+(1/3h²)[(-x_n+x)(-x+x_n+4h)]y_n+1 -(1/3h)[(-x_n+x)(-x+x_n+3h)]_fn+2?

i was trying to familiarize with maple 13 but encountered this. pls how may i go out from here?

> with(Student[NumericalAnalysis]);
> DE1 := d*y(t)/dt = y(t)-t^2+1;
                           d y(t)           2   
        ...