clear

clc

ng$=3$;

beta$=[5\mathrm{.}35\mathrm{.}55\mathrm{.}8]$;

gama$=[0\mathrm{.}0080\mathrm{.}0120\mathrm{.}018]$;

pd$=975$;

sinvgam$=$sum$(1./$gama$)$;

i$=1$;

imax$=10$;

dp$($i$)=100$; la$($i$)=6$;

while$($abs$($dp$)>1$e$-6$&&i