Let Iapp=1.33 Let Ei=0 Let Cs=1 Let Rs=0.001 Input "Cを入力してください。",CtL Input "Rを入力してください。",RtL For t = 1 to 1300 For n=1 to 100 Let a=n*PI() Let c=Cs/CtL DEF f(x)=TAN(x)-c/x DEF g(x)=COS(x)^(-2)+c*x^(-2) LET x=a LET E=0.000001 DO LET x1=x-f(x)/g(x) IF ABS(x1-x)<E THEN EXIT DO LET x=x1 LOOP S=0 S=S+cos(X1)*EXP(-t*X1^2/RtL/CtL)/(X1*((CtL+Cs)*sin(X1)+X1*CtL*cos(X1))) Next n V=Ei+1/1000*(Iapp*(Rs+2*t/(2*CtL+Cs)+RtL*CtL^2*(8*CtL+5*Cs)/(6*(2*CtL+Cs)^2)-2*RtL^2*S)) Print V Next t END
Print"Rs(単位 Ω)=(=γ*RtL)=",Rs Print"Cs(単位 F)=(=θ*CtL)=",Cs Input "CtLを入力してください(単位 F).",CtL Input "RtLを入力してください(単位 Ω).",RtL Print "Emeas/V" For t = 1 to 1300 n=0 While n <=100 n=n+1 Let a=n*PI() Let c=Cs/CtL DEF f(x)=TAN(x)-c/x DEF g(x)=COS(x)^(-2)+c*x^(-2) LET x=a LET E=0.000001 DO LET x1=x-f(x)/g(x) IF ABS(x1-x)<E THEN EXIT DO LET x=x1 LOOP S=0 S=S+cos(x1)*EXP(-t*X1^2/RtL/CtL)/(x1*((CtL+Cs)*sin(x1)+x1*CtL*cos(x1))) Wend Emeas=Ei+1/1000*(Iapp*(Rs+2*t/(2*CtL+Cs)+RtL*CtL^2*(8*CtL+5*Cs)/(6*(2*CtL+Cs)^2)-2*RtL^2*S)) Print a Next t END
>>234 DEF f(x) = x^2 - 2 DEF df(x) = 2*x LET EPS = 0.000001 LET IMAX = 1000 LET x = 2 FOR i=1 TO IMAX LET x1 = x - f(x) / df(x) IF ABS(x1 - x) < EPS THEN EXIT FOR LET x = x1 NEXT i IF i < IMAX THEN PRINT x1 END
236 名前:デフォルトの名無しさん [2007/02/02(金) 23:51:26 ]
10CLEAR :CLS 20DIM A(2),B(2) 30FOR J=0 TO 2 40A(J)= RND 6-1 50NEXT 60*MAIN 70LOCATE 0,0 80FOR J=0 TO 2 90IF B(J)=0 LET A(J)=(A(J)+1) MOD 6 100PRINT MID$ ("$*\7+&",A(J)+1,1); 110NEXT 120C= ASC INKEY$ -49 130IF C>-1 AND C<3 IF B(C)=0 LET D=D+1,B(C)=1 140IF D<3 GOTO *MAIN 150IF A(0)=A(1) AND A(1)=A(2) PRINT " GOOD" 160CALL 48381 170GOTO 10
>>238 100 CLS 3:DIM D(6):FOR I=1 TO 6:D(I)=0:NEXT I:RANDOMIZE 10 110 FOR I=1 TO 100:N=INT(RND*6)+1:D(N) = D(N) + 1:NEXT I 120 LINE (100,10)-(100,300),7:LINE(100,300)-(580,300),7 130 FOR I=1 TO 6:LOCATE 11+10*I,19:PRINT I;:NEXT I 140 FOR I=0 TO 6:LOCATE 9,18-I*3:PRINT I*5;:NEXT I 150 PSET(180,300-D(1)*10),7 160 FOR I=1 TO 6 170 LINE-(I*80+100,300-D(I)*10),7 180 LOCATE I*10+10,18-INT(D(I)*3/5):PRINT D(I); 190 NEXT I 200 LOCATE 0,20
10 REM TAIPINGU GE-MU 20 A$ = "RENSYUU SHITAI MOZIRETSU" 30 PRINT A$ 40 INPUT B$ 50 IF A$ = B$ THEN 80 60 PRINT "TIGAIMASU" 70 GOTO 30 80 PRINT "OWARI" 90 END
100 'NOURYOKU TORE-NINGU 110 SCREEN 0,1:CONSOLE ,,0,1:CLS 3:DEFINT A-Z 130 SIZ = 9:MAX = 100 : WAITCOUNT =10000 140 RANDOMIZE VAL(RIGHT$(TIME$,2)) 150 FOR I=0 TO SIZ 160 Q(I) = RND(MAX):ANS = ANS + Q(I) 170 NEXT 180 FOR I=0 TO SIZ 190 PRINT Q(I) 200 GOSUB *SWAIT 210 CLS 220 NEXT 230 INPUT "GOUKEI";G$ 240 IF VAL(G$) = ANS THEN PRINT "SEIKAI!":END 250 PRINT "ZANNEN!":END 300 *SWAIT 310 FOR S=0 TO WAITCOUNT :NEXT 320 RETURN
>>241 10 REM 4,6で左右移動 20 CONSOLE 0,14,0,1:CLS 3:X=320:Y=199:K=0:W=0 30 A$=INKEY$ 40 IF A$="4" THEN X=X-1:IF X<0 THEN X=0 50 IF A$="6" THEN X=X+1:IF X>639 THEN X=639 60 FOR I=1 TO 10:PSET(RND(1)*640,399),3:NEXT 70 ROLL 1 'ActiveBasicの場合は ROLL 0,-1 80 IF POINT(X,Y)=3 THEN PRINT "走行距離";K;"m":END ELSE K=K+1 90 PSET(X,Y),1:FOR I=1 TO W:NEXT:GOTO 30
100 REM じゃんけんゲーム 110 CLEAR:DIM A$(3),B$(3),H(3,3),S(3) 120 A$(1)="グー":A$(2)="チョキ":A$(3)="パー":S(1)=0:S(2)=0:S(3)=0 130 B$(1)="あいこ":B$(2)="あんたの勝ち":B$(3)="あんたの負け" 140 H(1,1)=1:H(1,2)=2:H(1,3)=3:H(2,1)=3:H(2,2)=1:H(2,3)=2:H(3,1)=2:H(3,2)=3:H(3,3)=1 150 CLS:PRINT "じゃんけんゲーム" 160 INPUT "(1:グー, 2:チョキ, 3:パー, 0:やめる)"; N1 170 IF N1=0 THEN GOTO 220 ELSE IF N1>=1 AND N1<=3 THEN GOTO 190 180 GOTO 160 190 N2=INT(RND(1)*3)*Z+N1-1:IF N2<1 THEN N2=3 200 K=H(N1,N2):S(K)=S(K)+1 210 PRINT "あんた:";A$(N1), "俺:"A$(N2)," ";B$(K):PRINT:GOTO 160 220 PRINT S(2);"勝", S(3);"敗", S(1);"引き分け":END
>>269 フランスとチェコの国旗(N88互換BASIC用) 10 CLS 20 LINE(15,100)-(115,300),1,BF 30 LINE(115,100)-(215,300),7,BF 40 LINE(215,100)-(315,300),2,BF 50 FOR Y=100 TO 200 60 XM=225+Y 70 LINE(325,Y)-(XM,Y),1 80 LINE(XM,Y)-(625,Y),7 90 NEXT Y 100 FOR Y=200 TO 300 120 XM=625-Y 120 LINE(325,Y)-(XM,Y),1 130 LINE(XM,Y)-(625,Y),2 140 NEXT Y 150 END
>>272 式を生成するんじゃなくて数値を計算するんだよね? こんな感じでどうかな。 n = 1 Do Until n = 0 Input "n (0で終了)"; n If n<> 0 Then Print A(n) Loop End Function A(n) If n = 1 Then A = Sqr(1/2) Else A = A(n-1) * Sqr(B(n)) End If End Function Function B(n) If n = 1 Then B = Sqr(1/2) Else B = 1/2 + 1/2 * B(n-1) End If End Function
maxcnt=0 theninit=0 for init=1 to 99999 n=init cnt=0 *loop if n=1 goto *out reminder=n mod 2: cnt=cnt+1 if reminder=1 then n=n*3+1: cnt=cnt+1: print n: goto *loop if reminder=0 then n=n/2: cnt=cnt+1: print n: goto *loop *out if maxcnt<cnt then maxcnt=cnt: theninit=init next print "計算処理回数がもっとも多い値="; theninit
100 PRINT 2 110 PRINT 3 120 N=5 130 IF (N MOD 3)=0 THEN GOTO 190 140 I=6 :M=INT(SQR(N))+1 150 IF I>M THEN PRINT N :GOTO 190 160 IF (N MOD (I-1))=0 THEN GOTO 190 170 IF (N MOD (I+1))=0 THEN GOTO 190 180 I=I+6 :GOTO 150 190 N=N+2 :GOTO 130
10 CLS 20 FOR I = 64 TO 576 STEP 64 30 X=I+12-24*RND : H = 50+RND(1)*100 : R = 50+RND(1)*50 40 CL = 1+RND(1)*7 80 LINE (X,400)-(X,H),CL 90 DR = R/3: R = DR : CF = 1+RND(1)*7 110 CIRCLE (X,H),R,CF 120 R = R + DR : CF = 1+RND(1)*7 130 CIRCLE (X,H),R,CF 150 R = R + DR : CF = 1+RND(1)*7 160 CIRCLE (X,H),R,CF 170 NEXT I 180 END
>>334 こういうことかな 1 rem 行列変数準備 2 dim a(3,3) 3 dim b(3,3) 4 dim c(3,3) 5 dim d(3,3) 6 dim e(3,3) 7 rem 行列A読み込み 8 for i = 1 to 3 9 for j = 1 to 3 10 read a(i,j) 11 next j 12 next i 13 rem 行列B読み込み 14 for i = 1 to 3 15 for j = 1 to 3 16 read b(i,j) 17 next j 18 next i
長いので続く
336 名前:デフォルトの名無しさん [2009/01/21(水) 16:42:21 ]
>>335 続き 19 rem 行列計算 20 for i = 1 to 3 21 for j = 1 to 3 22 c(i,j) = a(i,j) + b(i,j) 23 d(i,j) = a(i,j) - b(i,j) 24 e(i,j) = a(i,j) * b(i,j) 25 next j 26 next i 50 rem 行列A-data 51 data 1,2,3 52 data 4,5,6 53 data 7,8,9 60 rem 行列B-data 61 data 1,2,3 62 data 4,5,6 63 data 7,8,9
19 rem 行列計算 20 for i = 1 to 3 21 for j = 1 to 3 22 c(i,j) = a(i,j) + b(i,j) 23 d(i,j) = a(i,j) - b(i,j) 24 for k = 1 to 3 25 e(i,j) = e(i,j) + a(i,k) * b(k,j) 26 next k 27 next j 28 next i
1 p = 3 2 q = 3 3 k_hit = 100 4 m_max = 100 5 n_max = 100 6 for m = 0 to m_max 7 for n = 0 to n_max 8 k = m * p + n * q 9 if k <= k_hit then print "k=";k;" m=";m;" n=";n 10 next n 11 next m
342 名前:デフォルトの名無しさん [2009/01/21(水) 18:35:13 ]
>>341 微修正 2 q=7 9 if k=> 1 and k <= k_hit then print "k=";k;" m=";m;" n=";n
1 p = 3 2 q = 7 3 k_hit = 100 4 m_max = 100 5 n_max = 100 6 dim v$(100) 7 for m = 0 to m_max 8 for n = 0 to n_max 9 k = m * p + n * q 10 if k => 1 and k <= k_hit then v$(k) = "(";m;",";n;") " 11 next n 12 next m 13 for k = 1 to k_hit 14 print "k=";k;" ";v$(k) 15 next
345 名前:デフォルトの名無しさん [2009/01/21(水) 19:35:34 ]
>>344 6 dim v(100) 10 if k => 1 and k <= k_hit then v(k) = v(k) + 1 14 print "k=";k;" ";v(k)
>>353 十進数BASIC用 100 LET X=5 110 LET E=1E-5 120 DO 130 LET X0=X 140 LET FX0=X^4-2.03790*X^3-15.4245*X^2+15.6696*X+35.4936 150 LET FDX=4*X^3+3*( -2.03790) *X^2+2*( -15.4245) *X+15.6696 160 LET X=X0-FX0/FDX 170 LET FX=X^4-2.03790*X^3-15.4245*X^2+15.6696*X+35.4936 180 PRINT X ;FX 190 LET DFX = ABS(FX-FX0) 200 LOOP WHILE DFX > E 210 END