function Gat()
fplot(@(x)F(x,10,10,100,0.1,0.1),[0 100])
fminbnd(@(x)-F(x,10,10,100,0.1,0.1),0,100)

function y=F(x,I0,k,V,r,p)
    y=-(I0+k*x)-exp(-p*x)*V/r;
