常連

### 前提

以下の数列がある。
500,497,494,491,・・・・・
この数列の初項から第 100 項まで表示して、かつこの数列の合計（和）を表示
するプログラムを作成せよ。

このプログラムの書き方が分かりません。while文で表示するまではかけたのですが、かっこの数列の合計の仕方が分かりません。JavaScriptで分かる方がいたら教えてください。
プログラムを学び始めて2週間の初心者で至らぬ点も多いと思いますが、お願いします・・・。

### 実現したいこと

while文の数列の合計(和)と表示
※while文を使ってお願いします。出来れば下記のソースに「これを付け足せばok」とかだと非常に助かります。

### 発生している問題・エラーメッセージ


### 該当のソースコード
```
<html>
<head>
</head>
<body>
<script type = "text/JavaScript">

	var seisuu = -1
	var seisuu2 = 0
	var keika

	while ( seisuu < 100 )
	{
		seisuu = seisuu + 1
		seisuu2 = seisuu * 3
		keika = 500 - seisuu2
		document.write(keika,"<br>");
	}
	
</script>
</body>
</html>
```
### 試したこと

ここに問題に対して試したことを記載してください。

### 補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。


while文の質問。「この数列の初項から第 100 項まで表示して、かつこの数列の合計（和）を表示 するプログラムを作成せよ。」の解決方法。

### 実現したいこと

whileを使用した無限ループ中にトピックを受信したい

### 前提

UnityとROSの通信について勉強中です。
元々ROSのTkinterで出来ていたことをUnityで操作したくなり、結果苦戦しております。
Unityで送信したトピックをROSで受信する、というものを書きたいです。

### 発生している問題・エラーメッセージ

始めはUnityから送信した動作させるトピックでwhileに入り、
それをUnityから送信した停止させるトピックでwhileを抜け出したいのですが、
無限ループ中だとトピックが受信しませんでした。

UnityとROSの間で起きている問題なのかと思ったのですが、rostopic echoで確認するとトピックは受信しているようでして、無限ループが問題だとしてもROSで実行している/tfや/cmd_velは常に受信できているので無限ループが原因なのかもはっきりせず、ROSで実行しているトピックは無限ループ中でも受信できるのに対し、Unityで送信すると受信しないのはどうしてなのかわかりませんでした。
無限ループ中のUnityからのトピックを受信するためにはどうしたらいいのか、そしてなぜUnityからのトピックは無限ループ中受信しないのか教えてください。

```ここに言語を入力
#ROS#

class Test():
  def __init__(self):
    rospy.Subscriber('/cnd_vel', Twist, self.callback)
    rospy.Subscriber('/chatter', String, self.sub)
    rospy.spin()

  def callback(self, msg):
    print(msg.linear.x)　　#常に出力される

  def sub(self, s):
    if s.data == "Move":
      loop = "move"
    elif s.data == "Stop":
      loop = "stop"

    print(loop)　　　#無限ループ中は出力されない
    self.test(loop)

  def test(self, t)
    while t == "move":
      print("move")
      rospy.sleep(1)

      if t == "stop":　　　#無限ループ中は受信されないため、入らない
        break

    print("stop!")

if __name__ == '__main__':
  Test()
```

```ここに言語を入力
#Unity#

using System.Collections;
using System.Collections.Generic;
using UnityEngine;
using Unity.Robotics.ROSTCPConnector;
using RosMessageTypes.Std;

public class publisher : MonoBehaviour
{
　　ROSConnection ros;

　　void Start()
　　{
　　　　ros = ROSConnection.GetOrCreateInstance();
　　　　ros.RegisterPublisher<StringMsg>(“chatter”);
　　}
　　public void StartButtonClick()
　　{
　　　　StringMsg s = new StringMsg(“Move”);
　　　　ros.Publish(“chatter”, s);
　　}
　　public void StopButtonClick()
　　{
　　　　StringMsg s = new StringMsg(“Stop”);
　　　　ros.Publish(“chatter”, s);
　　}
}
```

### 試したこと

無限ループを数回実行して抜け出してみましたら、抜け出した後にトピックを受信してました。
Tkinterをそのまま使用してトピックでTkinterを操作してみましたが、無限ループ中は受信しませんでした。
しかし、Tkinterのボタンを押すと停止して抜け出せて、抜け出した後にトピックを受信してました。

### 補足情報（FW/ツールのバージョンなど）

ubuntu18
python2


unityからトピックを送信してrosでwhile文の無限ループ中にトピック受信して抜け出せなくて困っています

### 実現したいこと

ここに実現したいことを箇条書きで書いてください。
- [ ] ▲▲機能を動作するようにする

### 前提

ここに質問の内容を詳しく書いてください。
（例）
TypeScriptで●●なシステムを作っています。
■■な機能を実装中に以下のエラーメッセージが発生しました。
Pythonでwhile文を使い当選確率１％の１００連ガチャを当たりが出るまで
### 発生している問題・エラーメッセージ

```
エラーメッセージ
```

### 該当のソースコード

```ここに言語名を入力
ソースコード
```

### 試したこと

ここに問題に対して試したことを記載してください。

### 補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください


Pythonでwhile文を使い当たりが出るまでガチャを引く

### 実現したいこと

ここに実現したいことを箇条書きで書いてください。
- [ ] ▲▲機能を動作するようにする

### 前提

ここに質問の内容を詳しく書いてください。
（例）
TypeScriptで●●なシステムを作っています。
■■な機能を実装中に以下のエラーメッセージが発生しました。
Pythonでwhile文を使い当選確率１％のガチャを当たりが出るまで引き続けるプログラムを書きなさいという課題が出ました。１～１００のランダムな数字を使って、１がでたら当たり、それ意外なら外れと表示させることはできたのですが当たるまで繰り返す方法がわかりません
### 発生している問題・エラーメッセージ

```
エラーメッセージ
```

### 該当のソースコード

```ここに言語名を入力
ソースコード
```

### 試したこと

ここに問題に対して試したことを記載してください。

### 補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。


Pythonでwhile文を使ってガチャを作る方法

C言語の自作コンパイラを作成したいと思っています。以下のlex, yacc, cのコードを用いてwhile文を表現したいのですが、うまくいきません。**test.c の T_WHILE**の部分がおかしいのかなと思うのですが、どのように変更したらよいか教えていただきたいです。
```test.c
#include <stdio.h>
struct stab {
  int val;
  char name[20]; } stab[100];
int stabuse = 0;
struct node {
  int type, left, right; } ntab[400];
int ntabuse = 1;
#define T_STLIST 1
#define T_ASSIGN 2
#define T_READ 3
#define T_PRINT 4
#define T_ADD 5
#define T_SUB 6
#define T_MUL 7
#define T_DIV 8
#define T_REM 9
#define T_NUM 10
#define T_VAR 11
#define T_WHILE 12
#define T_IF 13
#define T_LT 14
#define T_GT 15
int lookup(char*);
int node(int, int, int);
void dotree(int);
extern char *yytext;
#include "y.tab.c"
#include "lex.yy.c"

int main() {
  yyparse();
  return 0;
}
int lookup(char *s) {
  int i;
  for(i = 0; i < stabuse; ++i)
    if(strcmp(stab[i].name, s) == 0) return i;
  if(stabuse >= 99) {
    printf("table overflow.\n"); exit(1); }
  strcpy(stab[stabuse].name, s); return stabuse++;
}
int node(int t, int l, int r) {
  int i = ntabuse++;
  ntab[i].type = t;
  ntab[i].left = l;
  ntab[i].right = r;
  return i;
}
void dotree(int i) {
  int stk;
  printf(" .section .rodata\n");
  printf(".Lprompt: .string\"> \"\n"); /* プロンプト */
  printf(".Lread: .string\"%%ld\"\n"); /* 読み取り用書式 */
  printf(".Lprint: .string\"%%ld\\n\"\n"); /* 書き出し用書式 */
  printf(" .text\n");
  printf(".global main\n");
  printf("main:\n");
  printf(" pushq %%rbp\n");
  printf(" movq %%rsp,%%rbp\n");
  stk = (8*stabuse + 15) / 16; /* 変数の個数*8 で 16 の倍数へ切り上げ */
  stk *= 16;
  printf(" subq $%d,%%rsp\n", stk);
  emittree(i);
  printf(" leave\n");
  printf(" ret\n");
}

void emittree(int i)
{
  static int labelno = 1;
  int l;
  switch(ntab[i].type) {
  case T_STLIST: if(ntab[i].left) emittree(ntab[i].left);
    emittree(ntab[i].right);
    break;
  case T_READ: printf(" movq $.Lprompt,%%rdi\n");
    printf(" movq $0,%%rax\n"); /* 浮動小数点レジスタを使わない */
    printf(" call printf\n");
    printf(" leaq %d(%%rbp),%%rsi\n", -(ntab[i].left+1)*8);
    printf(" movq $.Lread,%%rdi\n");
    printf(" movq $0,%%rax\n"); /* 浮動小数点レジスタを使わない */
    printf(" call scanf\n");
    break;
  case T_PRINT: emittree(ntab[i].left);
    printf(" popq %%rsi\n");
    printf(" movq $.Lprint,%%rdi\n");
    printf(" movq $0,%%rax\n"); /* 浮動小数点レジスタを使わない */
    printf(" call printf\n");
    break;
  case T_NUM: printf(" pushq $%d\n", ntab[i].left);
    break;
  case T_VAR: printf(" pushq %d(%%rbp)\n", -(ntab[i].left+1)*8);
    break;
  default: printf("NotImplemented: %d\n", ntab[i].type);
    break;
  case T_ASSIGN: emittree(ntab[i].right);
    printf(" popq %d(%%rbp)\n", -(ntab[i].left+1)*8);
    break;
  case T_ADD: emittree(ntab[i].left);
    emittree(ntab[i].right);
    printf(" popq %%rdx\n");
    printf(" popq %%rax\n");
    printf(" addq %%rdx,%%rax\n");
    printf(" pushq %%rax\n");
    break;
  case T_SUB: emittree(ntab[i].left);
    emittree(ntab[i].right);
    printf(" popq %%rdx\n");
    printf(" popq %%rax\n");
    printf(" subq %%rdx, %%rax\n");
    printf(" pushq %%rax\n");
    break;
  case T_MUL: emittree(ntab[i].left);
    emittree(ntab[i].right);
    printf(" popq %%rdx\n");
    printf(" popq %%rax\n");
    printf(" imulq %%rdx\n");
    printf(" pushq %%rax\n");
    break;
  case T_DIV: emittree(ntab[i].left);
    emittree(ntab[i].right);
    printf(" popq %%rsp\n");
    printf(" popq %%rax\n");
    printf(" movq $0,%%rdx\n");
    printf(" idivq %%rsp\n");
    printf(" pushq %%rax\n");
    break;
  case T_REM: emittree(ntab[i].left);
    emittree(ntab[i].right);
    printf(" popq %%rsp\n");
    printf(" popq %%rax\n");
    printf(" movq $0,%%rdx\n");
    printf(" idivq %%rsp\n");
    printf(" pushq %%rdx\n");
  case T_LT: emittree(ntab[i].left);
    emittree(ntab[i].right);
    printf(" popq %%rcx\n");
    printf(" popq %%rax\n");
    printf(" cmp %%rcx,%%rax\n");
    printf(" jge ");
    break;
  case T_GT: emittree(ntab[i].left);
    emittree(ntab[i].right);
    printf(" popq %%rcx\n");
    printf(" popq %%rax\n");
    printf(" cmp %%rcx,%%rax\n");
    printf(" jle ");
    break;
  case T_IF: l = labelno++;
    emittree(ntab[i].left);
    printf(".L%d\n", l);
    emittree(ntab[i].right);
    printf(".L%d:\n", l);
    break;
  case T_WHILE:  　　　　 // この部分
　 l = labelno++;　　// ラベル番号更新
    printf(".L%d:\n", l);   // ラベル l 生成
    emittree(ntab[i].left);　// 分岐条件判定
    printf(".L%d\n", l+1);　// 不成立なら l+1に移動
    emittree(ntab[i].right);　// 成立していたら、中身を実行
    printf(".L%d\n", l);　// 前のラベル l に戻る
    printf(".L%d:\n", l+1);　// ラベル l+1 生成
    break;
  }
}
```
```test.yacc
%token NUM;
%token IDENT;
%token READ;
%token PRINT;
%token WHILE;
%token IF;
%left '+' '-';
%left '*' '/';
%%
prog : IDENT '{' stlist'}' { dotree($3); return 0; }
;
stlist : { $$ = 0; }
| stlist stat { $$ = node(T_STLIST, $1, $2); }
;
stat : var '=' expr ';' { $$ = node(T_ASSIGN, $1, $3); }
| READ var ';' { $$ = node(T_READ, $2, 0); }
| PRINT expr ';' { $$ = node(T_PRINT, $2, 0); }
| WHILE '(' cond ')' stat { $$ = node(T_WHILE, $3, $5); }
| IF '(' cond ')' stat { $$ = node(T_IF, $3, $5); }
| '{' stlist '}' { $$ = $2; }
;

cond : expr '<' expr { $$ = node(T_LT, $1, $3); }
| expr '>' expr { $$ = node(T_GT, $1, $3); }
;
expr : term { $$ = $1; }
| expr '+' term { $$ = node(T_ADD, $1, $3); }
| expr '-' term { $$ = node(T_SUB, $1, $3); }
;
term: prim{ $$ = $1; }
| term '*' prim { $$ = node(T_MUL, $1, $3); }
| term '/' prim { $$ = node(T_DIV, $1, $3); }
| term '%' prim { $$ = node(T_REM, $1, $3); }
;

prim : NUM { $$ = node(T_NUM, atoi(yytext), 0); }
| var { $$ = node(T_VAR, $1, 0); }
| '(' expr ')' { $$ = $2; }
;
var : IDENT { $$ = lookup(yytext); }
;
```
```test.lex
alpha [a-zA-Z]
digit [0-9]
white [\n\t ]
%%
while { return WHILE; }
if { return IF; }
read { return READ; }
print { return PRINT; }
{alpha}({alpha}|{digit})* { return IDENT; }
{digit}+ { return NUM; }
[-+()=;{}<>*/%] { return yytext[0]; }
{white} { ; }
```

C言語の自作コンパイル　while文の表現

```Text
whileはif, gotoで代用できる。
よって、あなたはif, gotoのみを用いるべきである。
```

whileはif, gotoで代用できる。よって、あなたはif, gotoのみを用いるべきである。

### 実現したいこと

while文での繰り返し処理を、1枚のcsvでまとめて、出力したいです。
下記コードでは、最後の結果しか出力されません。
ご教示お願いしたいです。

作成したcsvファイルで、下記のように出力したいです。

1回目
2回目
・
・
5回目


### 該当のソースコード

```ここに言語名を入力
S = 0
B = "回目"
while S<30:
    S = S + 1
    C = str(S) + B
    print(C)

#csvの出力
f = open("ファイル名.csv", "w")
writecsv = csv.writer(f, lineterminator='\n')
writecsv.writerow(C)
f.close()
```

### 試したこと

while文の中にwritecsvを入れたら、エラー発生しています。

### 補足情報（FW/ツールのバージョンなど）
Python 3.11.0


【Python】while文の繰り返し処理の結果をまとめて1枚のcsvファイルで出力したい

### 前提
プログラムは量子物性
whileで同じ動作を繰り返したい。

### 実現したいこと

whileの中にある動作を一定回数繰り返してほしいが、コンパイルしても同じ値しか出てこない
### 発生している問題・エラーメッセージ

回数にかかわらず同じ値しか出てこない。

### 該当のソースコード

```ｃ
//グラフェンナノリボンエネルギーバンド　相互作用
#include <stdio.h>
#include <math.h>
#include <complex.h>
#include <stdlib.h>

#define M  6	// 存在するサイト数
#define N  12		//サイト数×２　A+Bの個数
#define L  24  	//対角化に必要　N*２
#define P  400		//リボン長　全電子数　P*L
#define EPS 0.000001  	//収束範囲
#define Q 10        //回数

double b_pi,N_c;
double complex t_1, t_2, t_3, t_4, t_5, t_6, t_7, t_8, T_1, T_2, T_3, T_4,T_int;		//,t_方向
double e[P+1],k_x[P+1];
double B_1[N][N], B_2[N][N],D[N][N];
double complex C[N][N];
double a[L][L],A[N];
double Energy[P+1][N];
double DOS[P+1][P+1][N];
double ddd[P+1][P+1],dd[P+1];
double fermi[P+1][N],fermisum_1[P+1],f[P+1];
double complex u[P+1][N][N];		//固有ベクトル
double u_2[P+1][N],u_3[P+1],u_4[P+1],u_5[P+1][N][N],u_6[P+1][N],u_7[N],u_8[N];
double gap[P+1],meow[P+1][N],min_gap;       //Eg,フェルミエネルギーバンド、最小Eg

int main(int argc, char** argv){
    b_pi=M_PI;
	N_c=(P+1)*N;
    double k_xx = -b_pi; 
    double ene=-3.0;
	double delta_E=2.75*pow(10,-2);		//エネルギー分解能
	double kt=0.00001*8.6171*pow(10,-5);		//-273℃＊ボルツマン定数(eV/K)
	double U=4.32;		//相互作用エネルギーの平均(eV)
	double t=2.75;		//飛び移り積分(eV)
	double n_up[N];		//代入する値
	double ribbon_length;		//リボン長

    FILE* fp;
    FILE* fp_1;
    FILE* fp_2;
	FILE* fp_3;

	if ((fp = fopen("zigintene.txt", "w")) == NULL)
	{
		printf("Cannot open the file\n");
		exit(1);
	}
	if ((fp_1 = fopen("zigintdos.txt", "w")) == NULL)
	{
		printf("Cannot open the file\n");
		exit(1);
	}
	if ((fp_2 = fopen("zigintnup_1.txt", "w")) == NULL)
	{
		printf("Cannot open the file\n");
		exit(1);
	}
	fp_3 = fopen("zigintnup.txt", "r");   
    if(fp_3 == NULL) {
        printf("ファイルを開くことが出来ませんでした．¥n");
    }
	for(int i=0;i<N;i++){
		fscanf(fp_3, "%lf", &(n_up[i]));  
	}
	double omori=0; //0 or 1
    
	//刻み幅
    for(int i=0;i<=P;i++){
        e[i]=ene+6.0/P*i;
        k_x[i]=k_xx+2.0*b_pi/P*i;
    }

	//リボン幅
	ribbon_length=0.246*pow(3,0.5)/2*(N-2)/2;

	//初期条件
	/*for(int i=0;i<N;i++){
		if(i%2==1){
			n_up[i]=0.4;
		}
		else if(i%4==3){
			n_up[i]=-0.3;
		}
		else{
			n_up[i]=0;
		}
	}*/
	
    int z=0;
    while(z<4){
	//行列計算
    for(int r=0;r<=P;r++){
        t_4 = 1;
	    t_5 = cexp(-1 * I * -1/ 2 * k_x[r]);			
        t_6 = cexp(-1 * I *1/2 * k_x[r]);		
        t_1 =1;
    	t_2 = cexp(I * -1 / 2 * k_x[r]);
	    t_3 = cexp(I * 1 / 2 * k_x[r]);
    	t_7 = -1 * t_4;
	    t_8 = -1 * t_1;
    	T_1 = t_5 + t_6;
     	T_2 = t_2 + t_3;
	    T_3 = -1 * T_1;
	    T_4 = -1 * T_2;

        for (int k = 0;k < N;k++) {
			for (int h = 0;h < N;h++) {
				T_int=-U/t*n_up[k];				//強磁性
				for (int i = 0;i < M - 2;i++) {
					if (h < M) {
						if (k < M) {
							C[h][k] = 0;
							C[k][k]=T_int;
						}
						else {
							C[0][M] = T_1;
							C[M - 1][N - 1] = T_1;
							C[M - 1][N - 2] = t_4;
							C[i + 1][M + i] = t_4;
							C[i + 1][M + i + 1] = T_1;
						}
					}
					else {
						if (k < M) {
							C[M][0] = T_2;
							C[M][1] = t_1;
							C[N - 1][M - 1] = T_2;
							C[M + i + 1][i + 1] = T_2;
							C[M + i + 1][i + 2] = t_1;

						}
						else {
							C[h][k] = 0;
							C[k][k]=T_int;
						}
					}
				}
			}
		}

		int k, h;
		for (k = 0;k < N;k++) {
			for (h = 0;h < N;h++) {
				D[h][k] = creal(C[h][k]);
				B_1[h][k] = -1 * cimag(C[h][k]);
				B_2[h][k] = cimag(C[h][k]);
			}
		}

		for (k = 0;k < L;k++) {
			for (h = 0;h < L;h++) {
				if (h < N) {
					if (k < N) {
						a[h][k] = D[h % N][k % N];
					}
					else {
						a[h][k] = B_1[h % N][k % N];
					}
				}
				else {
					if (k < N) {
						a[h][k] = B_2[h % N][k % N];
					}
					else {
						a[h][k] = D[h % N][k % N];
					}
				}
			}
		}

    	for(int i=0;i<N;i++){
			for(int j=0;j<N;j++){
				u[r][i][j]=0.0;
			}
		}
		for(int i=0;i<N;i++){
			u[r][i][i]=1.0;
		}

        while (1) {		

			double alpha, beta, gamma;
			double s, c, w;
			double wa, wb, wc;
			double max;
			int i, j, p, q, x, y;

			//最大要素の行と列を検索
			max = 0.0;
			for (i = 0;i < L - 1;i++) {
				for (j = i + 1;j < L;j++)
					if (fabs(a[i][j]) > max) {
						p = i;
						q = j;
						max = fabs(a[i][j]);
					}
			}
			//収束したら解答打ち出し
			if (max < EPS) break;

			//sin, cos　計算
			wa = a[p][p];
			wb = a[p][q];
			wc = a[q][q];
			alpha = -wb;
			beta = 0.5 * (wa - wc);
			gamma = fabs(beta) / sqrt(alpha * alpha + beta * beta);
			s = sqrt(0.5 * (1.0 - gamma));
			if (alpha * beta < 0) s = -s;
			c = sqrt(1.0 - s * s);
			//直行変換
			for (j = 0;j < L;j++) {
				w = a[p][j] * c - a[q][j] * s;
				a[q][j] = a[p][j] * s + a[q][j] * c;
				a[p][j] = w;
			}

			for (j = 0;j < L;j++) {
				a[j][p] = a[p][j];
				a[j][q] = a[q][j];
			}

			w = 2.0 * wb * s * c;
			a[p][p] = wa * c * c + wc * s * s - w;
			a[q][q] = wa * s * s + wc * c * c + w;
			a[p][q] = 0;
			a[q][p] = 0;

			for (int i = 0;i < N;i++)
			{
				A[i] = a[i][i];
			}

            double tmp;
		 	for (int i=0; i<N; i++) {
 	   			for (int j=i+1; j<N ;j++) {
    			  if (A[i] > A[j]) {
   	    	 		tmp =  A[i];
   	    	 		A[i] = A[j];
   	    	 		A[j] = tmp;
  	    			}
	   			}			
	 		}
			for(i=0;i<N;i++){
				w=u[r][i][p]*c-u[r][i][q]*s;
				u[r][i][q]=u[r][i][p]*s+u[r][i][q]*c;
				u[r][i][p]=w;
			}
        }

        for(int i=0;i<N;i++){
            Energy[r][i]=A[i];
        }
	}

	/*for (int i=0; i<=P; i++) {
		gap[i]=Energy[i][M]-Energy[i][M-1];
		//printf( "%2f %2f %2f\n",Energy[i][M-1],Energy[i][M],gap[i]);
	}
	for(int i = 0; i <=P; i++){
    	if(min_gap > gap[i])  
    	{
      		min_gap = gap[i];  
    	}
	}*/

	//フェルミエネルギー
	for(int i=0;i<=P;i++){
		for(int j=0;j<N;j++){
			meow[i][j]=(Energy[i][j]+Energy[i][N-1-j])/2;
		}
	}

	//状態密度
	for(int i=0;i<=P;i++){
		for(int j=0;j<=P;j++){
			for(int k=0;k<N;k++){
				DOS[i][j][k]=1/pow(b_pi*2,0.5)/delta_E/b_pi*exp(-1*pow(e[i]-Energy[j][k],2)/2/delta_E/delta_E);
			}
		}
	}
	for(int i=0;i<=P;i++){
		for(int j=0;j<=P;j++){
			for(int k=0;k<N;k++){
				ddd[i][j]+=DOS[i][j][k];
			}
		}
	}
	for(int i=0;i<=P;i++){
		for(int j=0;j<=P;j++){
			dd[i]+=ddd[i][j];
		}
	}

	//波動関数の二乗＝確率
	for(int l=0;l<=P;l++){
		for(int i=0;i<N;i++){
			for(int j=0;j<N;j++){
				u_2[l][i]+=pow(cabs(u[l][i][j]),2);
			}
		}
	}
	for(int l=0;l<=P;l++){
		for(int i=0;i<N;i++){
			u_3[l]+=u_2[l][i];
		}
	}
	for(int l=0;l<=P;l++){
		u_4[l]=u_3[l]/N;
		//printf( "%7.4lf\n", u_4[l]);	//１になって然るべき
	}

	//個数をフェルミ関数から計算
	for(int l=0;l<=P;l++){
		for(int i=0;i<N;i++){
			fermi[l][i]=1/(1+exp((Energy[l][i]-meow[l][i%M])/kt));
		}
		//printf( "%2f\n",fermi[l][0]);	
	}
	for(int l=0;l<=P;l++){
		for(int i=0;i<M;i++){
			f[l]+=fermi[l][i];			//Mと等しくなる
		}	
		//printf( "%2f\n",f[l]);
	}
	//printf( "%2f\n",f[0]);		//粒子数（サイト数）
	for(int i=0;i<N;i++){
		for(int l=0;l<=P;l++){
			for(int j=0;j<N;j++){
				u_5[l][i][j]=pow(cabs(u[l][i][j]),2)*fermi[l][j];
				u_6[l][i]+=u_5[l][i][j];
			}
			u_7[i]+=u_6[l][i];
		}
	}	
	for(int i=0;i<N;i++){
		if(i%4==omori+1){
			u_8[i]=n_up[i];				
			}
		else if(i%4==omori+2){
			u_8[i]=n_up[i];
		}
		else{
			u_8[i]=u_7[i];
		}
		u_8[i]=u_7[i];
		n_up[i]=u_8[i]/(P+1);
	}

    z++;
    }
	for(int i=0;i<=P;i++){
        for(int j=0;j<N;j++){
        	fprintf(fp, "%2f %lf %lf %lf\n",k_x[i],Energy[i][j]-meow[i][j],Energy[i][M-1]-meow[i][j],Energy[i][M]-meow[i][j]);
			//fprintf(fp, "%2f %lf %lf\n",k_x[i],Energy[i][M-1],Energy[i][M]);		//最低バンド
        }
    }
	for(int i=0;i<=P;i++){
		fprintf(fp_1, "%2f %lf\n",e[i],dd[i]);		//状態密度
	}
	for(int i=0;i<N;i++){	
		fprintf(fp_2,"%2f\n",u_8[i]/(P+1));	
	}

   	printf("リボン長：%2f\n", ribbon_length);
    printf("最小ギャップ：%f\n",min_gap);

	fclose(fp);
	fclose(fp_1);
	fclose(fp_2);
	fclose(fp_3);
	return 0;
}
```


