function [ynew,dyda] = fit_jenn(x,A)

k=4*pi/1e9;
q= [430 70 42];

ynew = zeros(length(x),1);
dyda = zeros(length(x),3);
dyda(:,2) = q(2)./x(:,2);
dyda(:,3) = q(3)./x(:,3);

dyda(:,1) = q(1)./x(:,1);
ynew = (q(1)*A(1)./x(:,1)+q(2)*A(2)./x(:,2)+q(3)*A(3)./x(:,3))./k;