LaguerreEig: f=leucemia3d(t,u) - File Exchange

Code covered by the BSD License

Highlights from
LaguerreEig

D=deriv(n)
P=param physiological values of the parameters
X=mult(n)
[Ps,Pu,L,V]=proj(A) calculates the projection operators, the eigenvalues
[x,A,f,kont]=jacob(x) calculates the critical points x and the corresponding Jacobian A
[x,coef,sol]=leuclag1(n,bas,P...
[x,w,T,iT,X,D]=lagtool(n,b) Tools for Laguerre spectral method
[x,w]=pd(n) Nodes and weights for Gauss-Laguerre quadrature
events_leuc( t,u )
f=leucemia3d(t,u)
laguer(a, x) http://dip.sun.ac.za/~weideman/research/mfiles/laguer.m
main_fermi
t=x2t(n,x)
airy_LagEig.m
anharmonic.m
hydrogen.m
laguerrestartup.m
main_leuc.m
marletta.m
pryce.m
siyyam.m
siyyam1.m
work2.m
View all files

from LaguerreEig by Damian Trif
Calculates numerical eigenvalues and eigenfunctions of Schroedinger problems on (0,infinity).

f=leucemia3d(t,u)

function f=leucemia3d(t,u)
% the r.h.s. f of the system u'=f(u)

P=param;f=zeros(3,1);
a=P(1,1);b=P(1,2);c=P(1,3);
A=P(2,1);B=P(2,2);C=P(2,3);
g=P(3,1);G=P(3,2);h=P(3,3);
f(1)=(a*(u(1)+u(2))/(1+b*(u(1)+u(2)+u(3)))/(u(1)+u(2)+g*u(3))-c)*u(1);
f(2)=(A*(u(1)+u(2))/(1+B*(u(1)+u(2)+u(3)))/(u(1)+u(2)+G*u(3))-C)*u(2);
f(3)=(a*u(3)/(1+b*(u(1)+u(2)+u(3)))/(u(3)+h*(u(1)+u(2)))-c)*u(3);

Highlights from LaguerreEig

Highlights from
LaguerreEig