In the example below, I find integers that solve the equation. That solution happens to be nonnegative. But in general, just plugging in zeros may not give me a nonnegative solution. What is a good way to find a *nonnegative* integer solution, if one exists?
> aa := x^(-4)*exp(-2*x);
> pmu:=product(muval[j]^a[j],j=1..kmax);
> vars := {a[j]$j=1..kmax};
> solve(identity(pmu=aa,x),vars) assuming x>100;
> isolve(%);
> subs({_Z1=0,_Z2=0},%);
>
This post was generated using the MaplePrimes File Manager
View 85_Jul13.mw on MapleNet or Download 85_Jul13.mw
View file details