function parms = GetParms(Ww,Lfd,Lo)
	%% Use lazy programmer's optimization to find best-fit parameters for inductance network *)
  
	% initial guess for model parameters
	q=(1:8)/2;
  
	% Error at initial guess for model parameters
	b = GetErr(Ww,Lfd,Lo,q);
  
	% Step size
	d = 0.1;
  
	% Counter for number of steps since progress was made
	n = 0;
    
    for m=0:100000
        qp=q;
        for k = 1:length(q)
            qp(k) = qp(k) + d*(2*rand - 1);
        end
        c = GetErr(Ww,Lfd,Lo,qp);
        if (c < b)
            n = 0;
            b = c;
            q = qp;
        else
            n=n+1;
        end
        if (n == 100)
            n = 0;
            d = d/2; 
%            fprintf('%i  %g\n',d,b);
        end
        if (d < 1e-5)
            break;
        end
    end
    
    L=zeros(1,5); L(1)=Lo;
    R=zeros(1,4);
    T=zeros(1,8);
    for k=1:8
        T(k) = 2*pi*10^q(k);
    end
    for k=1:4
        R(k)   = T(2*k-1)*L(k);
        L(k+1) = R(k)/T(2*k);
    end

    parms = [L,R];

  end