## Application 1
##
####data
global L=5;
global A=1;
global B=-30;
global pp=3;
global nel=6;
####mesh
global npe nodg ndof conecg conecdof coordg
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);
[vv obj info iter nf lambda]=sqp(vv0,@energy1,@restrictions);
k=plotfem1(conecg,coordg,pp,conecdof,vv)