Advanced Search
All Posts : MaplePrimes Commons General Technical Discussions

Please can someone help me figure out what is wrong with this procedure?

 

Thanks.

>  restart:
>  UseHardwareFloats := true:
>  with(SumTools):
>  digits:=10:
>  
>  UniformD:=proc(aN,alb,aub)
>  uses SumTools;
>  local i,j,l,m,F,A,temp,temp1;
>  global H,xroots;
>  H[0](x):=1:
>  for i from 1 to N do
>  F[i]:=A[i,k]*x^k:
>  temp:=sum(F[i],k=0..i):
>  A[i,i]:=1:
>  for m from 1 to i do
>  temp1[m]:=int(temp*H[m-1](x),x=lb..ub):
>  end do:
>  for j from 1 to i do
>       l:=i-j:
>       A[i,l]:=solve(temp1[l+1]=0):
>  end do:
>  H[i](x):=sum(F[i],k=0..i);
>  unapply(H[i](x),x);
>  #xroots[i]:=solve(H[i](x)=0);
>  end do:
>  end proc:
>
>  
>  N:=4:lb:=1.:ub:=10.:
>  UniformD(4,1,10);
>
Error, (in tools/add) too many levels of recursion

 

Please Wait...