function solve_n_eqs
clear all;clc
d=45;
D=150;
n=5000;   
x=linspace(-0.0005,0.001,n);
y=0.1248*x.^4.*( 1./(pi*d-0.2*x)-1./(pi*D-0.2*x)   )+0.2496*x.^3.*log(abs( (pi*D-0.2*x)./(pi*d-0.2*x)  ))+1.248*pi*x.^2*(D-d)-0.2*x.^3*(D^3-d^3)    ;
figure(1)
plot(x,y,'b',[min(x) max(x)],[0 0],'r'),axis([min(x) max(x) min(y) max(y)]);

format long

[x fval]=fsolve(@(x) 0.1248*x.^4.*( 1./(pi*d-0.2*x)-1./(pi*D-0.2*x)   )+0.2496*x.^3.*log(abs( (pi*D-0.2*x)./(pi*d-0.2*x)  ))+1.248*pi*x.^2*(D-d)-0.2*x.^3*(D^3-d^3),[1e-5 1e-3])