#
#
p=0.25;
for i=1:100
 deltap=(-2+0.5/p^0.2+0.5/p^0.5)/(0.1/p^1.2+0.25/p^1.5);
 p=deltap+p;
if  abs(deltap) < 1.e-8, break, end
end
disp(p);
