次のプログラムは基本的に 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...