vkf

 von karman-plate as a 10 component system including some regularization
 i.e., aux variables like u5 defined as u5=(1-del Lap)^(-1)\pa_x^2 u1

0001 function [c,a,f,b]=vkf(p,u,lam) 
0002 % von karman-plate as a 10 component system including some regularization
0003 % i.e., aux variables like u5 defined as u5=(1-del Lap)^(-1)\pa_x^2 u1
0004 del=0.05; % regularization parameters of aux varia
0005 c=zeros(400,1); 
0006 % first column of C as written in manual
0007 c(1)=1; c(4)=1; c(5)=lam; c(17)=1; c(24)=1; c(26)=0.5;c(27)=0.5; 
0008 c(45)=1; c(48)=1; % second column
0009 c(89)=1; c(92)=1; c(109)=1; c(116)=1; c(118)=0.5; c(119)=0.5; % third
0010 c(133)=1; c(136)=1;     % 4th
0011 c(177)=del; c(180)=del; % 5th
0012 c(221)=del; c(224)=del; 
0013 c(265)=del; c(268)=del; 
0014 c(309)=del; c(312)=del; 
0015 c(353)=del; c(356)=del; 
0016 c(397)=del; c(400)=del;   % to 10th
0017 a=zeros(100,1); a(11)=1; a(33)=1; a(45)=1; a(56)=1; 
0018 a(67)=1; a(78)=1; a(89)=1; a(100)=1;
0019 b=0; 
0020 u5=pdeintrp(p.points,p.tria,u(4*p.np+1:5*p.np)); 
0021 u6=pdeintrp(p.points,p.tria,u(5*p.np+1:6*p.np)); 
0022 u7=pdeintrp(p.points,p.tria,u(6*p.np+1:7*p.np)); 
0023 u8=pdeintrp(p.points,p.tria,u(7*p.np+1:8*p.np)); 
0024 u9=pdeintrp(p.points,p.tria,u(8*p.np+1:9*p.np)); 
0025 u10=pdeintrp(p.points,p.tria,u(9*p.np+1:10*p.np)); 
0026 f1=-(u5.*u9-2*u7.*u10+u6.*u8); 
0027 f2=u5.*u6-u7.^2;  zv=zeros(1,p.nt); 
0028 f=[zv;f1;zv;f2;zv;zv;zv;zv;zv;zv];