Displaying 1 result from an estimated 1 matches for "meshpts".
Did you mean:
meshes
2001 Nov 05
1
Why doesn't outer work?
...trix
I<-diag(n);
M<-matrix(0,n,n);
P<-matrix(0,n,n);
B<-matrix(0,n,n);
P<-t(outer(y,y,pxy))
B<-t(outer(y,y,fxy))
M<-P+B
K<-M;
M<-(U-L)*M/n;
P<-(U-L)*P/n;
B<-(U-L)*B/n;
return(list(matrix=M,kernel=K,meshpts=y,Pmatrix=P,Bmatrix=B,Imatrix=I));
}
bigmatrix<-function(n) {
# upper and lower integration limits
L<-minsize; U<-maxsize;
# boundary points b and mesh points y
b<-L+c(0:n)*(U-L)/n;
y<-0.5*(b[1:n]+b[2:(n+1)]);
# loop to construct the matrix
I<-diag(n);
M<-m...