These are our functions for f and q:

> #we are going to make a bunch of functions to be used as f's and q's...
#5 -> x
#6 -> y
#7 -> constant = 1
#8 -> harm(5,0,0)
#9 -> harm(.1,0,0)
#10 -> harm(100,0,0)
#11 -> (harm(5,0,0))^2
#12 -> quad-harm(-100,40,-130)
#16 -> tri-harm(-100,40,-130)
# function number 7 should be plugged in for funcno...

make_f_and_q:=proc(T,m,funcno)

make_points(T,m);

harmonic(T,m,1,1,1,funcno);

harmonic(T,m,5,0,0,funcno + 1);

harmonic(T,m,.1,0,0,funcno + 2);

harmonic(T,m,100,0,0,funcno + 3);

mult(T,m,funcno,funcno,funcno + 4);

T[[0],funcno+4]:=0:
T[[0],funcno+5]:=0:
T[[0],funcno+6]:=0:
T[[0],funcno+7]:=-100:
T[[1],funcno+4]:=0:
T[[1],funcno+5]:=0:
T[[1],funcno+6]:=0:
T[[1],funcno+7]:=40:
T[[2],funcno+4]:=0:
T[[2],funcno+5]:=0:
T[[2],funcno+6]:=0:
T[[2],funcno+7]:=-130:
n_harm(T,R,4,m,[0],[1],[2],funcno+5);

T[[0],funcno+8]:=0:
T[[0],funcno+9]:=0:
T[[0],funcno+10]:=-100:
T[[1],funcno+8]:=0:
T[[1],funcno+9]:=0:
T[[1],funcno+10]:=40:
T[[2],funcno+8]:=0:
T[[2],funcno+9]:=0:
T[[2],funcno+10]:=-130:
n_harm(T,R,3,m,[0],[1],[2],funcno+8);
end;