コメント一覧
まえへ つぎへ59 投稿者：名無しさん 2019/11/05 12:15:46
ここで一句。
60 投稿者：名無しさん 2019/11/05 12:15:55
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define N 4
int main(){
int i,j,k,p;
float A[N][N+1],x[N],q;
FILE *fp;
if( (fp=fopen("20191105-1.dat","r"))==NULL) {
printf("cannot open file.\n");
exit(0);
}
for(i=0;i<N;i++)
for(j=0;j<N+1;j++)
fscanf(fp,"%f",&A[i][j]);
fclose(fp);
for(i=0;i<N;i++){
for(j=0;j<N+1;j++) printf("%6.2f",A[i][j]);
printf("\n");
}
printf("\n");
/* forward elimination */
for(k=0;k<N-1;k++){
for(j=k+1;j<N;j++){
if(fabs(A[k][k])<1.0e-6){
printf("A[k][k] is zero.\n");
exit(0);
}
q=A[j][k]/A[k][k];
for(p=k;p<N+1;p++) A[j][p]-=q*A[k][p];
}
}
/*backward substitution*/
for(i=N-1;i>=0;i--){
x[i]=A[i][N];
for(j=i+1;j<N;j++) x[i]-=A[i][j]*x[j];
x[i]/=A[i][i];
}
/*output result*/
for(i=0;i<N;i++){
for(j=0;j<N+1;j++) printf("%6.2f",A[i][j]);
printf("\t%6.2f\n",x[i]);
}
}
61 投稿者：名無しさん 2019/11/05 12:16:05
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define N 3
int main(void){
int i,j,k;
float A[N][N],L[N][N],U[N][N],x[N],b[N],y[N];
FILE *fp;
if( (fp=fopen("20191105-2.dat","r"))=NULL){
printf("cannnot open file.\n");
exit(0);
}
for(i=0;i<N;i++){
for(j<0;j<K;j++) fscanf(fp, "%f",&A[i][j]);
fscanf(fp,"%f",&b[i]);
}
fclose(fp);
for(i=0;i<N;i++){
for(j=0;j<N+1;j++) printf("%6.2f",A[i][j]);
printf("\n");
}
printf("\n");
for(i=0;i<N;i++){
for(j=0;j<N;j++){
L[i][j]=0.0;
U[i][j]=0.0;
}
x[i]=0.0;
}
/*A=LU (U=LT) -> L,U*/
for(j=0;j<N;j++){
for(i=j;i<N;i++){
if(i==j){
L[i][j]=A[i][j];
for(k=0;k<i;k++)
L[i][j]-=L[i][k]*L[i][k];
L[i][j]=sqrt(L[i][j]);
} else {
L[i][j]=A[i][j];
for(k=0;k<j;k++)
L[i][j]-=L[j][k]*L[i][k];
L[i][j]/=L[j][j];
}
U[j][i]=L[i][j];
}
}
for(i=0;i<N;i++){
for(j=0;j<N;j++) printf("%6.2f",L[i][j]);
printf("\t");
for(j=0;j<N;j++) printf("%6.2f",U[i][j]);
printf("\n");
}
printf("\n");
/*Ly=b -> y*/
for(i=0;i<N;i++){
y[i]=b[i];
for(j=0;j<i;j++) y[i]-=L[i][j]*y[j];
f[i]/=L[i][i];
}
/*Ux->x*/
for(i=N-1;i>=0;i--){
x[i]=y[i];
for(j=i+1;j<N;j++) x[i]-=U[i][j]*x[j];
x[i]/=U[i][i];
}
/*result*/
for(i=0;i<N;i++)printf("%6.2f\n",x[i]);
}
62 投稿者：名無しさん 2019/11/05 12:20:04
ここで一句。
63 投稿者：名無しさん 2019/11/05 12:20:20
//20191105_2 MaedaSoichiro
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define N 3
int main(void){
int i,j,k,p;
float A[N][N],L[N][N],U[N][N],x[N],b[N],y[N];
FILE *fp;
if((fp=fopen("191105_2.dat","r"))==NULL){
printf("cannot open file.\n");
exit(0);
}
for(i=0;i<N;i++){
for(j=0;j<N+1;j++)
fscanf(fp,"%f",&A[i][j]);
fscanf(fp,"%f",&b[i]);
}
fclose(fp);
for(i=0;i<N;i++){
for(j=0;j<N+1;j++)printf("%6.2f",A[i][j]);
printf("\n");
}
printf("\n");
/**/
for(i=0;i<N;i++){
for(j=0;j<N;j++){
L[i][j]=0.0;
U[i][j]=0.0;
}
x[i]=0.0;
}
/*A=LU (U=LT) -> L,U */
for(j=0;j<N;j++){
for(i=j;i<N;i++){
if(i==j){
L[i][j]=A[i][j];
for(k=0;k<i;k++)
L[i][j]-=L[i][k]*L[i][k];
L[i][j]=sqrt(L[i][j]);
}else{
L[i][j]=A[i][j];
for(k=0;k<j;k++)
L[i][j]-=L[j][k]*L[i][k];
L[i][j]/=L[j][j];
}
U[j][i]=L[i][j];
}
}
for(i=0;i<N;i++){
for(j=0;j<N;j++)printf("%6.2f",L[i][j]);
printf("\t");
for(j=0;j<N;j++)printf("%6.2f",U[i][j]);
printf("\n");
}
printf("\n");
/*ly=b -> y */
for(i=0;i<N;i++){
y[i]=b[i];
for(j=0;j<i;j++)y[i]-=L[i][j]*y[j];
y[i]/=L[i][i];
}
/*Ux=y -> x */
for(i=N-1;i>=0;i--){
x[i]=y[i];
for(j=i+1;j<N;j++)x[i]-=U[i][j]*x[j];
x[i]/=U[i][i];
}
/*result*/
for(i=0;i<N;i++)printf("%6.2f\n",x[i]);
}
/*
4.00 -2.00 0.00 -2.00
-2.00 2.00 -1.00 0.00
0.00 -1.00 2.00 0.00
*/
64 投稿者：名無しさん 2019/12/05 12:57:14
65 投稿者：名無しさん 2019/12/05 12:58:02
66 投稿者：名無しさん 2019/12/09 09:54:25
#define _WINSOCK_DEPRECATED_NO_WARNINGS
#include <WinSock2.h>
#include <stdio.h>
#include<stdlib.h>
#include<string.h>
#include<conio.h>
int main(){
SOCKET listen_s;
SOCKET s;
WSADATA wsaData;
SOCKADDR_IN saddr, from;
char szBuf[1024];
int fromlen;
int nRcv;
u_short uport;
printf("ポート番号：");
gets_s(szBuf);
uport = atoi(szBuf);
if (WSAStartup(MAKEWORD(2, 2), &wsaData)){
printf("WSAStartup関数失敗\n");
WSACleanup();
return -1;
}
listen_s = socket(AF_INET, SOCK_STREAM, 0);
if (listen_s < 0){
printf("ソケットオープン失敗\n");
WSACleanup();
return -1;
}
listen_s = socket(AF_INET, SOCK_STREAM, 0);
if (listen_s < 0){
printf("ソケットオープンエラー\n");
WSACleanup();
return -2;
}
else{
printf("ソケットオープン完了");
}
memset(&saddr, 0, sizeof(SOCKADDR_IN));
saddr.sin_family = AF_INET;
saddr.sin_port = htons(uport);
saddr.sin_addr.s_addr = INADDR_ANY;
if (bind(listen_s, (struct sockaddr *)&saddr, sizeof(saddr)) == SOCKET_ERROR){
printf("バインディングエラー");
closesocket(listen_s);
WSACleanup();
return -3;
}
else{
printf("バインディング完了\n");
}
if (listen(listen_s, 0) == SOCKET_ERROR){
printf("Listen移行失敗\n");
closesocket(listen_s);
WSACleanup();
return -4;
}
else{
printf("Listen移行完了\nAccept状態へ遷移\n");
}
fromlen = (int)sizeof(from);
s = accept(listen_s, (SOCKADDR *)&from, &fromlen);
if (s == INVALID_SOCKET){
printf("接続失敗");
closesocket(listen_s);
WSACleanup();
return -5;
}
else{
printf("%sから接続\n", inet_ntoa(from.sin_addr));
printf("アクセプト完了\n送受信開始\n\n");
}
closesocket(listen_s);
while (1){
printf("\n受信待機...\n");
nRcv = recv(s, szBuf, sizeof(szBuf) - 1, 0);
if (nRcv == SOCKET_ERROR){
printf("受信エラー\n");
break;
}
szBuf[nRcv] = '\0';
if (strcmp(szBuf, "c_end") == 0){
printf("クライアント：接続切断\n");
break;
}
send(s, szBuf, (int)strlen(szBuf), 0);
}
shutdown(s, SD_BOTH);
closesocket(s);
WSACleanup();
printf("サーバー終了\n");
_getch();
return 0;
}
67 投稿者：名無しさん 2019/12/10 11:39:25
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define C 4
#define Cn C+1
#define Cnn 2*C+1
void yacobi(float A[Cn][Cn],float X[Cn],float B[Cn])
{
int i,j;
float tmp,err;
do {
err=0.0;
for(i=0;i<Cn;i++) {
tmp=B[i];
for(j=0;j<Cn;j++)
if(i!=j) tmp-=A[i][j]*X[j];
tmp/=A[i][i];
err+=fabs(tmp-X[i]);
X[i]=tmp;
}
} while(err>1.0e-6);
}
int main(void)
{
int i,j,N;
float sx[Cnn],sy[Cnn],tx,ty,t;
float A[Cn][Cn],X[Cn],B[Cn];
float x[100],y[100];
FILE *fp;
if((fp=fopen("191210_1.dat","r"))==NULL){
printf("cannot open file\n");
exit(0);
}
i=0;
while( fscanf(fp,"%f %f",&x[i],&y[i])!=EOF )
i++;
N=i;
for(i=0;i<Cnn;i++) sx[i]=sy[i]=0.0;
for(i=0;i<N;i++) {
for(j=0;j<Cnn;j++) {
if(j==0) tx=1.0; else tx*=x[i];
sx[j]+=tx;
}
for(j=0;j<Cn;j++) {
if(j==0) ty=y[i]; else ty*=x[i];
sy[j]+=ty;
}
}
for(i=0;i<Cn;i++) {
for(j=0;j<Cn;j++) A[i][j]=sx[i+j];
X[i]=1.0;
B[i]=sy[i];
}
yacobi(A,X,B);
for(i=0;i<Cn;i++) {
if(i==0) printf("y=%.1f",X[0]);
else {
if(X[i]>0.0) printf("+");
printf("%.1f*x",X[i]);
}
if(i>1) printf("^%d",i);
}
printf("\n");
for(i=0;i<N;i++) {
t=X[2];
for(j=1;j>=0;j--) t=t*x[i]+X[j];
printf("%f %f(%f)\n",x[i],y[i],t);
}
}
/*
[e152137@ei1-pc14 ~]$ a.out
y=8.0-23.9*x+21.9*x^2-8.0*x^3+1.0*x^4
0.000000 8.000000(7.993975)
0.500000 0.563000(1.504619)
1.000000 -1.000000(5.982754)
1.500000 -0.438000(21.428379)
2.000000 0.000000(47.841499)
2.500000 -0.438000(85.222107)
3.000000 -1.000000(133.570206)
3.500000 0.563000(192.885788)
4.000000 8.000000(263.168884)
*/
68 投稿者：名無しさん 2019/12/10 14:50:02
まず現代日本においてなにが古典で、古典の価値とは何かについてである。文学において古典とは古い書物であるだけではない。高い評価を受け長きにわたり読み継がれることによって「世の中を変えてきたもの」。これは「表現・美意識・民族の感性に影響を与えてきたもの」ということであり、現代の日本民族・文化においての伝統的思考や意識、哲学的思想そのものが記録されているものと解釈することができる。
さて日本には神道の根幹とも言える独特な神に対する視点があり、今回の話題はソレ「八百万の神」についてである。古くは「古事記」上巻の「天の岩戸」が初見であり、他にも同様の総称として八十諸神・八十万神・八十万群神などが「日本書紀」「万葉集」などに見られる。コレが日本文化に及ぼしてきた影響はとても大きく、妖怪を例にとってもその影をいろいろなところに見出すことができる。
宮崎駿の「もののけ姫」では長く生きた猪が「猪神」は「伊勢物語」に登場する「付喪神」、「祟り神」は「日本書紀」、アシタカ自身は鬼と呼ばれるなど、様々なものを垣間見ることができる。