14 lines
175 B
Text
14 lines
175 B
Text
mata
|
|
|
|
real scalar agm(real scalar a, real scalar b) {
|
|
real scalar c
|
|
do {
|
|
c=0.5*(a+b)
|
|
b=sqrt(a*b)
|
|
a=c
|
|
} while (a-b>1e-15*a)
|
|
return(0.5*(a+b))
|
|
}
|
|
|
|
agm(1,1/sqrt(2))
|
|
end
|