Home > demos > schnakenberg > schnakjac.m

schnakjac

PURPOSE

Jacobian for Schnakenberg

SYNOPSIS

function [c,fu,flam,b]=schnakjac(p,u,lam)

DESCRIPTION

 Jacobian for Schnakenberg

CROSS-REFERENCE INFORMATION

This function calls:

This function is called by:

schnakinit init-routine with parameters:

SOURCE CODE

0001 function [c,fu,flam,b]=schnakjac(p,u,lam) 
0002 % Jacobian for Schnakenberg
0003 u=pdeintrp(p.points,p.tria,u);d1=1; d2=60; c=[d1;0;0;d1;d2;0;0;d2]; 
0004 f1u=-1+2*u(1,:).*u(2,:);f1v=u(1,:).^2; f2u=-2*u(1,:).*u(2,:); f2v=-u(1,:).^2; 
0005 fu=[f1u; f2u; f1v; f2v]; flam=[0;1]; b=0;