pde_14_eq.c running
Given a4(x)*uxxxx(x) + a0(x)*u(x) = F(x)
  F(x) = a4(x)*96  a0(x)*(4 x^4 + 2 x  + 1) 
  a4(x) = x+1      a0(x) = x^2 -1
with 2 boundary conditions, at least one Dirichlet
xmin+hx<x<xmax-hx  inside Boundaries
Analytic solution u(x) = 4 x^4 + 2 x + 1
 
xmin=0.2, xmax=3.3
nx=10
x grid 
xg[0]=0.2, soln=1.4064
xg[1]=0.5444444444444444, soln=2.4403475689681446
xg[2]=0.8888888888888888, soln=5.274958085657674
xg[3]=1.2333333333333332, soln=12.721782716049377
xg[4]=1.5777777777777775, soln=28.943660783417148
xg[5]=1.9222222222222218, soln=59.454719768327955
xg[6]=2.2666666666666666, soln=111.12037530864195
xg[7]=2.6111111111111107, soln=192.15733119951219
xg[8]=2.9555555555555553, soln=312.133579393385
xg[9]=3.3, soln=481.9684
 
 
 
discreteize stiffness matrix and RHS
k computed stiffness matrix
f computed forcing function
kf build took 0.0 seconds
solve 10 equations in 10 unknowns
equations solved
simeq took 0.0010 seconds
 
check_soln on known solution
maxerr=6.361915438901633E-10, rmserr=2.9295132941740794E-10, avgerr=1.8701129533837957E-10
check_soln based on numeric approximation of PDE
maxerr=8.862116374075413E-9, rmserr=3.364380932288791E-9, avgerr=2.0319212978847645E-9
 
ug computed,  Ua analytic,  error  
ug[0]=1.4064,  Ua=1.4064,  err=0.0
ug[1]=2.440347706917862,  Ua=2.4403475689681446,  err=1.3794971742697726E-7
ug[2]=5.2749582878794286,  Ua=5.274958085657674,  err=2.0222175489692518E-7
ug[3]=12.721782926337312,  Ua=12.721782716049377,  err=2.1028793462107842E-7
ug[4]=28.943660963292643,  Ua=28.943660783417148,  err=1.798754958315385E-7
ug[5]=59.454719896388404,  Ua=59.454719768327955,  err=1.2806044935587124E-7
ug[6]=111.12037537919105,  Ua=111.12037530864195,  err=7.05490919017393E-8
ug[7]=192.15733122099945,  Ua=192.15733119951219,  err=2.1487267076736316E-8
ug[8]=312.13357938720833,  Ua=312.133579393385,  err=-6.17666273683426E-9
ug[9]=481.9684,  Ua=481.9684,  err=0.0
                         maxerr=2.1028793462107842E-7, avgerr=9.566083738477005E-8
total took 0.0090 seconds