next up previous
Next: 独り言 (何だか長くなっている) Up: A..5 BASIC Previous: A..5 BASIC

十進BASICプログラム例

次のプログラムは基本的に C プログラムと同じことをしている。

computesin.bas

DEF halfcos(COSX)=SQR((1+COSX)/2)
DEF halfsin(COSX)=SQR((1-COSX)/2)
LET  SIN18=(SQR(5)-1)/4
LET  cos18=(SQR((5+SQR(5))/2))/2
LET  SIN15=(SQR(3)-1)/(2*SQR(2))
LET  COS15=(SQR(3)+1)/(2*SQR(2))
LET  SIN3=COS15*SIN18-SIN15*COS18
LET  COS3=COS15*COS18+SIN15*SIN18
LET  SIN1p5=halfsin(COS3)
LET  COS1p5=halfcos(COS3)
LET  SIN0p75=halfsin(COS1p5)
LET  COS0p75=halfcos(COS1p5)
LET  SIN1=SIN0p75*4/3
PRINT "sin(18)=";SIN18
PRINT "cos(18)=";COS18
PRINT "sin(15)=";SIN15
PRINT "cos(15)=";COS15
PRINT "sin(3)=";SIN3
PRINT "cos(3)=";COS3
PRINT "sin(1)=";SIN1
END

[オプション]→[数値]で「10進1000桁」を選択して実行すると、 以下の結果が得られる。


sin(18)= .3090169943749474241022934171828190588601545899028814310677243113526302314094512248536036020946955687...
cos(18)= .9510565162951535721164393333793821434056986341257502224473056444301531700851935017187928109708113816...
sin(15)= .2588190451025207623488988376240483283490689013199305138140032073150569747488019969223679746942496655...
cos(15)= .9659258262890682867497431997288973676339048390084045504023430763104232139798555163475617418580704518...
sin(3)= .0523359562429438327221186296090784187310182539401649204835093815998571046417545468644645988118869398...
cos(3)= .9986295347545738737844920584394365805909522907677855324414412548316489737334783186353320282332754504...
sin(1)= .0174527940951259202537122796038029453580807447073735925504530431057687096100067416295782365116468161...


next up previous
Next: 独り言 (何だか長くなっている) Up: A..5 BASIC Previous: A..5 BASIC
Masashi Katsurada
平成17年7月19日