http://w.atwiki.jp/ferret_modoki/ 工房ferret ja 2009-06-04T12:35:42+09:00 1244086542 https://w.atwiki.jp/ferret_modoki/pages/81.html fb 2009-06-04T12:35:42+09:00 1244086542 https://w.atwiki.jp/ferret_modoki/pages/80.html i 2009-06-04T12:33:07+09:00 1244086387 https://w.atwiki.jp/ferret_modoki/pages/79.html t 2009-06-04T12:29:38+09:00 1244086178 https://w.atwiki.jp/ferret_modoki/pages/78.html i 2009-06-04T12:28:03+09:00 1244086083 https://w.atwiki.jp/ferret_modoki/pages/77.html #include <stdio.h> #include <stdlib.h> #include <math.h> main() { int i,j,k,a,b; double m,n,h,h2,o; static double L[700][701]; printf("分割数n="); scanf("%lf",&n); h = 1 / n; h2 = h*h-2; a = n+1; b = n+2; printf("n=%lf, h=%lf, h^2-2=%lf, a=%d, b=%d\n\n",n,h,h2,a,b); /*行列入力*/ L[0][0] = 1; L[n][n] = 1; for(i=1; i<a-1; i++){ L[i][i] = h2; if(i != 1){ L[i][i-1] = 1; } if(i != a-2){ L[i][i+1] = 1; } } for(i=1; i<n; i++){ L[i][a] = -1*h*h*h*i; } /*掃き出し実行*/ for(k=0; k<a; k++){ for(i=0; i<a; i++){ if(i != k){ m = L[i][k]; for(j=0; j<b; j++){ L[i][j] = L[i][j] - ((m / L[k][k]) * L[k][j]); } } } m = L[k][k]; for(j=k; j<b; j++){ L[k][j] = L[k][j] / m; } } /* 出力 */ for(i=0;i<a;i++){ m = (i*h); if( m==.2 || m==.4 || m==.6 || m==.8 || (.599999<m && m<.600001)) printf("u%lf = %lf \n",m,L[i][a]); } return 0; } /* 10 u0.200000 = 2009-05-29T14:55:45+09:00 1243576545 https://w.atwiki.jp/ferret_modoki/pages/76.html -[[掃き出しプログラム]] 2009-05-29T14:55:36+09:00 1243576536 https://w.atwiki.jp/ferret_modoki/pages/75.html /* ---------------------------------------- Image Processing pg1.c Image Ratio Processing / NDVI Gray Image Output (c)Ichinoseki National College Advanced Course ------------------------------------------ */ #include <stdio.h> #include <stdlib.h> #include <math.h> main(argc,argv) int argc; char *argv[]; { static int i,j,m; static double a1,a2,a3; static int pixel,line; FILE *fp1,*fp2,*fp3; /* Command parameter check */ if (argc!=6){ printf("Usage : command input_file1 input_file2 pixel line output_file\n"); return 1; } /* Image size parameter load */ sscanf(argv[3],"%d",&pixel); sscanf(argv[4],"%d",&line); /* Input image file1 check */ if (( fp1=fopen(argv[1],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Input image file1 check */ if (( fp2=fopen(argv[2],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Output image file check */ if (( fp3=fopen(argv[5],"w") )==NULL){ 2009-05-28T12:00:11+09:00 1243479611 https://w.atwiki.jp/ferret_modoki/pages/74.html /* ---------------------------------------- Image Processing pg1.c Image Ratio Processing / NDVI Gray Image Output (c)Ichinoseki National College Advanced Course ------------------------------------------ */ #include <stdio.h> #include <stdlib.h> #include <math.h> main(argc,argv) int argc; char *argv[]; { static int i,j,m; static double a1,a2,a3; static int pixel,line; FILE *fp1,*fp2,*fp3; /* Command parameter check */ if (argc!=6){ printf("Usage : command input_file1 input_file2 pixel line output_file\n"); return 1; } /* Image size parameter load */ sscanf(argv[3],"%d",&pixel); sscanf(argv[4],"%d",&line); /* Input image file1 check */ if (( fp1=fopen(argv[1],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Input image file1 check */ if (( fp2=fopen(argv[2],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Output image file check */ if (( fp3=fopen(argv[5],"w") )==NULL){ 2009-05-28T11:59:41+09:00 1243479581 https://w.atwiki.jp/ferret_modoki/pages/73.html /* ---------------------------------------- Image Processing pg1.c Image Ratio Processing / NDVI Gray Image Output (c)Ichinoseki National College Advanced Course ------------------------------------------ */ #include <stdio.h> #include <stdlib.h> #include <math.h> main(argc,argv) int argc; char *argv[]; { static int i,j,m; static double a1,a2,a3; static int pixel,line; FILE *fp1,*fp2,*fp3; /* Command parameter check */ if (argc!=6){ printf("Usage : command input_file1 input_file2 pixel line output_file\n"); return 1; } /* Image size parameter load */ sscanf(argv[3],"%d",&pixel); sscanf(argv[4],"%d",&line); /* Input image file1 check */ if (( fp1=fopen(argv[1],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Input image file1 check */ if (( fp2=fopen(argv[2],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Output image file check */ if (( fp3=fopen(argv[5],"w") )==NULL){ 2009-05-28T11:59:07+09:00 1243479547 https://w.atwiki.jp/ferret_modoki/pages/72.html [[pg1]] [[pg2]] [[pg3]] 2009-05-28T11:58:41+09:00 1243479521