PRINT
PRINT "３点を通る放物線"
PRINT "     y = a x^2 + b x + c"
PRINT
INPUT PROMPT "第１点 （ x1 ， y1 ）  ":x1,y1
INPUT PROMPT "第２点 （ x2 ， y2 ）  ":x2,y2
INPUT PROMPT "第３点 （ x3 ， y3 ）  ":x3,y3
PRINT
DIM Z(3,3),J(3),K(3,3),L(3)
LET Z(1,1)=x1^2
LET Z(1,2)=x1
LET Z(1,3)=1
LET Z(2,1)=x2^2
LET Z(2,2)=x2
LET Z(2,3)=1
LET Z(3,1)=x3^2
LET Z(3,2)=x3
LET Z(3,3)=1
LET J(1)=y1
LET J(2)=y2
LET J(3)=y3
IF DET(Z)=0 THEN GOTO 10
MAT K=INV(Z)
MAT L=K*J
IF L(1)=0 THEN GOTO 10
PRINT "y = (";ROUND(L(1),5);") x^2 + (";ROUND(L(2),5);") x + (";ROUND(L(3),5);")"
PRINT
PRINT "y1 : ";y1;"=";ROUND(L(1)*x1^2+L(2)*x1+L(3),5)
PRINT "y2 : ";y2;"=";ROUND(L(1)*x2^2+L(2)*x2+L(3),5)
PRINT "y3 : ";y3;"=";ROUND(L(1)*x3^2+L(2)*x3+L(3),5)
GOTO 20
10 PRINT "そんな２次関数は存在しません"
20 PRINT
END