function asdfg
clear
clc
format long
[x0,y0] = ode45(@myfun,[0.1581:1e-5:0.67],0);
figure(1)
plot(y0,x0)
grid on
hold on
t = 1:300;
x = spline(y0,x0,t);
plot(t,x,'r.')
A = 1.515*0.2842*x./(1+0.2842*x);
B = (0.02105*x.^0.9208)./(1-x).^0.9208;
[A',B']
function y = myfun(x,y)
y = 0.0025*(0.9015./(1+0.2842*x)-0.2562*x./(1+0.2842*x).^2+0.2298./(x.^0.0792.*(1-x).^0.9208)...
    +0.2298*x.^0.9208./(1-x).^1.9208+0.017965)./(630*(8208*10^(-15)*3168).*(0.9-x));