どのようなコードを書けば良いのでしょうか?よろしくお願いします。
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define M 4
#define N 100
#define delta 1.0
int main(){
int i,j,k,n;
int b[M][M] = {0};
double u[2*M+1][N];
for(i=0;i<M;i++){
b[i][i]=1;
}
for(i=0;i<M-1;i++){
b[i][i+1]=1;
}
for(n=0;n<N;n++){
u[0][n]=0;
u[2*M][n]=0;
}
for(k=1;k<2*M;k++){
u[k][0] = 1.0/(1+u[k-1][0]);
}
for(n=0;n<N-1;n++){
for(k=1;k<2*M;k++){
u[k][n+1] = (u[k][n]*(1+delta*u[k+1][n])) / (1+delta*u[k-1][n+1]);
}
}
for(n=0;n<N;n++){
for(k=1;k<2*M;k++){
u[k][n] = sqrt(u[k][n]);
}
}
for(n=0;n<N;n++){
for(k=0;k<M;k++){
printf("%.9f ",u[2*k+1][n]);
}
printf("\n");
}
return 0;
}