#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[1],0,
seq(print(x>=Lx[i] and x<Lx[i+1],Lf[i]),i=1..nops(Lf)) ,
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