clear
clc
u=solve('x-b*x^2*log(x/a)=0','x');

[a,b] = meshgrid(1:1:10);
k=subs(u);
contour3(a,b,k,30)
xlabel('a')
ylabel('b')
zlabel('k')

k=subs(u,'[a,b]',[7 5]);
a=2.2;
b=5;
while k>2
    a=a-.005;
    k=subs(u,'[a,b]',[a,b])
end
fprintf('desirable value : a=%6.2f and b=%6.2f \n',a,b)