#Sowmya Srinivasan hw23.txt 17 April 2014 #Pw(Lf,x,Lx) should translate a piecewise function from list format to piecewise Pw:=proc(Lf,x,Lx) local i,a1,a2,a3 : if nops(Lf)<>nops(Lx)-1 then RETURN(FAIL) : fi : piecewise(x=Lx[i] and x=Lx[-1],0) : end : #Compare(p,q,v,x,x1) Compare:=proc(p,q,v,x,x1) local u,f ,i ,S,L ,N ,k : u:=v*x*(1-x) : f:=-diff(p*diff(u,x),x)+q*u : S:={10,20,30,40} : for N in S do for i from 1 to N-1 do if i/N<=x1 and x1<=(i+1)/N then L:=[op(L),subs(x=x1,RR(p,x,f,x,[seq(i/N,i=0..N)])[i+1])]: fi : fi : od : od : L: end : #This code doesn't work yet because I'm not sure how to choose the right #element of the output of RR