Boundary+Processing 源代码

function [phix,phiy,neff] = wgmodes (lambda, guess, nmodes, dx, dy, varargin); % This function computes the two transverse magn…

function L=linsht(F1,F2,a,b,alpha,beta,M) %Input - F1 and F2 are the systems of first-order equations % representing …

function S=csfit(X,Y,dx0,dxn) %Input - X is the 1xn abscissa vector % - Y is the 1xn ordinate vector % - dxo =…

function B = BoundMirrorExpand(A) % Expand the matrix using mirror boundary condition % % for example % % A = [ % 1 2 3 1…

function B = BoundMirrorEnsure(A) % Ensure mirror boundary condition % % The number of rows and columns of A must be greater than …

% add some comments here % describe clearly why we need the matrix 'index' % the disadvantange of this method here is to % waste…

#include "cbox5.h" void Input5(Xi,Yi,X,Y,Code,Bc) float Xi[21],Yi[21],X[102],Y[102],Bc[101]; int Code[101]; { FILE *infi…

#include "cbox1.h" void Input1(X,Y,Xi,Yi,Code,Bc) float X[102], Y[102],Xi[21],Yi[21],Bc[101]; int Code[101]; { FILE *infile, *o…

function [t, y] = bvp2ode(g0, g1, g2, tspan, bc, n) % Numerical solution y of the boundary value problem % y'' = g0(t) + g1(t)…

function [t, y] = bvp2ode(g0, g1, g2, tspan, bc, n) % Numerical solution y of the boundary value problem % y'' = g0(t) + g1(t)…