kikosiak

Untitled

Dec 3rd, 2019
76
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.82 KB | None | 0 0
  1. object uklady extends App{
  2. def gauss(n:Int,ab:Array[Array[Double]]):Array[Double] = {
  3. var m:Double = 0
  4. var s:Double = 0
  5. val wynik:Array[Double]=new Array[Double](n)
  6. for(i<-0 until n-1){
  7. for(j<-i+1 until n){
  8. m= -ab(j)(i)/ab(i)(i)
  9. for(k<-i to n){
  10. ab(j)(k)+=m*ab(i)(k)
  11. }
  12. }
  13. }
  14. wynik(n-1)=ab(n-1)(n)/ab(n-1)(n-1)
  15. for(i<-n-1 to 0 by -1){
  16. s=ab(i)(n)
  17. for(j<-n-1 to i+1 by -1){
  18. s -= ab(i)(j) * wynik(j)
  19. wynik(i)=s/ab(i)(i)
  20. }
  21. }
  22. wynik
  23. }
  24.  
  25. var ab = Array.ofDim[Double](3,4)
  26. 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))
  27.  
  28. gauss(3,ab).foreach(x => println(x))
  29.  
  30. }
  31.  
  32. ------------wynik--------
  33. 0.10284839476952914
  34. 0.2582002978866951
  35. 0.5385245138208938
Add Comment
Please, Sign In to add comment