連立方程式

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26

#module
    //2×2の行列式の計算マクロ
    #define ctype Determinant2(%1,%2,%3,%4)  ((%1*%4)-(%2*%3))
    //3×3の行列式の計算マクロ
    #define ctype Determinant3(%1,%2,%3,%4,%5,%6,%7,%8,%9)  ((%1*%5*%9)+(%2*%6*%7)+(%3*%4*%8)-(%3*%5*%7)-(%2*%4*%9)-(%1*%6*%8))
 
    #deffunc SimulEq2 array  sol , array exp1,array exp2
        a1 = double(exp1(0)) : b1 = double(exp1(1)) : c1 = double(exp1(2))
        a2 = double(exp2(0)) : b2 = double(exp2(1)) : c2 = double(exp2(2))
        D = Determinant2(a1,b1,a2,b2)
        if (D==0):return 0
        sol(0) = Determinant2(c1,b1,c2,b2)/D
        sol(1) = Determinant2(a1,c1,a2,c2)/D
        return 1
 
    #deffunc SimulEq3 array  sol , array exp1,array exp2,array exp3
        a1 = double(exp1(0)) : b1 = double(exp1(1)) : c1 = double(exp1(2)) :d1 = double(exp1(3))
        a2 = double(exp2(0)) : b2 = double(exp2(1)) : c2 = double(exp2(2)) :d2 = double(exp2(3))
        a3 = double(exp3(0)) : b3 = double(exp3(1)) : c3 = double(exp3(2)) :d3 = double(exp3(3))
        D = Determinant3(a1,b1,c1,a2,b2,c2,a3,b3,c3)
        if (D==0):return 0
        sol(0) = Determinant3(d1,b1,c1,d2,b2,c2,d3,b3,c3)/D
        sol(1) = Determinant3(a1,d1,c1,a2,d2,c2,a3,d3,c3)/D
        sol(2) = Determinant3(a1,b1,d1,a2,b2,d2,a3,b3,d3)/D
        return 1
#global

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
-
|
|
-
|
|
!
 
 

#include "Equation.hsp"
#packopt name "21Equation"
title "2元1次連立方程式を解くプログラム"
width 450,320
pos 0,25
font msgothic,125
print "{"
ddim  exp1,3:ddim  exp2,3:ddim sol,2
font msgothic,16
objmode 2
navit="x+","y="
repeat 3
    pos 75+130*cnt,50
    Input  exp1(cnt),100,30
    pos 180+130*cnt,50
    if cnt<2:print navit(cnt)
    loop
 
repeat 3
    pos 75+130*cnt,100
    Input  exp2(cnt),100,30
    pos 180+130*cnt,100
    if cnt<2:print navit(cnt)
    loop
 
pos 5,200
objsize 200,30
button  "解く",*Solve
pos 5,240:print "x=" : pos 30,240 : Input Sol(0),100,30
pos 5,280:print "y=" : pos 30,280 : Input Sol(1),100,30
stop
*Solve
 
 SimulEq2 Sol,exp1,exp2
 if stat{
     objprm 7,str(Sol(0))
     objprm 8,str(Sol(1))
 }else{
     objprm 7,"解なし"
     objprm 8,"解なし"
 }
 
 stop

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
-
|
|
|
|
-
|
|
|
|
!
 
 

#include "Equation.hsp"
#packopt name "31Equation"
title "3元1次連立方程式を解くプログラム"
width 600,370
pos 0,30
font msgothic,175
print "{"
ddim  exp1,4:ddim  exp2,4:ddim  exp3,4:ddim sol,3
font msgothic,16
objmode 2
navit="x+","y+","z="
repeat 4
    pos 75+130*cnt,50
    Input  exp1(cnt),100,30
    pos 180+130*cnt,50
    if cnt<3:print navit(cnt)
    loop
 
repeat 4
    pos 75+130*cnt,100
    Input  exp2(cnt),100,30
    pos 180+130*cnt,100
    if cnt<3:print navit(cnt)
    loop
repeat 4
    pos 75+130*cnt,150
    Input  exp3(cnt),100,30
    pos 180+130*cnt,150
    if cnt<3:print navit(cnt)
    loop
 
 
pos 5,210
objsize 200,30
button  "解く",*Solve
pos 5,250:print "x=" : pos 30,250 : Input Sol(0),100,30
pos 5,290:print "y=" : pos 30,290 : Input Sol(1),100,30
pos 5,330:print "z=" : pos 30,330 : Input Sol(2),100,30
 
stop
*Solve
 
 SimulEq3 Sol,exp1,exp2,exp3
 if stat{
     objprm 13,str(Sol(0))
     objprm 14,str(Sol(1))
     objprm 15,str(Sol(2))
     
 }else{
     objprm 13,"解なし"
     objprm 14,"解なし"
     objprm 15,"解なし"
     
 }
 
 stop