Pure Iron V4 potential

Pure Iron V3 potential

Pure Iron V2 potential

Pure Iron V1 potential

v4 potential predict cross-slip at Si precipitate

Screw dislocation cross slip in front of Ni Precipitate

It cross slip also for Pure iron

纯铁的变形

不好的计算结果

Deformation simulation

Some MD results using my eam potential

deform script

body {background-image: url(http://4.bp.blogspot.com/_A-VFi-3eGjI/SMJUvjPxbmI/AAAAAAAAALA/hC7Sc2-CiV0/s1600-h/1220688472156-space.jpg);
background-repeat: repeat;
background-position: center top;
}

# Test of YYM potential for A508 alloy system

units metal
boundary p p p

lattice sc 1.0
atom_style atomic

read_data yymScrew.lat
region upper block 100. INF INF INF INF INF
region lower block INF 20. INF INF INF INF
group up region upper
group low region lower
group act subtract all up low
group bnd union up low

pair_style eam/fs
pair_coeff * * yym.eam Si P Mn Fe Ni Cu

neighbor 0.3 bin

compute 1 all coord/atom 2.48
compute 2 all pe/atom
compute 3 all stress/atom

dump 1 all custom 1000 outyym tag type x y z c_1 c_2 c_3[5]

min_style cg
minimize 1.0e-6 0.001 10000 100000

velocity bnd set 0. 0. 0.
fix 1 act nvt 600.0 600.0 100.0
fix 2 bnd setforce 0. 0. 0.
restart 5000 poly.%.1 poly.%.2
timestep 0.001
run 5000

displace_atoms up move 0. 0. 1. units box
displace_atoms low move 0. 0. -1. units box
run 5000
displace_atoms up move 0. 0. 1. units box
displace_atoms low move 0. 0. -1. units box
run 5000
displace_atoms up move 0. 0. 1. units box
displace_atoms low move 0. 0. -1. units box
run 5000
displace_atoms up move 0. 0. 1. units box
displace_atoms low move 0. 0. -1. units box
run 5000
displace_atoms up move 0. 0. 1. units box
displace_atoms low move 0. 0. -1. units box
run 5000

write_restart restart.equil

cascade script,

# Test of YYM potential for A508 alloy system

units metal
boundary p p p

lattice none
atom_style atomic

read_data modScrew.lat
pair_style eam/fs
pair_coeff * * yym.eam Si P Mn Fe Ni Cu

neighbor 0.3 bin

compute 1 all coord/atom 2.48
compute 2 all pe/atom
compute 3 all stress/atom
fix 1 all nve

dump 1 all custom 1000 outyym tag type x y z c_1 c_2 c_3[5]

min_style cg
minimize 1.0e-6 0.001 10000 100000
set atom 1 vx 130000.7643
timestep 0.00001
run 5000

google hacking

google hacking其实并算不上什么新东西,当时并没有重视这技术,认为webshell什么的,并无太大实际用途.google hacking其实并非



如此简单...



google hacking的简单实现

使用google中的一些语法可以提供给我们更多的信息(当然也提供给那些习惯攻击的人更多他们所想要的.),下面就来介绍一些常用的语法.

intext:

这个就是把网页中的正文内容中的某个字符做为搜索条件.例如在google里输入:intext:动网.将返回所有在网页正文部分包含"动网"的网页



.allintext:使用方法和intext类似.



intitle:

和上面那个intext差不多,搜索网页标题中是否有我们所要找的字符.例如搜索:intitle:安全天使.将返回所有网页标题中包含"安全天使"的网



页.同理allintitle:也同intitle类似.



cache:

搜索google里关于某些内容的缓存,有时候也许能找到一些好东西哦.



define:

搜索某个词语的定义,搜索:define:hacker,将返回关于hacker的定义.



filetype:

这个我要重点推荐一下,无论是撒网式攻击还是我们后面要说的对搜索指定类型的文件.例如输入



:filetype:doc.将返回所有以doc结尾的文件URL.当然如果你找.bak、.mdb或.inc也是可以的,获得的信息也许会更丰富



info:

查找指定站点的一些基本信息.



inurl:

搜索我们指定的字符是否存在于URL中.例如输入:inurl:admin,将返回N个类似于这样的连接:http://www.xxx.com/xxx/admin,用来找管理员登陆的URL不错.allinurl也同inurl类似,可指定多个字符.



link:

例如搜索:inurl:www.4ngel.net可以返回所有和www.4ngel.net做了链接的URL.



site:

这个也很有用,例如:site:www.4ngel.net.将返回所有和4ngel.net这个站有关的URL.



对了还有一些*作符也是很有用的:

+ 把google可能忽略的字列如查询范围

- 把某个字忽略

~ 同意词

. 单一的通配符

* 通配符，可代表多个字母

"" 精确查询



下面开始说说实际应用



以下内容均在google上搜索,对于一个居心叵测的攻击者来说,可能他最感兴趣的就是密码文件了.而google正因为其强大的搜索能力往往会



把一些敏感信息透露给他们.用google搜索以下内容:

intitle:"index of" etc

intitle:"Index of" .sh_history

intitle:"Index of" .bash_history

intitle:"index of" passwd

intitle:"index of" people.lst

intitle:"index of" pwd.db

intitle:"index of" etc/shadow

intitle:"index of" spwd

intitle:"index of" master.passwd

intitle:"index of" htpasswd

"# -FrontPage-" inurl:service.pwd

有时候因为各种各样的原因一些重要的密码文件被毫无保护的暴露在网络上,如果被别有用心的人获得,那么危害是很大的



同样可以用google来搜索一些具有漏洞的程序,例如ZeroBoard前段时间发现个文件代码泄露漏洞,可以用google来找网上使用这套程序的站



点:

intext:ZeroBoard filetype:php

或者使用:

inurlutlogin.php?_zb_path= site:.jp

来寻找我们所需要的页面.phpmyadmin是一套功能强大的数据库*作软件,一些站点由于配置失误,导致我们可以不使用密码直接对phpmyadmin进行*作.我们可以用google搜索存在这样漏洞的程序URL:

intitle:phpmyadmin intext:Create new database



还记http://www.xxx.com/_vti_bin/..%5 ... ystem32/cmd.exe?dir吗?用google找找，你也许还可以找到很多



古董级的机器。同样我们可以用这个找找有其他cgi漏洞的页面。

allinurl：winnt system32



前面已经简单的说过可以用google来搜索数据库文件,用上一些语法来精确查找能够获得更多东西(access的数据库,mssql、mysql的连接文件等等).举个例子示例一下:

allinurl:bbs data

filetype:mdb inurl:database

filetype:inc conn

inurl:data filetype:mdb

intitle:"index of" data //在一些配置不正确的apache+win32的服务器上经常出现这种情况,和上面的原理一样,我们还可以用google来找后台.



利用google完全是可以对一个站点进行信息收集和渗透的，下面我们用google对特定站点进行一次测试。

首先用google先看这个站点的一些基本情况(一些细节部分就略去了):

site:xxxx.com

从返回的信息中，找到几个该校的几个系院的域名：

http://a1.xxxx.com

http://a2.xxxx.com

http://a3.xxxx.com

http://a4.xxxx.com

顺便ping了一下，应该是在不同的服务器.学校一般都会有不少好的资料，先看看有什么好东西没



site:xxxx.com filetype:doc

得到N个不错的doc。



先找找网站的管理后台地址：

site:xxxx.com intext:管理

site:xxxx.com inurl:login

site:xxxx.com intitle:管理

超过获得2个管理后台地址：

http://a2.xxxx.com/sys/admin_login.asp

http://a3.xxxx.com:88/_admin/login_in.asp



还算不错，看看服务器上跑的是什么程序：

site:a2.xxxx.com filetype:asp

site:a2.xxxx.com filetype:php

site:a2.xxxx.com filetype:aspx

site:a3.xxxx.com filetype:asp

site:.......

......

a2服务器用的应该是IIS，上面用的是asp的整站程序，还有一个php的论坛

a3服务器也是IIS，aspx+asp。web程序都应该是自己开发的。有论坛那就看看能不能遇见什么公共的FTP帐号什么的：

site:a2.xxxx.com intext:ftp://*:*

没找到什么有价值的东西。再看看有没有上传一类的漏洞：

site:a2.xxxx.com inurl:file

site:a3.xxxx.com inurl:load

在a2上发现一个上传文件的页面：

http://a2.xxxx.com/sys/uploadfile.asp

用IE看了一下，没权限访问。试试注射，

site:a2.xxxx.com filetype:asp

得到N个asp页面的地址，体力活就让软件做吧，这套程序明显没有对注射做什么防范，dbowner权限，虽然不高但已足矣，back a shell不太喜欢，而且看起来数据库的个头就不小，直接把web管理员的密码暴出来再说，MD5加密过。一般学校的站点的密码都比较有规律，通常都是域名+电话一类的变形，用google搞定吧。

site:xxxx.com //得到N个二级域名

site:xxxx.com intext:*@xxxx.com //得到N个邮件地址，还有邮箱的主人的名字什么的

site:xxxx.com intext:电话 //N个电话

把什么的信息做个字典吧，挂上慢慢跑。过了一段时间就跑出4个帐号，2个是学生会的，1个管理员，还有一个可能是老师的帐号。登陆上去：

name：网站管理员

pass：a2xxxx7619 //说了吧，就是域名+4个数字

要再怎么提权那就不属于本文讨论访问了，呵呵，到此为止。



这段时间在国外的一些google hack的研究站点看了看，其实也都差不多是一些基本语法的灵活运用，或者配合某个脚本漏洞，主要还是靠个人的灵活思维。国外对于google hack方面的防范也并不是很多，所以大家还是点到为止，不要去搞破坏拉，呵呵。对于一些在win上跑apache的网管们应该多注意一下这方面，一个 intitle:index of就差不多都出来了

1.查找利用php webshell



intitle:"php shell*" "Enable stderr" filetype:php



(注: intitle—网页标题 Enable stderr—UNIX标准输出和标准错误的缩写filetype—文件类型)。搜索结果中，你能找到很多直接在机器上执



行命令的web shell来。如果找到的PHPSHELL不会利用，如果你不熟悉UNIX，可以直接看看LIST，这里就不详细说了，有很多利用价值。要说明



的是，我们这里搜索出来的一些国外的PHPSHELL上都要使用UNIX命令，都是system调用出来的函数(其实用百度及其他搜索引擎都可以，只是填写搜索的内容不同)。这个PHPWEBSHELL是可以直接Echo(Unix常用命令)。一句话就把首页搞定了:



echo "召唤" > index.jsp



现在看看首页，已经被我们改成: "召唤" 了。



我们也可以用WGET上传一个文件上去(比如你要替换的叶子吧)。然后execute Command输入 cat file > index.html or echo "" > file



echo "test" >> file



这样一条条打出来，站点首页就成功被替换了。同样的也可以



uname -a;cat /etc/passwd



不过有点要注意，有些WEBSHELL程序有问题，执行不了的，

velocity of a 50eV iron

( 1.e10/1.e12 )*sqrt( 2*50*1.602e-19/1.673e-27/56 )=130.7643

the unit is Angstroms/picosecond

upload example

Send this file:

1.5GHz的VIA C7-D处理器, 60$

觉得OLPC或者Eee PC还是太贵？这里有更便宜的。国外网络商店已经开始销售一款名为“gOS开发套装”的产品，包含主板、CPU和操作系统的售价仅为60美元。

该套装包括VIA PC2500E主板，搭载了1.5GHz的VIA C7-D处理器。主板内置Realtek ALC655 6声道音频，10/100M网卡，2个PATA接口，2个SATA接口，支持VGA视频输出，4个USB接口。只要添加内存、硬盘和电源就能成为一套完整的系统。而其附赠的gOS系统基于Ubuntu，预装了几乎所有的Google软件服务以及一些必备的网络软件，是一套非常实用的网络化操作系统。

虽然目前并没有gOS得到Google官方支持的任何消息，但无疑这种超低价的Linux软硬件套装能够吸引大量的开发人员投入到gOS的开发工作中来。在开发人才的争夺战中，Google不费一兵一卒再下一城。


60美元 VIA主板+CPU+gOS操作系统


60美元 VIA主板+CPU+gOS操作系统

MCAD

class="org.alfresco.module.demoscripts.SimpleWebScript"
parent="webscript">


The World's Simplest Webscript
Hands back a little bit of JSON
/demo/simple
none
argument

VNC set

Virtual Network Computing (VNC) allows a computer to be seen and accessed remotely from other computers, who can see the screen and control the system using the keyboard and mouse. VNC is useful for remote technical support or remote access to personal/work computers. tip: If you need to access a Mac remotely, read Apple Remote Desktop instead.
目录
[隐藏]

* 1 Terminology
o 1.1 General Security
* 2 Enabling VNC connections ~ Server setup
o 2.1 Using the default vnc servers
+ 2.1.1 using GNOME / Ubuntu (vino)
+ 2.1.2 using KDE / Kubuntu
+ 2.1.3 using XFCE / Xubuntu
o 2.2 VNC Server with Login Screen via GDM
o 2.3 Tightvncserver Server with Login Screen Via GDM
o 2.4 FreeNX
o 2.5 Tunnel VNC through SSH
* 3 Enabling VNC connections ~ Client setup
o 3.1 Ubuntu clients
+ 3.1.1 Terminal Server Client
+ 3.1.2 VNC Viewer
+ 3.1.3 Logging into a Tight VNC server
o 3.2 Windows Clients
* 4 VNC Access over the Internet
o 4.1 Server IP
o 4.2 Router
o 4.3 Firewall
o 4.4 Method 1 ~ Using vncviewer from the command line
o 4.5 Method 2 ~ Ubuntu clients ~ Tunnel over ssh directly
o 4.6 Method 3 ~ Windows or Ubuntu ~ Tunnel over ssh manually
o 4.7 Method 4 ~ Via a web browser (firefox for example)
+ 4.7.1 Server setup
+ 4.7.2 Client setup
+ 4.7.3 Connect
o 4.8 Connecting with a Windows XP client

Terminology

VNC Server = the computer you will connect to (log onto remotely). VNC Client = the computer you will connect with (use to log on to a server). With the default vnc server (vino) or x11vnc you will log into a shared desktop. With tightvnc server or with configuration of GDM you will have a separate session.

* If you are logging into a shared desktop you will not be able to connect if there is no user logged into the system.
* With a separate session you can log in even if either no one is logged into the server, or even if X is not running on the server at all.

Just a quick note about terminology. X sessions or desktops are numbered starting with 0 and is referred to on the command line as :0.

* So the desktop you have when you first log in can be referred to as
* localhost:0
* 172.0.0.1:0
* Your lan ip (ie 192.168.1.25:0)
* Your internet IP (128.220.223.246:0)

Not to be confused with ports. Servers listen on your computer via ports and vnc connections uses 5900 + by default. So the first vnc session is port 5900. The second vnc server will use 5901.

* The default ports can of course be changed.

When you activate your vnc server via System -> Preferences -> Remote Desktop and tic off the "Allow other users to view your desktop" you are :

* Activating a server on your desktop which allows others to connect to your desktop (see "security issues" below).
* Your vnc server is listening for connections on your ip:5900 (ie 192.168.1.25:5900) or your internet IP :5900 (ie 128.220.223.246:5900).

When you start a new vnc server with tightvncsserver (see below) this second session starts on :1

* So the second vnc session is on localhost:1 (192.168.1.25:1 and on ...)
* And the port is :5901 (localhost:5901 or 192.168.1.25:5901)

When you connect to a vnc server you will use the server ip address and :0 or :1 to refer to the vnc session and 192.168.1.25:5900 or 192.168.1.25:5901 for configuration of routers/firewalls.
General Security

Security issuesBy default Ubuntu will allow all connections so at a minimum set a password when you activate Remote desktop. Without one people can watch your desktop from your LAN without any password. So if you are connected directly to the internet, without a router, if you do not set a password, anyone can access your computer. You may also consider:

* Purchasing a router.
* Create a user(s) specifically for vnc sessions without admin (sudo) access.
* Either logging out or locking your desktop if not in use (System -> Preferences -> Screensaver tic off the "Lock screen when screensaver is active"). If you screen is locked your log in password will be required to unlock the screen.
o Kde Configure Desktop -> screen Saver tab -> tic off the "Require password to stop." option.
* Installing firestarter to help configure your firewall.
* Using ssh to log into vnc sessions VNC over SSH

Enabling VNC connections ~ Server setup

To allow other computers to access your desktop sessions, when your logged in, perform the following steps. Note: This is like MS Remote Assistance and only works when your are logged into the computer sharing that instance of your Xserv. See below for creating Xservs on the fly and allowing multiple loggings and XDM/GDM.
Using the default vnc servers
using GNOME / Ubuntu (vino)

System > Preferences > Remote Desktop 'Check' the first two boxes to activate the service: Allow other users to view your desktop (view only) Allow other users to control your desktop (view & control). Below you can set security. The two options are: Ask you for confirmation (ie; someone at the machine must click OK to grant remote access. This will be a problem if you plan on accessing your home machine from work or visa versa, as no one may be there to grant you access.) Require the user to enter this password: This will require a password from anyone trying to remotely connect to your machine. This is ALWAYS a good idea.
using KDE / Kubuntu

System Settings > Sharing > Desktop Sharing > Create & Manage Invitations ...Choose 'New Personal Invitation...' (you give the invitation by whichever means you prefer (Email, Instant Messaging, Written-Down Note) or 'New Email Invitation...'. The newly created invitation will last a default of 1 hour. tip If you don't see this option ensure that the package 'krfb' (Desktop Sharing for KDE) is installed. Sometimes it is installed but doesn't appear on the menu. If it is installed, type Alt-F2 and enter krfb or type it in the Konsole. You can manually add menu entries by righ clicking the KDE menu icon.
using XFCE / Xubuntu

There is not vnc server installed by default in Xubuntu. The default gnome vnc server is vino and you can install this package, X11vnc, or tightvncserver.
VNC Server with Login Screen via GDM

This method is somewhat more complicated, but when connecting this way you get a login prompt and begin a new session. This also works when no user is logged in and allows multiple parallel loggings. Step 1 Append the following line to /etc/services

$> gksudo gedit /etc/services

vnc 5901/tcp # VNC with GDM

Step 2 Create the following file /etc/xinetd.d/vnc

$> sudo pico /etc/xinetd.d/vnc

service vnc
{
disable = no
socket_type = stream
protocol = tcp
wait = no
user = nobody
server = /usr/bin/Xvnc
server_args = -inetd :1 -query localhost -broadcast -once -fp /usr/share/X11/fonts/misc/ -securitytypes=none -desktop=vnc://MyDesktop/
}

* Note: In 6.10 the default fount's were not found. The -fp /usr/share/X11/fonts/misc/ line should resolve this error.
* Note: -desktop=vnc://MyDesktop/ is the title that appears to the user when connecting.
* Note: -query localhost is optional, it may fix a problem when clients get "connection unexpectedly closed"

Step 3 Enable XDMCP in your login configuration to allow remote logins to GDM (the gnome login screen). Edit /etc/gdm/gdm.conf

$> sudo pico /etc/gdm/gdm.conf

find the section [xdmcp] and set the enable to true:

[xdmcp]
....
Enable=true

uncomment this line:

RemoteGreeter=/usr/lib/gdm/gdmlogin

Step 4 Stop and restart Xinetd

$> sudo /etc/init.d/xinetd restart

Problems

* If you cant connect check your router/port forwarding, firewall, or try running the following to start a session without Xinetd to verify Xvnc is working.

$> Xvnc :1 -fp /usr/share/fonts/X11/misc/

* Make sure you connect to the proper port, in this case vnc://localhost:5901. In some clients this is set by choosing display 1.

Tightvncserver Server with Login Screen Via GDM

This method will give you an independent desktop once you log in (ie independent of session started when you logged in via GDM). Step 1 - Install tightvncserver

sudo apt-get install vnc-common tightvncserver

Step 2 - Edit vncserver script

sudo vim /usr/bin/vncserver

First, you want a valid X11 font path for Xvnc. Add these lines for Feisty/Gutsy (you will see a fonts section with a number of font paths commented out):

$fontPath = join ',',qw(
/usr/share/fonts/X11/misc
/usr/share/fonts/X11/100dpi/:unscaled
/usr/share/fonts/X11/75dpi/:unscaled
/usr/share/fonts/X11/Type1
/usr/share/fonts/X11/100dpi
/usr/share/fonts/X11/75dpi
);

You may add additional font paths as needed.

* see this link for font paths with earlier versions of Ubuntu

Optional: Set display size and color depth (make this whatever you want):

$geometry = "1280x1024";
$depth = 16;

Optional: Uncomment the line:

$colorPath = "/usr/lib/X11/rgb";

Save file and close vim:

:wq

Step 3 - Run vncserver for the first time

vncserver :1

The first time you run the server you will be asked to provide a name and password :

ubuntu@ubuntu:~$ vncserver :1

You will require a password to access your desktops.

Password: #Enter your desired password here
Verify: #Confirm Password

New 'X' desktop is ubuntu:1

Starting applications specified in /etc/X11/Xsession
Log file is /home/ubuntu/.vnc/ubuntu:1.log

ubuntu@ubuntu:~$

This will create a new directory in your home directory ~/.vnc

* To change you password later use vncpasswd

vncpasswd ~/.vnc/passwd

Step 4 - Edit your VNC startup script We may want to edit the VNC (X) startup script in ~/.vnc/xstartup You can use these if you like : Gnome

xrdb $HOME/.Xresources
xsetroot -solid navy # Choose your color
x-window-manager &
gnome-panel 2> /dev/null &
xterm &

XFCE

xrdb $HOME/.Xresources
xfwm4 2> /dev/null &
xfce4-panel 2> /dev/null &
xfce4-terminal &

KDE

xrdb $HOME/.Xresources
xsetroot -solid navy # Choose your color
x-terminal-emulator -geometry 80x24+10+1- -ls -title "$VNCDESKTOP Desktop" &
x-window-manager &
kicker 2> /dev/null &

Step 5 - Restart the VNC server

killall Xtightvnc
vncserver

Note that you can vary the screen size, depth, and number when starting vncserver:

vncserver -geometry 1288x1024 -depth 24 :3

* See vncserver manpage for additional options

FreeNX

See Ubuntu Wiki How to FreeNX
Tunnel VNC through SSH

If you wish to tunnel over ssh you need to install, setup, and secure the openssh server.

* Ubuntu Wiki How to SSH
* Ubuntu Wiki, Advanced SSH (security) settings

Enabling VNC connections ~ Client setup
Ubuntu clients
Terminal Server Client

This method works with both the default vino server and the tightvnc server This is the default method in Ubuntu and uses a gui (graphical interface). Use this if you are averse to the command line. Go to Applications -> Internet -> Terminal Server client In the General tab :

* Put the server ip in the "Computer" box (ie 192.168.1.25:0 for the default vino server or 192.168.2.25:1 for a tightvnc server)
* Select VNC from the pull down menu in "Protocol"

Hit Connect Enter the password you set on your server.
VNC Viewer

This method uses the command line. Open a terminal an enter vncviewer and you will be asked to enter the ip address and password. If you know the ip address you can use the ip in conjunction, like this :

vncviewer 192.168.1.25:0

Or if you want to get fancy, copy ~/.vnc/passwd from the server to the client (saving it in ~/.vnc/passwd on both the server and client). Now you can connect directly with :

vncviewer 192.168.2.25:0 -passwd ~/.vnc/passwd

* If you like, you can re-name the ~/.vnc/passwd to any name you like and keep one file for each server (each with a unique name).

Logging into a Tight VNC server

If desired you will need to configure your desktop. Ubuntu I do not know how to set the background image on the tightvnc server, but the gnome panel works. XFCE Applications -> Settings -> Desktop Settings

* Tic off the "Allow Xfce to manage the desktop"

KDE The kicker works fine, as with gnome I could not get set the backgound image.
Windows Clients

This method works with Windows 2000 and XP. I am not sure about other versions. Connecting with a windows client is fairly straight forward, all you need to do is download the tight vnc viewer for windows. You then run the viewer and enter the server address and password very similar to connecting from Ubuntu. Tight VNC viewer for windows

* The tight vnc viewer is available without installation ( tightvnc-1.3.9_x86_viewer.zip "Viewer executable, does not require installation")

I have also used Ultra and Real vnc viewers

* You can set up tight, ultra, or real vnc servers on windows and connect just as easily from ubuntu.

VNC Access over the Internet

Strongly consider tunneling over ssh VNC Over SSH , using a router (rather then connecting your server directly to the internet), and configuring your firewall Firestarter. The most difficult part of internet access is configuring your router and firewall ...
Server IP

The IP address of your VNC server is different on a LAN vs an internet connection. The internet IP address is assigned by your internet provider.

* You can check your IP address here (or elsewhere).

The problem can occur if you use DHCP (rather then a static IP address) the internet IP address can change from time to time. The solution is to register at dyndns or other providers. dyndns will provide free service. You can then determine your vnc server address via ping.
Router

You must configure your router to forward the ports. The details vary by router. tip: if you are connecting through a router you'll need to forward port 5900 to the machine you need to connect to. (VNC can use other ports as well. If you have multiple machines you would like to connect to you can forward 5900 to the first, 5901 to the second, 5902 to the third, etc.) If you are using a non-standard port (ie; other than 5900) you will need to specify the port in the connection command. You can simplify this step some what by tunneling though ssh (which also increases security). See the ssh section below.
Firewall

This is very easy to do via firestarter, a gui front end to IP Tables.

* Be sure to configure firestarter to allow pings.

Open the firestarter gui (Applications -> Internet -> Firestarter) In the "Policy" tab, under the "Allow service" section, right click anywhere in the white space. Select "Add rule"

* Under "Port" enter the ports you want to enable (5900 and/or 5901)

Allow "Anyone"
Method 1 ~ Using vncviewer from the command line

* Using tightvncserver - See the Tightvncserver section above.
* If you are using the default vnc server, vino, for a shared desktop, use 5900 or :0

Open a terminal and start vncviewer with this command

vncviewer

Enter the ip address :1 Example:

192.168.1.25:1

Enter the vnc password
Method 2 ~ Ubuntu clients ~ Tunnel over ssh directly

You can use this method with Ubuntu clients. Use the -via flag -via = use ssh authentication. vncviewer -via

vncviewer -via 192.168.1.25 ubuntu:1

Enter ssh password, enter vnc password

* You are given the name of the vncserver by tight vnc when you start it up, see the tight vnc server section above

Method 3 ~ Windows or Ubuntu ~ Tunnel over ssh manually

This is a quick guide and assumes you have a ssh server set up on the vnc server. See this link for a more detailed description : VNC Over SSH The trick is to forward the ports over ssh. In this example I will use 5900 , the default path for the defalut VNC server vino. If you use tightvnc you will need to change the forwarded port to 5901. Step 1 Make the ssh connection :

ssh -fCNT user@192.168.1.25 -L 5901:127.0.0.1:5901

-f = Allows ssh to close after the connection is established.
-C = Use Compression
-N = No commands will be issued
-T = No terminal session will be started

-L = Port forwarding. The terminology is : the trick is we are using 127.0.0.1: for the client. 127.0.0.1 must be used (not localhost or the client ip address)

Step 2 Make the vnc connection. Now we make the vnc connection, but now we use localhost:1 as the server ip.

vncviewer localhost:1

Enter your password. Step 3 To disconnect, close the vnc viewer, and enter killall ssh in the terminal.
Method 4 ~ Via a web browser (firefox for example)

This does not work with the default vnc server, vino. You will need to install and configure the tight vnc server as above. You will need to enable the commercial repositories on both the server and client (or download the debs) :

deb http://archive.canonical.com/ubuntu gutsy commercial
deb-src http://archive.canonical.com/ubuntu gutsy commercial

Server setup

Install by any means vnc-common, tightvncserver, and tightvnc-java

sudo apt-get install vnc-common tightvncserver tightvnc-java

Configure the tight vnc server as above. You will likely want to reduce the resolution as the java applet will run in a firefox window :

vncserver -geometry 800x600 -depth 24 :1

* If you have a large monitor you may be able to increase the server resolution.

The java server will start automatically
Client setup

On an Ubuntu client install by any means sun-java6-jre and sun-java-6-plugin

sudo aptitude install sun-java6-jre sun-java-6-plugin

Allow java : In Firefox Edit -> Preferences Select the "Content" tab, tic off the "Load images automatically" "Enable JavaScript" and "Enable Java" boxes.
Connect

Open Firefox, in the address bar type vnc server ip : 5801 Example:

192.168.1.25:5801

The java applet will start automatically. If you use NoScript (or other java blockers) you will need to allow 192.168.1.25:5801 Click the connect button.

* If you reload the firefox window you will need to log in again.

Connecting with a Windows XP client

You may use any of the above methods with Windows.

* VNC - See the windows section above.
* SSH - For windows I have used both Putty and Cygwin with success. See here for details vnc over ssh Windows clients
* For firefox just be sure to install java and the java plugins on the windows client.

Note: If you tunnel ssh connections on a windows client, you must first ssh method (ie start ssh first, then connect with tightvnc viewer). Note: Cygwin offers a ssh server so you could potentially tunnel vnc connections from a Windows vnc server.
取自"http://wiki.ubuntu.org.cn/index.php?title=UbuntuHelp:VNC&variant=zh-cn"

SSH Tunneling Made Easy

By Frank Wiles

I was surprised to how long it took me to find a good HOWTO on setting up a simple SSH tunnel that I wanted to write up this Quick-Tip.

Using OpenSSH on a Linux/Unix system you can tunnel all of the traffic from your local box to a remote box that you have an account on.

For example I tunnel all of my outbound E-mail traffic back to my personal server to avoid having to change SMTP servers, use SMTP-AUTH, etc. when I am behind firewalls. I find that hotel firewalls, wireless access points, and the other various NATing devices you end up behind while traveling often do not play nice.

To do this I use the following:

ssh -f user@personal-server.com -L 2000:personal-server.com:25 -N

The -f tells ssh to go into the background just before it executes the command. This is followed by the username and server you are logging into. The -L 2000:personal-server.com:25 is in the form of -L local-port:host:remote-port. Finally the -N instructs OpenSSH to not execute a command on the remote system.

This essentially forwards the local port 2000 to port 25 on personal-server.com over, with nice benefit of being encrypted. I then simply point my E-mail client to use localhost:2000 as the SMTP server and we're off to the races.

Another useful feature of port forwarding is for getting around pesky firewall restrictions. For example, a firewall I was behind recently did not allow outbound Jabber protocol traffic to talk.google.com. With this command:

ssh -f -L 3000:talk.google.com:5222 home -N

I was able to send my Google Talk traffic encrypted through the firewall back to my server at home and then out to Google. All I had to do was reconfigure my Jabber client to use localhost as the server and the port 3000 that I had configured.

Hopefully this helps you to better understand SSH tunneling. If you found this page useful, you may also be interested in how to make your SSH connections faster. If you find any errors or have any suggestions regarding this please feel free to E-mail me at frank@revsys.com.

example of xorg.conf

cat /etc/X11/xorg.conf
# xorg.conf (xorg X Window System server configuration file)
#
# This file was generated by dexconf, the Debian X Configuration tool, using
# values from the debconf database.
#
# Edit this file with caution, and see the xorg.conf manual page.
# (Type "man xorg.conf" at the shell prompt.)
#
# This file is automatically updated on xserver-xorg package upgrades *only*
# if it has not been modified since the last upgrade of the xserver-xorg
# package.
#
# If you have edited this file but would like it to be automatically updated
# again, run the following command:
# sudo dpkg-reconfigure -phigh xserver-xorg

Section "Files"
EndSection

Section "InputDevice"
Identifier "Generic Keyboard"
Driver "kbd"
Option "CoreKeyboard"
Option "XkbRules" "xorg"
Option "XkbModel" "pc105"
Option "XkbLayout" "us"
Option "XkbVariant" "dvorak"
Option "XkbOptions" "lv3:ralt_switch"
EndSection

Section "InputDevice"
Identifier "Configured Mouse"
Driver "mouse"
Option "CorePointer"
Option "Device" "/dev/input/mice"
Option "Protocol" "ImPS/2"
Option "ZAxisMapping" "4 5"
EndSection

Section "InputDevice"
Driver "wacom"
Identifier "stylus"
Option "Device" "/dev/input/wacom"
Option "Type" "stylus"
Option "ForceDevice" "ISDV4" # Tablet PC ONLY
EndSection

Section "InputDevice"
Driver "wacom"
Identifier "eraser"
Option "Device" "/dev/input/wacom"
Option "Type" "eraser"
Option "ForceDevice" "ISDV4" # Tablet PC ONLY
EndSection

Section "InputDevice"
Driver "wacom"
Identifier "cursor"
Option "Device" "/dev/input/wacom"
Option "Type" "cursor"
Option "ForceDevice" "ISDV4" # Tablet PC ONLY
EndSection

Section "Device"
Identifier "Intel Corporation 82845G/GL[Brookdale-G]/GE Chipset Integrated Graphics Device"
Driver "intel"
BusID "PCI:0:2:0"
EndSection

Section "Monitor"
Identifier "HKC775XX"
Option "DPMS"
HorizSync 30-75
VertRefresh 60-85
EndSection

Section "Screen"
Identifier "Default Screen"
Device "Intel Corporation 82845G/GL[Brookdale-G]/GE Chipset Integrated Graphics Device"
Monitor "HKC775XX"
DefaultDepth 24
SubSection "Display"
Modes "1280x1024" "1024x768" "800x600" "640x480"
EndSubSection
EndSection

Section "ServerLayout"
Identifier "Default Layout"
Screen "Default Screen"
InputDevice "Generic Keyboard"
InputDevice "Configured Mouse"

# Uncomment if you have a wacom tablet
# InputDevice "stylus" "SendCoreEvents"
# InputDevice "cursor" "SendCoreEvents"
# InputDevice "eraser" "SendCoreEvents"
EndSection

linux ssh access control

This is the setup guide for ssh access for linux hosts

Creating a simple tunnel

Creating a simple tunnel

At the beginning let’s assume that SSH connections are not blocked, but we can receive e-mails only from our local office account. This means that ports 25 and 110 (SMTP and POP) are open only for our local network. But what if we want to be able to receive mail from external servers (EXT_POP_SVR) and send it using external accounts as well (EXT_SMTP_SVR)? Here is where the tunnels enter the stage. We just need to add a parameter while connecting to the REMOTE_HOST via SSH. So, we’re digging the tunnel from our work machine to home with the following command:

worker@LOCAL_HOST:~$ ssh user@IP_NUMBER \\
-L 10025:EXT_SMTP_SVR:25 -L 10110:EXT_POP_SVR:110

Let’s try to decode this. The option -L can be read as “listen on local host”. After a space bar we put the port on which SSH has to listen (in our case respectively 10025 and 10110, but it’s just an example; we can choose different ones but in order to tunnel ports with numbers lower than 1024 we need root privileges). After first colon we point where our REMOTE_HOST should forward the connection. After second colon we state on which port the final server/computer awaits our connection. I emphasize that just after the first colon we give an address relative to REMOTE_HOST. If we want to connect via tunnel with REMOTE_HOST on port 25 (and not further than there) we should use something like:

worker@LOCAL_HOST:~$ ssh user@IP_NUMBER -L 10025:localhost:25

The relation between REMOTE_HOST and… REMOTE_HOST is the following: this is the very same host ;), so it’s localhost for itself.

Let’s get back to our tunnels. We can now configure an e-mail client to connect to an external server. As we’ve mentioned earlier, our LOCAL_HOST listening on port 10025 is now “EXT_SMTP_SVR server listening on port 25”. By analogy, we can say that localhost:10110 is an EXT_POP_SVR:110. In our e-mail client we shall set an incoming server to localhost and port to 10110 - SSH will carry on further communication. Naturally, localhost set on LOCAL_HOST will point to itself. And again by analogy outgoing server should be set to localhost and port to 10025.

In this quite simple (I hope!) way we can bypass not very complex restrictions only. But our network administrator is not sleeping so…

Note that you don’t have to use your IP_NUMBER all the time. If you want to use REMOTE_HOST by its name you need to put a line like this

IP_NUMBER REMOTE_HOST

to your /etc/hosts file (on the machine at work), so that your work machine translates the REMOTE_HOST to IP_NUMBER on the fly. You don’t need to do this only if the REMOTE_HOST is a domain name (or if you’re fine with using your IP_NUMBER all the time). In this text we use IP_NUMBER nevertheless.

Let’s swing into a higher gear

Our present situation has changed. It seems that somebody in our work abuses usage of WWW. Access via HTTP has been cut to the minimum allowed by the boss. Right now we are not interested in how it has been done. Port 22 is still open, but we can’t say the same about port 80 (responsible for HTTP communication). Now when we know how the tunnels work, this is not a big problem for us.

It seems impossible but even the Google search has been blocked. So, to solve this unfortunate situation, we sit down and log onto our REMOTE_HOST like that:

worker@LOCAL_HOST:~$ ssh user@IP_NUMBER \\
-L 10080:www.google.com:80

Just a few clicks in Firefox: Edit -> Preferences -> Connection Settings -> Manual proxy configuration. Here we can put in the HTTP proxy server field: localhost and 10080 for port. If there has been some other values in those fields, we shall write them down. We may need to use them later on. I’ll call those extra settings PROXY_SVR and PROXY_PORT.

Now, we can type http://www.google.com in the address bar and… hip hip hooray – the WWW is wide open again. Don’t be so happy though. Try http://www.altavista.com. And what you’ve got – Google search again. Well, what else can we expect if we’re forwarding all traffic to www.google.com via SSH! To visit AltaVista we should reconnect again with:

worker@LOCAL_HOST:~$ ssh user@IP_NUMBER \\
-L 10080:www.altavista.com:80

This is rather a pretty non-effective way to work with the Internet, don’t you think?

PARA16: 0.3972 8.6032 3.4948 4.7131

bulk modulus=1.690768e+02GPa,Target=1.700000e+02
Lattice Constan=2.8665, Target=2.8665

>> CalcForce(para,2.82,10,1)
1, -0.15025, -0.15025, -0.15025, -0.18126, -0.18126, -0.18126
2, 0.00000, 0.00000, 0.10042, -0.00000, -0.00000, 0.10740
3, -0.01467, -0.01467, 0.00000, -0.01445, -0.01445, 0.00000
4, 0.00000, 0.00000, -0.10042, -0.00000, -0.00000, -0.10740
5, -0.15025, -0.15025, 0.15025, -0.18126, -0.18126, 0.18126
6, 0.00000, 0.10042, 0.00000, -0.00000, 0.10740, -0.00000
7, -0.01467, 0.00000, -0.01467, -0.01445, 0.00000, -0.01445
8, 0.00000, -0.01072, -0.01072, -0.00000, -0.03716, -0.03716
9, -0.01754, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01072, 0.01072, -0.00000, -0.03716, 0.03716
11, -0.01467, 0.00000, 0.01467, -0.01445, 0.00000, 0.01446
12, 0.00000, -0.10042, 0.00000, -0.00000, -0.10740, -0.00000
13, -0.15025, 0.15025, -0.15025, -0.18126, 0.18126, -0.18126
14, 0.00000, 0.01072, -0.01072, -0.00000, 0.03716, -0.03716
15, -0.01467, 0.01467, 0.00000, -0.01445, 0.01446, 0.00000
16, 0.00000, 0.01072, 0.01072, -0.00000, 0.03716, 0.03716
17, -0.15025, 0.15025, 0.15025, -0.18126, 0.18126, 0.18126
18, 0.10042, 0.00000, 0.00000, 0.10740, -0.00000, -0.00000
19, 0.00000, -0.01467, -0.01467, 0.00000, -0.01445, -0.01445
20, -0.01072, 0.00000, -0.01072, -0.03716, -0.00000, -0.03716
21, 0.00000, -0.01754, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01072, 0.00000, 0.01072, -0.03716, -0.00000, 0.03716
23, 0.00000, -0.01467, 0.01467, 0.00000, -0.01445, 0.01446
24, -0.01072, -0.01072, 0.00000, -0.03716, -0.03716, -0.00000
25, 0.00000, 0.00000, -0.01754, 0.00000, 0.00000, -0.00000
26, 0.02401, 0.02401, 0.02401, -0.01556, -0.01556, -0.01556
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.02401, 0.02401, -0.02401, -0.01556, -0.01556, 0.01556
29, 0.00000, 0.00000, 0.01754, 0.00000, 0.00000, 0.00000
30, -0.01072, 0.01072, 0.00000, -0.03716, 0.03716, -0.00000
31, 0.00000, 0.01467, -0.01467, 0.00000, 0.01446, -0.01445
32, 0.02401, -0.02401, 0.02401, -0.01556, 0.01556, -0.01556
33, 0.00000, 0.01754, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.02401, -0.02401, -0.02401, -0.01556, 0.01556, 0.01556
35, 0.00000, 0.01467, 0.01467, 0.00000, 0.01446, 0.01446
36, -0.10042, 0.00000, 0.00000, -0.10740, -0.00000, -0.00000
37, 0.15025, -0.15025, -0.15025, 0.18126, -0.18126, -0.18126
38, 0.01072, 0.00000, -0.01072, 0.03716, -0.00000, -0.03716
39, 0.01467, -0.01467, 0.00000, 0.01446, -0.01445, 0.00000
40, 0.01072, 0.00000, 0.01072, 0.03716, -0.00000, 0.03716
41, 0.15025, -0.15025, 0.15025, 0.18126, -0.18126, 0.18126
42, 0.01072, -0.01072, 0.00000, 0.03716, -0.03716, -0.00000
43, 0.01467, 0.00000, -0.01467, 0.01446, 0.00000, -0.01445
44, -0.02401, 0.02401, 0.02401, 0.01556, -0.01556, -0.01556
45, 0.01754, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.02401, 0.02401, -0.02401, 0.01556, -0.01556, 0.01556
47, 0.01467, 0.00000, 0.01467, 0.01446, 0.00000, 0.01446
48, 0.01072, 0.01072, 0.00000, 0.03716, 0.03716, -0.00000
49, 0.15025, 0.15025, -0.15025, 0.18126, 0.18126, -0.18126
50, -0.02401, -0.02401, 0.02401, 0.01556, 0.01556, -0.01556
51, 0.01467, 0.01467, 0.00000, 0.01446, 0.01446, 0.00000
52, -0.02401, -0.02401, -0.02401, 0.01556, 0.01556, 0.01556
53, 0.15025, 0.15025, 0.15025, 0.18126, 0.18126, 0.18126
>> CalcForce(para,2.84,11,1)
1, -0.13232, -0.13232, -0.13232, -0.14331, -0.14331, -0.14331
2, 0.00000, 0.00000, 0.10109, 0.00000, 0.00000, 0.12209
3, -0.00586, -0.00586, 0.00000, -0.01417, -0.01417, -0.00000
4, 0.00000, 0.00000, -0.10109, 0.00000, 0.00000, -0.12209
5, -0.13232, -0.13232, 0.13232, -0.14331, -0.14331, 0.14331
6, 0.00000, 0.10109, 0.00000, 0.00000, 0.12209, 0.00000
7, -0.00586, 0.00000, -0.00586, -0.01417, -0.00000, -0.01417
8, 0.00000, -0.01627, -0.01627, 0.00000, -0.03656, -0.03656
9, -0.00210, 0.00000, 0.00000, -0.00000, -0.00000, -0.00000
10, 0.00000, -0.01627, 0.01627, 0.00000, -0.03656, 0.03656
11, -0.00586, 0.00000, 0.00586, -0.01417, -0.00000, 0.01417
12, 0.00000, -0.10109, 0.00000, 0.00000, -0.12209, 0.00000
13, -0.13232, 0.13232, -0.13232, -0.14331, 0.14331, -0.14331
14, 0.00000, 0.01627, -0.01627, 0.00000, 0.03656, -0.03656
15, -0.00586, 0.00586, 0.00000, -0.01417, 0.01417, -0.00000
16, 0.00000, 0.01627, 0.01627, 0.00000, 0.03656, 0.03656
17, -0.13232, 0.13232, 0.13232, -0.14331, 0.14331, 0.14331
18, 0.10109, 0.00000, 0.00000, 0.12209, 0.00000, 0.00000
19, 0.00000, -0.00586, -0.00586, -0.00000, -0.01417, -0.01417
20, -0.01627, 0.00000, -0.01627, -0.03656, 0.00000, -0.03656
21, 0.00000, -0.00210, 0.00000, -0.00000, -0.00000, -0.00000
22, -0.01627, 0.00000, 0.01627, -0.03656, 0.00000, 0.03656
23, 0.00000, -0.00586, 0.00586, -0.00000, -0.01417, 0.01417
24, -0.01627, -0.01627, 0.00000, -0.03656, -0.03656, 0.00000
25, 0.00000, 0.00000, -0.00210, -0.00000, -0.00000, -0.00000
26, 0.02959, 0.02959, 0.02959, -0.01533, -0.01533, -0.01533
27, 0.00000, 0.00000, 0.00000, -0.00000, -0.00000, -0.00000
28, 0.02959, 0.02959, -0.02959, -0.01533, -0.01533, 0.01533
29, 0.00000, 0.00000, 0.00210, -0.00000, -0.00000, 0.00000
30, -0.01627, 0.01627, 0.00000, -0.03656, 0.03656, 0.00000
31, 0.00000, 0.00586, -0.00586, -0.00000, 0.01417, -0.01417
32, 0.02959, -0.02959, 0.02959, -0.01533, 0.01533, -0.01533
33, 0.00000, 0.00210, 0.00000, -0.00000, 0.00000, -0.00000
34, 0.02959, -0.02959, -0.02959, -0.01533, 0.01533, 0.01533
35, 0.00000, 0.00586, 0.00586, -0.00000, 0.01417, 0.01417
36, -0.10109, 0.00000, 0.00000, -0.12209, 0.00000, 0.00000
37, 0.13232, -0.13232, -0.13232, 0.14331, -0.14331, -0.14331
38, 0.01627, 0.00000, -0.01627, 0.03656, 0.00000, -0.03656
39, 0.00586, -0.00586, 0.00000, 0.01417, -0.01417, -0.00000
40, 0.01627, 0.00000, 0.01627, 0.03656, 0.00000, 0.03656
41, 0.13232, -0.13232, 0.13232, 0.14331, -0.14331, 0.14331
42, 0.01627, -0.01627, 0.00000, 0.03656, -0.03656, 0.00000
43, 0.00586, 0.00000, -0.00586, 0.01417, -0.00000, -0.01417
44, -0.02959, 0.02959, 0.02959, 0.01533, -0.01533, -0.01533
45, 0.00210, 0.00000, 0.00000, 0.00000, -0.00000, -0.00000
46, -0.02959, 0.02959, -0.02959, 0.01533, -0.01533, 0.01533
47, 0.00586, 0.00000, 0.00586, 0.01417, -0.00000, 0.01417
48, 0.01627, 0.01627, 0.00000, 0.03656, 0.03656, 0.00000
49, 0.13232, 0.13232, -0.13232, 0.14331, 0.14331, -0.14331
50, -0.02959, -0.02959, 0.02959, 0.01533, 0.01533, -0.01533
51, 0.00586, 0.00586, 0.00000, 0.01417, 0.01417, -0.00000
52, -0.02959, -0.02959, -0.02959, 0.01533, 0.01533, 0.01533
53, 0.13232, 0.13232, 0.13232, 0.14331, 0.14331, 0.14331

>> CalcForce(para,2.8553,12,1)
1, -0.12101, -0.12101, -0.12101, -0.11610, -0.11610, -0.11610
2, 0.00000, 0.00000, 0.12340, -0.00001, -0.00001, 0.13232
3, -0.01059, -0.01059, 0.00000, -0.01396, -0.01396, 0.00000
4, 0.00000, 0.00000, -0.12340, -0.00001, -0.00001, -0.13231
5, -0.12101, -0.12101, 0.12101, -0.11610, -0.11610, 0.11611
6, 0.00000, 0.12340, 0.00000, -0.00001, 0.13232, -0.00001
7, -0.01059, 0.00000, -0.01059, -0.01396, 0.00000, -0.01396
8, 0.00000, -0.01398, -0.01398, -0.00001, -0.03611, -0.03611
9, -0.00361, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01398, 0.01398, -0.00001, -0.03611, 0.03611
11, -0.01059, 0.00000, 0.01059, -0.01396, 0.00000, 0.01396
12, 0.00000, -0.12340, 0.00000, -0.00001, -0.13231, -0.00001
13, -0.12101, 0.12101, -0.12101, -0.11610, 0.11611, -0.11610
14, 0.00000, 0.01398, -0.01398, -0.00001, 0.03611, -0.03611
15, -0.01059, 0.01059, 0.00000, -0.01396, 0.01396, 0.00000
16, 0.00000, 0.01398, 0.01398, -0.00001, 0.03611, 0.03611
17, -0.12101, 0.12101, 0.12101, -0.11610, 0.11611, 0.11611
18, 0.12340, 0.00000, 0.00000, 0.13232, -0.00001, -0.00001
19, 0.00000, -0.01059, -0.01059, 0.00000, -0.01396, -0.01396
20, -0.01398, 0.00000, -0.01398, -0.03611, -0.00001, -0.03611
21, 0.00000, -0.00361, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01398, 0.00000, 0.01398, -0.03611, -0.00001, 0.03611
23, 0.00000, -0.01059, 0.01059, 0.00000, -0.01396, 0.01396
24, -0.01398, -0.01398, 0.00000, -0.03611, -0.03611, -0.00001
25, 0.00000, 0.00000, -0.00361, 0.00000, 0.00000, -0.00000
26, 0.02708, 0.02708, 0.02708, -0.01516, -0.01516, -0.01516
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.02708, 0.02708, -0.02708, -0.01516, -0.01516, 0.01517
29, 0.00000, 0.00000, 0.00361, 0.00000, 0.00000, 0.00001
30, -0.01398, 0.01398, 0.00000, -0.03611, 0.03611, -0.00001
31, 0.00000, 0.01059, -0.01059, 0.00000, 0.01396, -0.01396
32, 0.02708, -0.02708, 0.02708, -0.01516, 0.01517, -0.01516
33, 0.00000, 0.00361, 0.00000, 0.00000, 0.00001, 0.00000
34, 0.02708, -0.02708, -0.02708, -0.01516, 0.01517, 0.01517
35, 0.00000, 0.01059, 0.01059, 0.00000, 0.01396, 0.01396
36, -0.12340, 0.00000, 0.00000, -0.13231, -0.00001, -0.00001
37, 0.12101, -0.12101, -0.12101, 0.11611, -0.11610, -0.11610
38, 0.01398, 0.00000, -0.01398, 0.03611, -0.00001, -0.03611
39, 0.01059, -0.01059, 0.00000, 0.01396, -0.01396, 0.00000
40, 0.01398, 0.00000, 0.01398, 0.03611, -0.00001, 0.03611
41, 0.12101, -0.12101, 0.12101, 0.11611, -0.11610, 0.11611
42, 0.01398, -0.01398, 0.00000, 0.03611, -0.03611, -0.00001
43, 0.01059, 0.00000, -0.01059, 0.01396, 0.00000, -0.01396
44, -0.02708, 0.02708, 0.02708, 0.01517, -0.01516, -0.01516
45, 0.00361, 0.00000, 0.00000, 0.00001, 0.00000, 0.00000
46, -0.02708, 0.02708, -0.02708, 0.01517, -0.01516, 0.01517
47, 0.01059, 0.00000, 0.01059, 0.01396, 0.00000, 0.01396
48, 0.01398, 0.01398, 0.00000, 0.03611, 0.03611, -0.00001
49, 0.12101, 0.12101, -0.12101, 0.11611, 0.11611, -0.11610
50, -0.02708, -0.02708, 0.02708, 0.01517, 0.01517, -0.01516
51, 0.01059, 0.01059, 0.00000, 0.01396, 0.01396, 0.00000
52, -0.02708, -0.02708, -0.02708, 0.01517, 0.01517, 0.01517
53, 0.12101, 0.12101, 0.12101, 0.11611, 0.11611, 0.11611
>> CalcForce(para,2.8665,13,1)
1, -0.11336, -0.11336, -0.11336, -0.09713, -0.09713, -0.09713
2, 0.00000, 0.00000, 0.12556, -0.00000, -0.00000, 0.13928
3, -0.00937, -0.00937, 0.00000, -0.01380, -0.01380, 0.00000
4, 0.00000, 0.00000, -0.12556, -0.00000, -0.00000, -0.13928
5, -0.11336, -0.11336, 0.11336, -0.09713, -0.09713, 0.09713
6, 0.00000, 0.12556, 0.00000, -0.00000, 0.13928, -0.00000
7, -0.00937, 0.00000, -0.00937, -0.01380, 0.00000, -0.01380
8, 0.00000, -0.01110, -0.01110, -0.00000, -0.03578, -0.03578
9, -0.01130, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01110, 0.01110, -0.00000, -0.03578, 0.03578
11, -0.00937, 0.00000, 0.00937, -0.01380, 0.00000, 0.01380
12, 0.00000, -0.12556, 0.00000, -0.00000, -0.13928, -0.00000
13, -0.11336, 0.11336, -0.11336, -0.09713, 0.09713, -0.09713
14, 0.00000, 0.01110, -0.01110, -0.00000, 0.03578, -0.03578
15, -0.00937, 0.00937, 0.00000, -0.01380, 0.01380, 0.00000
16, 0.00000, 0.01110, 0.01110, -0.00000, 0.03578, 0.03578
17, -0.11336, 0.11336, 0.11336, -0.09713, 0.09713, 0.09713
18, 0.12556, 0.00000, 0.00000, 0.13928, -0.00000, -0.00000
19, 0.00000, -0.00937, -0.00937, 0.00000, -0.01380, -0.01380
20, -0.01110, 0.00000, -0.01110, -0.03578, -0.00000, -0.03578
21, 0.00000, -0.01130, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01110, 0.00000, 0.01110, -0.03578, -0.00000, 0.03578
23, 0.00000, -0.00937, 0.00937, 0.00000, -0.01380, 0.01380
24, -0.01110, -0.01110, 0.00000, -0.03578, -0.03578, -0.00000
25, 0.00000, 0.00000, -0.01130, 0.00000, 0.00000, -0.00000
26, 0.01773, 0.01773, 0.01773, -0.01504, -0.01504, -0.01504
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.01773, 0.01773, -0.01773, -0.01504, -0.01504, 0.01504
29, 0.00000, 0.00000, 0.01130, 0.00000, 0.00000, 0.00000
30, -0.01110, 0.01110, 0.00000, -0.03578, 0.03578, -0.00000
31, 0.00000, 0.00937, -0.00937, 0.00000, 0.01380, -0.01380
32, 0.01773, -0.01773, 0.01773, -0.01504, 0.01504, -0.01504
33, 0.00000, 0.01130, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.01773, -0.01773, -0.01773, -0.01504, 0.01504, 0.01504
35, 0.00000, 0.00937, 0.00937, 0.00000, 0.01380, 0.01380
36, -0.12556, 0.00000, 0.00000, -0.13928, -0.00000, -0.00000
37, 0.11336, -0.11336, -0.11336, 0.09713, -0.09713, -0.09713
38, 0.01110, 0.00000, -0.01110, 0.03578, -0.00000, -0.03578
39, 0.00937, -0.00937, 0.00000, 0.01380, -0.01380, 0.00000
40, 0.01110, 0.00000, 0.01110, 0.03578, -0.00000, 0.03578
41, 0.11336, -0.11336, 0.11336, 0.09713, -0.09713, 0.09713
42, 0.01110, -0.01110, 0.00000, 0.03578, -0.03578, -0.00000
43, 0.00937, 0.00000, -0.00937, 0.01380, 0.00000, -0.01380
44, -0.01773, 0.01773, 0.01773, 0.01504, -0.01504, -0.01504
45, 0.01130, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.01773, 0.01773, -0.01773, 0.01504, -0.01504, 0.01504
47, 0.00937, 0.00000, 0.00937, 0.01380, 0.00000, 0.01380
48, 0.01110, 0.01110, 0.00000, 0.03578, 0.03578, -0.00000
49, 0.11336, 0.11336, -0.11336, 0.09713, 0.09713, -0.09713
50, -0.01773, -0.01773, 0.01773, 0.01504, 0.01504, -0.01504
51, 0.00937, 0.00937, 0.00000, 0.01380, 0.01380, 0.00000
52, -0.01773, -0.01773, -0.01773, 0.01504, 0.01504, 0.01504
53, 0.11336, 0.11336, 0.11336, 0.09713, 0.09713, 0.09713

PARA6(6 para): 3.7877 8.5591 4.7566 2.0454 1.7929 1.7630

>> CalcForce(para,2.82,10,1)
1, -0.20033, -0.20033, -0.20033, -0.18927, -0.18927, -0.18927
2, 0.00000, 0.00000, 0.13389, -0.00000, -0.00000, 0.16864
3, -0.01955, -0.01955, 0.00000, -0.00810, -0.00810, 0.00000
4, 0.00000, 0.00000, -0.13389, -0.00000, -0.00000, -0.16864
5, -0.20033, -0.20033, 0.20033, -0.18927, -0.18927, 0.18927
6, 0.00000, 0.13389, 0.00000, -0.00000, 0.16864, -0.00000
7, -0.01955, 0.00000, -0.01955, -0.00810, 0.00000, -0.00810
8, 0.00000, -0.01429, -0.01429, -0.00000, -0.01738, -0.01738
9, -0.02339, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01429, 0.01429, -0.00000, -0.01738, 0.01738
11, -0.01955, 0.00000, 0.01955, -0.00810, 0.00000, 0.00811
12, 0.00000, -0.13389, 0.00000, -0.00000, -0.16864, -0.00000
13, -0.20033, 0.20033, -0.20033, -0.18927, 0.18927, -0.18927
14, 0.00000, 0.01429, -0.01429, -0.00000, 0.01738, -0.01738
15, -0.01955, 0.01955, 0.00000, -0.00810, 0.00811, 0.00000
16, 0.00000, 0.01429, 0.01429, -0.00000, 0.01738, 0.01738
17, -0.20033, 0.20033, 0.20033, -0.18927, 0.18927, 0.18927
18, 0.13389, 0.00000, 0.00000, 0.16864, -0.00000, -0.00000
19, 0.00000, -0.01955, -0.01955, 0.00000, -0.00810, -0.00810
20, -0.01429, 0.00000, -0.01429, -0.01738, -0.00000, -0.01738
21, 0.00000, -0.02339, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01429, 0.00000, 0.01429, -0.01738, -0.00000, 0.01738
23, 0.00000, -0.01955, 0.01955, 0.00000, -0.00810, 0.00811
24, -0.01429, -0.01429, 0.00000, -0.01738, -0.01738, -0.00000
25, 0.00000, 0.00000, -0.02339, 0.00000, 0.00000, -0.00000
26, 0.03201, 0.03201, 0.03201, -0.00643, -0.00643, -0.00643
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.03201, 0.03201, -0.03201, -0.00643, -0.00643, 0.00643
29, 0.00000, 0.00000, 0.02339, 0.00000, 0.00000, 0.00000
30, -0.01429, 0.01429, 0.00000, -0.01738, 0.01738, -0.00000
31, 0.00000, 0.01955, -0.01955, 0.00000, 0.00811, -0.00810
32, 0.03201, -0.03201, 0.03201, -0.00643, 0.00643, -0.00643
33, 0.00000, 0.02339, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.03201, -0.03201, -0.03201, -0.00643, 0.00643, 0.00643
35, 0.00000, 0.01955, 0.01955, 0.00000, 0.00811, 0.00811
36, -0.13389, 0.00000, 0.00000, -0.16864, -0.00000, -0.00000
37, 0.20033, -0.20033, -0.20033, 0.18927, -0.18927, -0.18927
38, 0.01429, 0.00000, -0.01429, 0.01738, -0.00000, -0.01738
39, 0.01955, -0.01955, 0.00000, 0.00811, -0.00810, 0.00000
40, 0.01429, 0.00000, 0.01429, 0.01738, -0.00000, 0.01738
41, 0.20033, -0.20033, 0.20033, 0.18927, -0.18927, 0.18927
42, 0.01429, -0.01429, 0.00000, 0.01738, -0.01738, -0.00000
43, 0.01955, 0.00000, -0.01955, 0.00811, 0.00000, -0.00810
44, -0.03201, 0.03201, 0.03201, 0.00643, -0.00643, -0.00643
45, 0.02339, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.03201, 0.03201, -0.03201, 0.00643, -0.00643, 0.00643
47, 0.01955, 0.00000, 0.01955, 0.00811, 0.00000, 0.00811
48, 0.01429, 0.01429, 0.00000, 0.01738, 0.01738, -0.00000
49, 0.20033, 0.20033, -0.20033, 0.18927, 0.18927, -0.18927
50, -0.03201, -0.03201, 0.03201, 0.00643, 0.00643, -0.00643
51, 0.01955, 0.01955, 0.00000, 0.00811, 0.00811, 0.00000
52, -0.03201, -0.03201, -0.03201, 0.00643, 0.00643, 0.00643
53, 0.20033, 0.20033, 0.20033, 0.18927, 0.18927, 0.18927
>> CalcForce(para,2.8665,13,1)
1, -0.15115, -0.15115, -0.15115, -0.11310, -0.11310, -0.11310
2, 0.00000, 0.00000, 0.16742, -0.00000, -0.00000, 0.18628
3, -0.01250, -0.01250, 0.00000, -0.00789, -0.00789, 0.00000
4, 0.00000, 0.00000, -0.16742, -0.00000, -0.00000, -0.18628
5, -0.15115, -0.15115, 0.15115, -0.11310, -0.11310, 0.11310
6, 0.00000, 0.16742, 0.00000, -0.00000, 0.18628, -0.00000
7, -0.01250, 0.00000, -0.01250, -0.00789, 0.00000, -0.00789
8, 0.00000, -0.01480, -0.01480, -0.00000, -0.01698, -0.01698
9, -0.01507, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01480, 0.01480, -0.00000, -0.01698, 0.01698
11, -0.01250, 0.00000, 0.01250, -0.00789, 0.00000, 0.00789
12, 0.00000, -0.16742, 0.00000, -0.00000, -0.18628, -0.00000
13, -0.15115, 0.15115, -0.15115, -0.11310, 0.11310, -0.11310
14, 0.00000, 0.01480, -0.01480, -0.00000, 0.01698, -0.01698
15, -0.01250, 0.01250, 0.00000, -0.00789, 0.00789, 0.00000
16, 0.00000, 0.01480, 0.01480, -0.00000, 0.01698, 0.01698
17, -0.15115, 0.15115, 0.15115, -0.11310, 0.11310, 0.11310
18, 0.16742, 0.00000, 0.00000, 0.18628, -0.00000, -0.00000
19, 0.00000, -0.01250, -0.01250, 0.00000, -0.00789, -0.00789
20, -0.01480, 0.00000, -0.01480, -0.01698, -0.00000, -0.01698
21, 0.00000, -0.01507, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01480, 0.00000, 0.01480, -0.01698, -0.00000, 0.01698
23, 0.00000, -0.01250, 0.01250, 0.00000, -0.00789, 0.00789
24, -0.01480, -0.01480, 0.00000, -0.01698, -0.01698, -0.00000
25, 0.00000, 0.00000, -0.01507, 0.00000, 0.00000, -0.00000
26, 0.02364, 0.02364, 0.02364, -0.00630, -0.00630, -0.00630
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.02364, 0.02364, -0.02364, -0.00630, -0.00630, 0.00630
29, 0.00000, 0.00000, 0.01507, 0.00000, 0.00000, 0.00000
30, -0.01480, 0.01480, 0.00000, -0.01698, 0.01698, -0.00000
31, 0.00000, 0.01250, -0.01250, 0.00000, 0.00789, -0.00789
32, 0.02364, -0.02364, 0.02364, -0.00630, 0.00630, -0.00630
33, 0.00000, 0.01507, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.02364, -0.02364, -0.02364, -0.00630, 0.00630, 0.00630
35, 0.00000, 0.01250, 0.01250, 0.00000, 0.00789, 0.00789
36, -0.16742, 0.00000, 0.00000, -0.18628, -0.00000, -0.00000
37, 0.15115, -0.15115, -0.15115, 0.11310, -0.11310, -0.11310
38, 0.01480, 0.00000, -0.01480, 0.01698, -0.00000, -0.01698
39, 0.01250, -0.01250, 0.00000, 0.00789, -0.00789, 0.00000
40, 0.01480, 0.00000, 0.01480, 0.01698, -0.00000, 0.01698
41, 0.15115, -0.15115, 0.15115, 0.11310, -0.11310, 0.11310
42, 0.01480, -0.01480, 0.00000, 0.01698, -0.01698, -0.00000
43, 0.01250, 0.00000, -0.01250, 0.00789, 0.00000, -0.00789
44, -0.02364, 0.02364, 0.02364, 0.00630, -0.00630, -0.00630
45, 0.01507, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.02364, 0.02364, -0.02364, 0.00630, -0.00630, 0.00630
47, 0.01250, 0.00000, 0.01250, 0.00789, 0.00000, 0.00789
48, 0.01480, 0.01480, 0.00000, 0.01698, 0.01698, -0.00000
49, 0.15115, 0.15115, -0.15115, 0.11310, 0.11310, -0.11310
50, -0.02364, -0.02364, 0.02364, 0.00630, 0.00630, -0.00630
51, 0.01250, 0.01250, 0.00000, 0.00789, 0.00789, 0.00000
52, -0.02364, -0.02364, -0.02364, 0.00630, 0.00630, 0.00630
53, 0.15115, 0.15115, 0.15115, 0.11310, 0.11310, 0.11310

PARA5(6 para): 3.7301 8.3466 4.7830 2.1489 1.7884 1.7894

>> CalcForce(para,2.82,10,1)
1, -0.20033, -0.20033, -0.20033, -0.18205, -0.18205, -0.18205
2, 0.00000, 0.00000, 0.13389, -0.00000, -0.00000, 0.16910
3, -0.01955, -0.01955, 0.00000, -0.00836, -0.00836, 0.00000
4, 0.00000, 0.00000, -0.13389, -0.00000, -0.00000, -0.16910
5, -0.20033, -0.20033, 0.20033, -0.18205, -0.18205, 0.18205
6, 0.00000, 0.13389, 0.00000, -0.00000, 0.16910, -0.00000
7, -0.01955, 0.00000, -0.01955, -0.00836, 0.00000, -0.00836
8, 0.00000, -0.01429, -0.01429, -0.00000, -0.01808, -0.01808
9, -0.02339, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01429, 0.01429, -0.00000, -0.01808, 0.01808
11, -0.01955, 0.00000, 0.01955, -0.00836, 0.00000, 0.00836
12, 0.00000, -0.13389, 0.00000, -0.00000, -0.16910, -0.00000
13, -0.20033, 0.20033, -0.20033, -0.18205, 0.18205, -0.18205
14, 0.00000, 0.01429, -0.01429, -0.00000, 0.01808, -0.01808
15, -0.01955, 0.01955, 0.00000, -0.00836, 0.00836, 0.00000
16, 0.00000, 0.01429, 0.01429, -0.00000, 0.01808, 0.01808
17, -0.20033, 0.20033, 0.20033, -0.18205, 0.18205, 0.18205
18, 0.13389, 0.00000, 0.00000, 0.16910, -0.00000, -0.00000
19, 0.00000, -0.01955, -0.01955, 0.00000, -0.00836, -0.00836
20, -0.01429, 0.00000, -0.01429, -0.01808, -0.00000, -0.01808
21, 0.00000, -0.02339, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01429, 0.00000, 0.01429, -0.01808, -0.00000, 0.01808
23, 0.00000, -0.01955, 0.01955, 0.00000, -0.00836, 0.00836
24, -0.01429, -0.01429, 0.00000, -0.01808, -0.01808, -0.00000
25, 0.00000, 0.00000, -0.02339, 0.00000, 0.00000, -0.00000
26, 0.03201, 0.03201, 0.03201, -0.00674, -0.00674, -0.00674
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.03201, 0.03201, -0.03201, -0.00674, -0.00674, 0.00674
29, 0.00000, 0.00000, 0.02339, 0.00000, 0.00000, 0.00000
30, -0.01429, 0.01429, 0.00000, -0.01808, 0.01808, -0.00000
31, 0.00000, 0.01955, -0.01955, 0.00000, 0.00836, -0.00836
32, 0.03201, -0.03201, 0.03201, -0.00674, 0.00674, -0.00674
33, 0.00000, 0.02339, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.03201, -0.03201, -0.03201, -0.00674, 0.00674, 0.00674
35, 0.00000, 0.01955, 0.01955, 0.00000, 0.00836, 0.00836
36, -0.13389, 0.00000, 0.00000, -0.16910, -0.00000, -0.00000
37, 0.20033, -0.20033, -0.20033, 0.18205, -0.18205, -0.18205
38, 0.01429, 0.00000, -0.01429, 0.01808, -0.00000, -0.01808
39, 0.01955, -0.01955, 0.00000, 0.00836, -0.00836, 0.00000
40, 0.01429, 0.00000, 0.01429, 0.01808, -0.00000, 0.01808
41, 0.20033, -0.20033, 0.20033, 0.18205, -0.18205, 0.18205
42, 0.01429, -0.01429, 0.00000, 0.01808, -0.01808, -0.00000
43, 0.01955, 0.00000, -0.01955, 0.00836, 0.00000, -0.00836
44, -0.03201, 0.03201, 0.03201, 0.00674, -0.00674, -0.00674
45, 0.02339, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.03201, 0.03201, -0.03201, 0.00674, -0.00674, 0.00674
47, 0.01955, 0.00000, 0.01955, 0.00836, 0.00000, 0.00836
48, 0.01429, 0.01429, 0.00000, 0.01808, 0.01808, -0.00000
49, 0.20033, 0.20033, -0.20033, 0.18205, 0.18205, -0.18205
50, -0.03201, -0.03201, 0.03201, 0.00674, 0.00674, -0.00674
51, 0.01955, 0.01955, 0.00000, 0.00836, 0.00836, 0.00000
52, -0.03201, -0.03201, -0.03201, 0.00674, 0.00674, 0.00674
53, 0.20033, 0.20033, 0.20033, 0.18205, 0.18205, 0.18205
>> CalcForce(para,2.8665,13,1)
1, -0.15115, -0.15115, -0.15115, -0.10720, -0.10720, -0.10720
2, 0.00000, 0.00000, 0.16742, -0.00000, -0.00000, 0.18702
3, -0.01250, -0.01250, 0.00000, -0.00813, -0.00813, 0.00000
4, 0.00000, 0.00000, -0.16742, -0.00000, -0.00000, -0.18702
5, -0.15115, -0.15115, 0.15115, -0.10720, -0.10720, 0.10720
6, 0.00000, 0.16742, 0.00000, -0.00000, 0.18702, -0.00000
7, -0.01250, 0.00000, -0.01250, -0.00813, 0.00000, -0.00813
8, 0.00000, -0.01480, -0.01480, -0.00000, -0.01764, -0.01764
9, -0.01507, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01480, 0.01480, -0.00000, -0.01764, 0.01764
11, -0.01250, 0.00000, 0.01250, -0.00813, 0.00000, 0.00813
12, 0.00000, -0.16742, 0.00000, -0.00000, -0.18702, -0.00000
13, -0.15115, 0.15115, -0.15115, -0.10720, 0.10720, -0.10720
14, 0.00000, 0.01480, -0.01480, -0.00000, 0.01764, -0.01764
15, -0.01250, 0.01250, 0.00000, -0.00813, 0.00813, 0.00000
16, 0.00000, 0.01480, 0.01480, -0.00000, 0.01764, 0.01764
17, -0.15115, 0.15115, 0.15115, -0.10720, 0.10720, 0.10720
18, 0.16742, 0.00000, 0.00000, 0.18702, -0.00000, -0.00000
19, 0.00000, -0.01250, -0.01250, 0.00000, -0.00813, -0.00813
20, -0.01480, 0.00000, -0.01480, -0.01764, -0.00000, -0.01764
21, 0.00000, -0.01507, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01480, 0.00000, 0.01480, -0.01764, -0.00000, 0.01764
23, 0.00000, -0.01250, 0.01250, 0.00000, -0.00813, 0.00813
24, -0.01480, -0.01480, 0.00000, -0.01764, -0.01764, -0.00000
25, 0.00000, 0.00000, -0.01507, 0.00000, 0.00000, -0.00000
26, 0.02364, 0.02364, 0.02364, -0.00660, -0.00660, -0.00660
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.02364, 0.02364, -0.02364, -0.00660, -0.00660, 0.00660
29, 0.00000, 0.00000, 0.01507, 0.00000, 0.00000, 0.00000
30, -0.01480, 0.01480, 0.00000, -0.01764, 0.01764, -0.00000
31, 0.00000, 0.01250, -0.01250, 0.00000, 0.00813, -0.00813
32, 0.02364, -0.02364, 0.02364, -0.00660, 0.00660, -0.00660
33, 0.00000, 0.01507, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.02364, -0.02364, -0.02364, -0.00660, 0.00660, 0.00660
35, 0.00000, 0.01250, 0.01250, 0.00000, 0.00813, 0.00813
36, -0.16742, 0.00000, 0.00000, -0.18702, -0.00000, -0.00000
37, 0.15115, -0.15115, -0.15115, 0.10720, -0.10720, -0.10720
38, 0.01480, 0.00000, -0.01480, 0.01764, -0.00000, -0.01764
39, 0.01250, -0.01250, 0.00000, 0.00813, -0.00813, 0.00000
40, 0.01480, 0.00000, 0.01480, 0.01764, -0.00000, 0.01764
41, 0.15115, -0.15115, 0.15115, 0.10720, -0.10720, 0.10720
42, 0.01480, -0.01480, 0.00000, 0.01764, -0.01764, -0.00000
43, 0.01250, 0.00000, -0.01250, 0.00813, 0.00000, -0.00813
44, -0.02364, 0.02364, 0.02364, 0.00660, -0.00660, -0.00660
45, 0.01507, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.02364, 0.02364, -0.02364, 0.00660, -0.00660, 0.00660
47, 0.01250, 0.00000, 0.01250, 0.00813, 0.00000, 0.00813
48, 0.01480, 0.01480, 0.00000, 0.01764, 0.01764, -0.00000
49, 0.15115, 0.15115, -0.15115, 0.10720, 0.10720, -0.10720
50, -0.02364, -0.02364, 0.02364, 0.00660, 0.00660, -0.00660
51, 0.01250, 0.01250, 0.00000, 0.00813, 0.00813, 0.00000
52, -0.02364, -0.02364, -0.02364, 0.00660, 0.00660, 0.00660
53, 0.15115, 0.15115, 0.15115, 0.10720, 0.10720, 0.10720

PARA4(6 para): 3.6149 7.9288 4.7911 2.1121 1.7760 1.8761

>> CalcForce(para,2.82,10,1)
1, -0.20033, -0.20033, -0.20033, -0.18146, -0.18146, -0.18146
2, 0.00000, 0.00000, 0.13389, -0.00000, -0.00000, 0.16944
3, -0.01955, -0.01955, 0.00000, -0.00857, -0.00857, 0.00000
4, 0.00000, 0.00000, -0.13389, -0.00000, -0.00000, -0.16944
5, -0.20033, -0.20033, 0.20033, -0.18146, -0.18146, 0.18146
6, 0.00000, 0.13389, 0.00000, -0.00000, 0.16944, -0.00000
7, -0.01955, 0.00000, -0.01955, -0.00857, 0.00000, -0.00857
8, 0.00000, -0.01429, -0.01429, -0.00000, -0.01826, -0.01826
9, -0.02339, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01429, 0.01429, -0.00000, -0.01826, 0.01826
11, -0.01955, 0.00000, 0.01955, -0.00857, 0.00000, 0.00857
12, 0.00000, -0.13389, 0.00000, -0.00000, -0.16944, -0.00000
13, -0.20033, 0.20033, -0.20033, -0.18146, 0.18146, -0.18146
14, 0.00000, 0.01429, -0.01429, -0.00000, 0.01826, -0.01826
15, -0.01955, 0.01955, 0.00000, -0.00857, 0.00857, 0.00000
16, 0.00000, 0.01429, 0.01429, -0.00000, 0.01826, 0.01826
17, -0.20033, 0.20033, 0.20033, -0.18146, 0.18146, 0.18146
18, 0.13389, 0.00000, 0.00000, 0.16944, -0.00000, -0.00000
19, 0.00000, -0.01955, -0.01955, 0.00000, -0.00857, -0.00857
20, -0.01429, 0.00000, -0.01429, -0.01826, -0.00000, -0.01826
21, 0.00000, -0.02339, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01429, 0.00000, 0.01429, -0.01826, -0.00000, 0.01826
23, 0.00000, -0.01955, 0.01955, 0.00000, -0.00857, 0.00857
24, -0.01429, -0.01429, 0.00000, -0.01826, -0.01826, -0.00000
25, 0.00000, 0.00000, -0.02339, 0.00000, 0.00000, -0.00000
26, 0.03201, 0.03201, 0.03201, -0.00670, -0.00670, -0.00670
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.03201, 0.03201, -0.03201, -0.00670, -0.00670, 0.00670
29, 0.00000, 0.00000, 0.02339, 0.00000, 0.00000, 0.00000
30, -0.01429, 0.01429, 0.00000, -0.01826, 0.01826, -0.00000
31, 0.00000, 0.01955, -0.01955, 0.00000, 0.00857, -0.00857
32, 0.03201, -0.03201, 0.03201, -0.00670, 0.00670, -0.00670
33, 0.00000, 0.02339, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.03201, -0.03201, -0.03201, -0.00670, 0.00670, 0.00670
35, 0.00000, 0.01955, 0.01955, 0.00000, 0.00857, 0.00857
36, -0.13389, 0.00000, 0.00000, -0.16944, -0.00000, -0.00000
37, 0.20033, -0.20033, -0.20033, 0.18146, -0.18146, -0.18146
38, 0.01429, 0.00000, -0.01429, 0.01826, -0.00000, -0.01826
39, 0.01955, -0.01955, 0.00000, 0.00857, -0.00857, 0.00000
40, 0.01429, 0.00000, 0.01429, 0.01826, -0.00000, 0.01826
41, 0.20033, -0.20033, 0.20033, 0.18146, -0.18146, 0.18146
42, 0.01429, -0.01429, 0.00000, 0.01826, -0.01826, -0.00000
43, 0.01955, 0.00000, -0.01955, 0.00857, 0.00000, -0.00857
44, -0.03201, 0.03201, 0.03201, 0.00670, -0.00670, -0.00670
45, 0.02339, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.03201, 0.03201, -0.03201, 0.00670, -0.00670, 0.00670
47, 0.01955, 0.00000, 0.01955, 0.00857, 0.00000, 0.00857
48, 0.01429, 0.01429, 0.00000, 0.01826, 0.01826, -0.00000
49, 0.20033, 0.20033, -0.20033, 0.18146, 0.18146, -0.18146
50, -0.03201, -0.03201, 0.03201, 0.00670, 0.00670, -0.00670
51, 0.01955, 0.01955, 0.00000, 0.00857, 0.00857, 0.00000
52, -0.03201, -0.03201, -0.03201, 0.00670, 0.00670, 0.00670
53, 0.20033, 0.20033, 0.20033, 0.18146, 0.18146, 0.18146
>> CalcForce(para,2.8665,13,1)
1, -0.15115, -0.15115, -0.15115, -0.10935, -0.10935, -0.10935
2, 0.00000, 0.00000, 0.16742, -0.00000, -0.00000, 0.18869
3, -0.01250, -0.01250, 0.00000, -0.00835, -0.00835, 0.00000
4, 0.00000, 0.00000, -0.16742, -0.00000, -0.00000, -0.18869
5, -0.15115, -0.15115, 0.15115, -0.10935, -0.10935, 0.10935
6, 0.00000, 0.16742, 0.00000, -0.00000, 0.18869, -0.00000
7, -0.01250, 0.00000, -0.01250, -0.00835, 0.00000, -0.00835
8, 0.00000, -0.01480, -0.01480, -0.00000, -0.01785, -0.01785
9, -0.01507, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01480, 0.01480, -0.00000, -0.01785, 0.01785
11, -0.01250, 0.00000, 0.01250, -0.00835, 0.00000, 0.00835
12, 0.00000, -0.16742, 0.00000, -0.00000, -0.18869, -0.00000
13, -0.15115, 0.15115, -0.15115, -0.10935, 0.10935, -0.10935
14, 0.00000, 0.01480, -0.01480, -0.00000, 0.01785, -0.01785
15, -0.01250, 0.01250, 0.00000, -0.00835, 0.00835, 0.00000
16, 0.00000, 0.01480, 0.01480, -0.00000, 0.01785, 0.01785
17, -0.15115, 0.15115, 0.15115, -0.10935, 0.10935, 0.10935
18, 0.16742, 0.00000, 0.00000, 0.18869, -0.00000, -0.00000
19, 0.00000, -0.01250, -0.01250, 0.00000, -0.00835, -0.00835
20, -0.01480, 0.00000, -0.01480, -0.01785, -0.00000, -0.01785
21, 0.00000, -0.01507, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01480, 0.00000, 0.01480, -0.01785, -0.00000, 0.01785
23, 0.00000, -0.01250, 0.01250, 0.00000, -0.00835, 0.00835
24, -0.01480, -0.01480, 0.00000, -0.01785, -0.01785, -0.00000
25, 0.00000, 0.00000, -0.01507, 0.00000, 0.00000, -0.00000
26, 0.02364, 0.02364, 0.02364, -0.00658, -0.00658, -0.00658
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.02364, 0.02364, -0.02364, -0.00658, -0.00658, 0.00658
29, 0.00000, 0.00000, 0.01507, 0.00000, 0.00000, 0.00000
30, -0.01480, 0.01480, 0.00000, -0.01785, 0.01785, -0.00000
31, 0.00000, 0.01250, -0.01250, 0.00000, 0.00835, -0.00835
32, 0.02364, -0.02364, 0.02364, -0.00658, 0.00658, -0.00658
33, 0.00000, 0.01507, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.02364, -0.02364, -0.02364, -0.00658, 0.00658, 0.00658
35, 0.00000, 0.01250, 0.01250, 0.00000, 0.00835, 0.00835
36, -0.16742, 0.00000, 0.00000, -0.18869, -0.00000, -0.00000
37, 0.15115, -0.15115, -0.15115, 0.10935, -0.10935, -0.10935
38, 0.01480, 0.00000, -0.01480, 0.01785, -0.00000, -0.01785
39, 0.01250, -0.01250, 0.00000, 0.00835, -0.00835, 0.00000
40, 0.01480, 0.00000, 0.01480, 0.01785, -0.00000, 0.01785
41, 0.15115, -0.15115, 0.15115, 0.10935, -0.10935, 0.10935
42, 0.01480, -0.01480, 0.00000, 0.01785, -0.01785, -0.00000
43, 0.01250, 0.00000, -0.01250, 0.00835, 0.00000, -0.00835
44, -0.02364, 0.02364, 0.02364, 0.00658, -0.00658, -0.00658
45, 0.01507, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.02364, 0.02364, -0.02364, 0.00658, -0.00658, 0.00658
47, 0.01250, 0.00000, 0.01250, 0.00835, 0.00000, 0.00835
48, 0.01480, 0.01480, 0.00000, 0.01785, 0.01785, -0.00000
49, 0.15115, 0.15115, -0.15115, 0.10935, 0.10935, -0.10935
50, -0.02364, -0.02364, 0.02364, 0.00658, 0.00658, -0.00658
51, 0.01250, 0.01250, 0.00000, 0.00835, 0.00835, 0.00000
52, -0.02364, -0.02364, -0.02364, 0.00658, 0.00658, 0.00658
53, 0.15115, 0.15115, 0.15115, 0.10935, 0.10935, 0.10935

PARA3(6 para): 3.7705 8.5241 4.7915 2.2370 1.7951 1.8203

>> CalcForce(para,2.82,10,1)
1, -0.20033, -0.20033, -0.20033, -0.17443, -0.17443, -0.17443
2, 0.00000, 0.00000, 0.13389, -0.00000, -0.00000, 0.17984
3, -0.01955, -0.01955, 0.00000, -0.00860, -0.00860, 0.00000
4, 0.00000, 0.00000, -0.13389, -0.00000, -0.00000, -0.17984
5, -0.20033, -0.20033, 0.20033, -0.17443, -0.17443, 0.17444
6, 0.00000, 0.13389, 0.00000, -0.00000, 0.17984, -0.00000
7, -0.01955, 0.00000, -0.01955, -0.00860, 0.00000, -0.00860
8, 0.00000, -0.01429, -0.01429, -0.00000, -0.01870, -0.01870
9, -0.02339, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01429, 0.01429, -0.00000, -0.01870, 0.01870
11, -0.01955, 0.00000, 0.01955, -0.00860, 0.00000, 0.00860
12, 0.00000, -0.13389, 0.00000, -0.00000, -0.17984, -0.00000
13, -0.20033, 0.20033, -0.20033, -0.17443, 0.17444, -0.17443
14, 0.00000, 0.01429, -0.01429, -0.00000, 0.01870, -0.01870
15, -0.01955, 0.01955, 0.00000, -0.00860, 0.00860, 0.00000
16, 0.00000, 0.01429, 0.01429, -0.00000, 0.01870, 0.01870
17, -0.20033, 0.20033, 0.20033, -0.17443, 0.17444, 0.17444
18, 0.13389, 0.00000, 0.00000, 0.17984, -0.00000, -0.00000
19, 0.00000, -0.01955, -0.01955, 0.00000, -0.00860, -0.00860
20, -0.01429, 0.00000, -0.01429, -0.01870, -0.00000, -0.01870
21, 0.00000, -0.02339, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01429, 0.00000, 0.01429, -0.01870, -0.00000, 0.01870
23, 0.00000, -0.01955, 0.01955, 0.00000, -0.00860, 0.00860
24, -0.01429, -0.01429, 0.00000, -0.01870, -0.01870, -0.00000
25, 0.00000, 0.00000, -0.02339, 0.00000, 0.00000, -0.00000
26, 0.03201, 0.03201, 0.03201, -0.00701, -0.00701, -0.00701
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.03201, 0.03201, -0.03201, -0.00701, -0.00701, 0.00701
29, 0.00000, 0.00000, 0.02339, 0.00000, 0.00000, 0.00000
30, -0.01429, 0.01429, 0.00000, -0.01870, 0.01870, -0.00000
31, 0.00000, 0.01955, -0.01955, 0.00000, 0.00860, -0.00860
32, 0.03201, -0.03201, 0.03201, -0.00701, 0.00701, -0.00701
33, 0.00000, 0.02339, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.03201, -0.03201, -0.03201, -0.00701, 0.00701, 0.00701
35, 0.00000, 0.01955, 0.01955, 0.00000, 0.00860, 0.00860
36, -0.13389, 0.00000, 0.00000, -0.17984, -0.00000, -0.00000
37, 0.20033, -0.20033, -0.20033, 0.17444, -0.17443, -0.17443
38, 0.01429, 0.00000, -0.01429, 0.01870, -0.00000, -0.01870
39, 0.01955, -0.01955, 0.00000, 0.00860, -0.00860, 0.00000
40, 0.01429, 0.00000, 0.01429, 0.01870, -0.00000, 0.01870
41, 0.20033, -0.20033, 0.20033, 0.17444, -0.17443, 0.17444
42, 0.01429, -0.01429, 0.00000, 0.01870, -0.01870, -0.00000
43, 0.01955, 0.00000, -0.01955, 0.00860, 0.00000, -0.00860
44, -0.03201, 0.03201, 0.03201, 0.00701, -0.00701, -0.00701
45, 0.02339, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.03201, 0.03201, -0.03201, 0.00701, -0.00701, 0.00701
47, 0.01955, 0.00000, 0.01955, 0.00860, 0.00000, 0.00860
48, 0.01429, 0.01429, 0.00000, 0.01870, 0.01870, -0.00000
49, 0.20033, 0.20033, -0.20033, 0.17444, 0.17444, -0.17443
50, -0.03201, -0.03201, 0.03201, 0.00701, 0.00701, -0.00701
51, 0.01955, 0.01955, 0.00000, 0.00860, 0.00860, 0.00000
52, -0.03201, -0.03201, -0.03201, 0.00701, 0.00701, 0.00701
53, 0.20033, 0.20033, 0.20033, 0.17443, 0.17443, 0.17443
>> CalcForce(para,2.8665,13,1)
1, -0.15115, -0.15115, -0.15115, -0.09815, -0.09815, -0.09815
2, 0.00000, 0.00000, 0.16742, -0.00000, -0.00000, 0.19706
3, -0.01250, -0.01250, 0.00000, -0.00835, -0.00835, 0.00000
4, 0.00000, 0.00000, -0.16742, -0.00000, -0.00000, -0.19706
5, -0.15115, -0.15115, 0.15115, -0.09815, -0.09815, 0.09815
6, 0.00000, 0.16742, 0.00000, -0.00000, 0.19706, -0.00000
7, -0.01250, 0.00000, -0.01250, -0.00835, 0.00000, -0.00835
8, 0.00000, -0.01480, -0.01480, -0.00000, -0.01824, -0.01824
9, -0.01507, 0.00000, 0.00000, -0.00000, 0.00000, 0.00000
10, 0.00000, -0.01480, 0.01480, -0.00000, -0.01824, 0.01824
11, -0.01250, 0.00000, 0.01250, -0.00835, 0.00000, 0.00835
12, 0.00000, -0.16742, 0.00000, -0.00000, -0.19706, -0.00000
13, -0.15115, 0.15115, -0.15115, -0.09815, 0.09815, -0.09815
14, 0.00000, 0.01480, -0.01480, -0.00000, 0.01824, -0.01824
15, -0.01250, 0.01250, 0.00000, -0.00835, 0.00835, 0.00000
16, 0.00000, 0.01480, 0.01480, -0.00000, 0.01824, 0.01824
17, -0.15115, 0.15115, 0.15115, -0.09815, 0.09815, 0.09815
18, 0.16742, 0.00000, 0.00000, 0.19706, -0.00000, -0.00000
19, 0.00000, -0.01250, -0.01250, 0.00000, -0.00835, -0.00835
20, -0.01480, 0.00000, -0.01480, -0.01824, -0.00000, -0.01824
21, 0.00000, -0.01507, 0.00000, 0.00000, -0.00000, 0.00000
22, -0.01480, 0.00000, 0.01480, -0.01824, -0.00000, 0.01824
23, 0.00000, -0.01250, 0.01250, 0.00000, -0.00835, 0.00835
24, -0.01480, -0.01480, 0.00000, -0.01824, -0.01824, -0.00000
25, 0.00000, 0.00000, -0.01507, 0.00000, 0.00000, -0.00000
26, 0.02364, 0.02364, 0.02364, -0.00686, -0.00686, -0.00686
27, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
28, 0.02364, 0.02364, -0.02364, -0.00686, -0.00686, 0.00686
29, 0.00000, 0.00000, 0.01507, 0.00000, 0.00000, 0.00000
30, -0.01480, 0.01480, 0.00000, -0.01824, 0.01824, -0.00000
31, 0.00000, 0.01250, -0.01250, 0.00000, 0.00835, -0.00835
32, 0.02364, -0.02364, 0.02364, -0.00686, 0.00686, -0.00686
33, 0.00000, 0.01507, 0.00000, 0.00000, 0.00000, 0.00000
34, 0.02364, -0.02364, -0.02364, -0.00686, 0.00686, 0.00686
35, 0.00000, 0.01250, 0.01250, 0.00000, 0.00835, 0.00835
36, -0.16742, 0.00000, 0.00000, -0.19706, -0.00000, -0.00000
37, 0.15115, -0.15115, -0.15115, 0.09815, -0.09815, -0.09815
38, 0.01480, 0.00000, -0.01480, 0.01824, -0.00000, -0.01824
39, 0.01250, -0.01250, 0.00000, 0.00835, -0.00835, 0.00000
40, 0.01480, 0.00000, 0.01480, 0.01824, -0.00000, 0.01824
41, 0.15115, -0.15115, 0.15115, 0.09815, -0.09815, 0.09815
42, 0.01480, -0.01480, 0.00000, 0.01824, -0.01824, -0.00000
43, 0.01250, 0.00000, -0.01250, 0.00835, 0.00000, -0.00835
44, -0.02364, 0.02364, 0.02364, 0.00686, -0.00686, -0.00686
45, 0.01507, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000
46, -0.02364, 0.02364, -0.02364, 0.00686, -0.00686, 0.00686
47, 0.01250, 0.00000, 0.01250, 0.00835, 0.00000, 0.00835
48, 0.01480, 0.01480, 0.00000, 0.01824, 0.01824, -0.00000
49, 0.15115, 0.15115, -0.15115, 0.09815, 0.09815, -0.09815
50, -0.02364, -0.02364, 0.02364, 0.00686, 0.00686, -0.00686
51, 0.01250, 0.01250, 0.00000, 0.00835, 0.00835, 0.00000
52, -0.02364, -0.02364, -0.02364, 0.00686, 0.00686, 0.00686
53, 0.15115, 0.15115, 0.15115, 0.09815, 0.09815, 0.09815