function dy=g1alin(x,y)

k=input('Introduceti nr k= ');
g=input('Introduceti nr g= ');

a=k+(1-k)/g;
m=(1-k)*(1-a);
n=a*(1-k);
p=k*(1-a);
sc=log((1-k)/(1-a))-(k/(1-k))*log(a/k);
%c=n/p
%c=Momentul critic







dy=(m*y^2+n*y)/(n-p*y);

