File: pendc3.sci

package info (click to toggle)
scilab 2.6-4
  • links: PTS
  • area: non-free
  • in suites: woody
  • size: 54,632 kB
  • ctags: 40,267
  • sloc: ansic: 267,851; fortran: 166,549; sh: 10,005; makefile: 4,119; tcl: 1,070; cpp: 233; csh: 143; asm: 135; perl: 130; java: 39
file content (27 lines) | stat: -rw-r--r-- 664 bytes parent folder | download | duplicates (3) 
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
 
// Copyright INRIA

// x=l(u)=u-2*k3*u*sin(k1*sqrt(1+4*k3^2*u^2))/sqrt(1+4*k3^2*u^2)/k2
// y=m(u)=k3*u^2+sin(k1*sqrt(1+4*k3^2*u^2))/sqrt(1+4*k3^2*u^2)/k2

// xx(u)=l(u)-x
function res=xx(u)
S=sqrt(1+4*k3^2*u^2);
res=u-2*k3*u*sin(k1*S)/S/k2-vvx;

function res=fx(x,y)
vvx=x; u=fsolve(uux0,xx,0.01);
S=sqrt(1+4*k3^2*u^2);
res=-S^3*k2/(-S^3*k2+2*k3*sin(k1*S)*S^2+8*k3^3*u^2*cos(k1*S)*k1*S-..
    8*k3^3*u^2*sin(k1*S));
uux0=res;

// yy(u)=m(u)-y
function res=yy(u)
S=sqrt(1+4*k3^2*u^2);
res=k3*u^2+sin(k1*S)/S/k2-vvy;

function res=fy(x,y)
vvy=y; u=fsolve(uuy0,yy,0.01);
S=sqrt(1+4*k3^2*u^2);
res=-S^3*k2/k3/u/(S^3*k2+2*cos(k1*S)*k1*k3*S-2*k3*sin(k1*S))/2;
uuy0=res;