/usr/share/openturns/validation/Cas2.txt is in openturns-validation 1.5-7build2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 | > restart:
> N:=3;
> G:=x[1]-sum(x[i]^2/i,i=2..N):
N := 3
> Gu:=expand(subs(x[1]=1/2+1/10*u[1],seq(x[i]=1/5+1/10*u[i],i=2..N),G)):
> L:=1/2*sum(u[i]^2,i=1..N)+lambda*Gu:
> sol:=fsolve({seq(diff(L,u[i]),i=1..N),diff(L,lambda)},{seq(u[i],i=1..N
> ),lambda});
sol := {u[2] = 1.423848987, u[1] = -4.158620876,
u[3] = 0.7671756198, lambda = 41.58620876}
> subs(sol,[seq(u[i],i=1..N),lambda]);
[-4.158620876, 1.423848987, 0.7671756198, 41.58620876]
> with(plots):
> with(plottools):
> p1:=implicitplot3d(-Gu,u[1]=-5..0,u[2]=-10..10,u[3]=-10..10,grid=[31,3
> 1,31],style=PATCHNOGRID,axes=framed,scaling=constrained):
> p2:=pointplot3d([subs(sol,[seq(u[i],i=1..N)]),[seq(0,i=1..N)]],symbol=
> circle,symbolsize=20,color=red):
> p3:=spacecurve(subs(sol,[seq(u[i]*t,i=1..N)]),t=0..1,color=black):
> p4:=sphere([seq(0,i=1..N)],sqrt(subs(sol,sum(u[i]^2,i=1..N)))):
> display3d([p1,p2,p3,p4]);
Warning, the name arrow has been redefined
Warning, the name arrow has been redefined
>
|