Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- object uklady extends App{
- def gauss(n:Int,ab:Array[Array[Double]]):Array[Double] = {
- var m:Double = 0
- var s:Double = 0
- val wynik:Array[Double]=new Array[Double](n)
- for(i<-0 until n-1){
- for(j<-i+1 until n){
- m= -ab(j)(i)/ab(i)(i)
- for(k<-i to n){
- ab(j)(k)+=m*ab(i)(k)
- }
- }
- }
- wynik(n-1)=ab(n-1)(n)/ab(n-1)(n-1)
- for(i<-n-1 to 0 by -1){
- s=ab(i)(n)
- for(j<-n-1 to i+1 by -1){
- s -= ab(i)(j) * wynik(j)
- wynik(i)=s/ab(i)(i)
- }
- }
- wynik
- }
- var ab = Array.ofDim[Double](3,4)
- ab = Array(Array(2.63,0.15,0.15,0.39), Array(0.115,0.1,0.71,0.42), Array(0.73,0.74,0.1,0.32))
- gauss(3,ab).foreach(x => println(x))
- }
- ------------wynik--------
- 0.10284839476952914
- 0.2582002978866951
- 0.5385245138208938
Add Comment
Please, Sign In to add comment