## Application 1
##
####data
global L=5;
global A=0;
global B=0;
global pp=2;
global nel=4;
####mesh
global npe nodg ndof conecg conecdof coordg rr
npe=2;
nodg=nel+1;
ndof=nel*pp+1;
conecg=[1 2];
conecdof=[1:pp+1];
for i=2:nel
 conecg=[conecg;i,i+1];
 conecdof=[conecdof;(i-1)*pp+1:i*pp+1];
end
coordg=zeros(1,nodg);
hsize=L/nel;
for i=2:nodg
 coordg(i)=(i-1)*hsize;
end
####
vv0=zeros(ndof,1);
for rr=2:-0.05:1.2
[vv obj info iter nf lambda]=sqp(vv0,@energy1nl,@restrictions);
k=plotfem1(conecg,coordg,pp,conecdof,vv/max(vv));
vv0=vv;
end