QR decomposition: Difference between revisions

[double[]]$u = $R[$h..$pm] | foreach{$_[$h]}

[double]$nu = $u | foreach {[double]$sq = 0} {$sq += $_*$_} {[Math]::Sqrt($sq)}

[double]$nu = $u | foreach {$sq = 0} {$sq += $_*$_} {[Math]::Sqrt($sq)}

[double[]]$u = $u | foreach { $_/$nu}