int main(int argc, char **argv){ static xyc *Z; static nde *N; FILE *fp; int i; initop(argc, argv); fp = argf(argc, argv); fp2msh(fp,&Z, &N); pltmsh(stdout,Z,N); }
main(int argc, char **argv){ int i, n = 512; static MX *A; static double *x, *b; initop(argc,argv); initmx(A,n+1,8); ary1(x, n+1); ary1(b, n+1); for ( i =1; i<=n; i++){ mx(A,i,i) = 2.0; b[i] = 1.0; } for (i=1;i<n;i++) mx(A,i,i+1) = -1.0; for (i=1;i<n;i++) mx(A,i+1,i) = -1.0; solver(A,x,b); for (i=1;i<=n;i++) printf("%f\n", x[i] ); }
int main(int argc, char **argv) { static xyc *Z; static nde *N; static double **A, *u; initop(argc, argv); fp2mesh(stdfp(),Z,N); ary2(A,dim1(Z)+1, dim1(Z)+1); ary1(u,dim1(Z)+1); set_A(Z,N,A); set_u(Z,u); esolver(A,u); plt(NULL,NULL,Z,N,u); sleep(1000); return 0; }