lcombustion.m

来自「LiScNLS is a Matlab application for the 」· M 代码 · 共 11 行

function [p,dp,q,r,alp,dom,x,wint,wipr]=Lcombustion(n)
dom=[0,1];alpha=(dom(2)-dom(1))/2;beta=(dom(2)+dom(1))/2;M=6;
e1=ones(n,1)/2;e2=e1;e2(2)=1;X=spdiags([e1,e2],[-1,1],n,n);
X=alpha*X+beta*eye(n);
p=expm(-M*X);
dp=-M*p;
q=0;
r=p;
alp=[M,-1;1,0];
[x,wint]=pd(n,dom);
wipr=exp(-M*x);