贝贝 的个人资料凡尔赛玫瑰照片日志列表更多 ![]() | 帮助 |
|
2009/7/6 kaist博士申请结果--录取等待了2个月的结果。
很兴奋, 收到录取通知的那一天。 终于可以有理由逃离KIST了。
有了更好的学习环境。
刚刚给那边的教授写信, 希望能在他的实验室里继续haptics专业的学习。
希望KAIST的环境会比KIST要轻松。
再累都没有关系, 可是我真的很讨厌在巨大的压力下却不知道要做什么。
Give me a goal and a direction, I can get it.
2009/2/11 我的2008很久没有写日志了, 赶在2009新年到来前记录一下最近的生活. 首先最值得高兴的一件事情是我的第二个小侄子出生了, 呵呵. 知道消息以后乐了一个晚上. 我又做姑姑了, 真想回家过年, 等我明年回家的时候估计他都可以说话了吧. 2个周以前才开始工作, 过年只有规定的3天假期, 本来想多请几天家回家过年, 可是我老大没有同意, 所以我也没问姜博士免得挨批. 反正9月份刚刚回去过, 再回去的话被有可能被当成间谍啊. 昨天听办公室他们说这边的邮件和电话都是被NIS监听的, 不过估计是逗我呢. 还有一个原因就是签证的问题, 换签证需要2个周的时候, 过完年以后才能拿到.... 第一次在外面过年. 好在这边的外国朋友很多, 也有几个中国朋友留下来过年, 到时候可以一起看看春节晚会, 买点儿만두煮一下, ^_^, 也很不错啊 帖几张照片把2008做个总结,希望新的一年开开心心. 1.首先是宿舍生活,2008年住进了kut的国际学生宿舍并做英语助教. 在KUT国际宿舍住的这一年是自高中开始的宿舍生活中最充实最开心也最忙碌的一年. 在这里有来自菲律宾的5个交换学生, Ron, Mark, Sara, Rvon, Weana, 每个人都是个乐天派, 也把菲律宾人的热情带到了国际宿舍的生活中.他们能歌善舞总是把大家带的很活跃. 来自澳大利亚的管理员夫妇John和elizabeth. John总是开一些澳大利亚人特有的幽默, 不过我的英语水平还是不够高很多时候不知道他在说什么^_^. 再有就是来自俄罗斯的ildar, khalis, igor, 还有igor的女朋友, 因为igor继续读博才来韩国学习的Anastasia. 然后是越南的hai, hai的女朋友Antra, 和hai_bac. 其实kut的爱情故事还是不少的, 我就曾经被hai和女朋友感到过, 希望他们一直幸福下去. 来自埃及的Ghaly, KUT留学生当中唯一的一个父亲级的, 总是像大哥似的对其他人很照顾. 还有包括我在内的5个中国人.
国际宿舍staff2008年末合影
宿舍活动合影 2. 实验室生活 呆在韩国的这两年80%的时间是在实验室度过的, 几乎跟家一样. 没有什么事情也喜欢呆在实验室里, 上上网, 和hai或者其他的朋友聊天. 最值得感谢的一个人是??, 因为有她在实验室而很开心. 她给予我的不仅是帮助还有朋友的关心. 本来我比她早一个学期毕业, 可是她先找到了工作-在国家航天局, 很替她高兴,不过她要走的时候大家很伤心, 互相安慰说以后还是可以每天打电话的啊, 周末也可以见面. 可是工作以后才知道相互有多忙, 别说见面了打电话的时候都很少. 不过想起她来还是很温暖的... 两年来我最亲密的朋友. 要一直幸福哦....
07年夏天在济州岛开会时的合影
汉拿山顶峰合影. 我们花了将近5个小时才登上山顶,下山又用了3个多小时. 登上山顶真的很开心啊
08年春天实验室滑雪
滑雪后聚餐, 烤肉+烧酒
08年夏天登山 值得纪念的事情好多.... KUT实验室的生活已经告一段落, 不过在韩国的生活还会继续一段时间, 继续像以前快乐的生活... 3. 研究生毕业 08年夏天研究生毕业了. 毕业好难, 都不敢轻易开始博士生活了, 很担心毕不了业啊. 从山东大学一起来韩国的朋友一起毕业了. 半年前还凑在一起商量着怎么写论文, 以后要做什么, 现在都有了各自不同的生活了, 有的回国工作了, 有在韩国工作的, 还有去了日本继续读博的, 还有的留在了KUT继续博士生活. 希望大家一切顺利
电子系的同学们
活宝3人团合影 3. ICCAS2008论文发表 毕业以后继续在实验室留了半年, 做实验写论文. 只是半年的成果因为投稿失误全部丧失了. 幸运的是夏天写的一篇论文被ICCAS录取了, 于是有了我的第一次国际会议发表的机会. 发表之前的两个星期很紧张, 教授盯着检查PPT, 还要在实验内部seminar上做演练, 好几次都没有通过, 直到发表前一天还在seminar室里对着教授又练习了一遍, 不过发表的时候倒不紧张了, 呵呵. 主席还把我的名字误以为是男生.
感谢教授在学术上的指导 4. KIST 工作 面试后等了两个月终于在09年元旦以后开始了在KIST的工作. 一切也基本安定下来了, 剩下的就是努力工作啦. 在这边也很不错, 不过刚开始还不是很适应, 很想念在KUT有朋友在一起的日子. 这边比较好的一点是可以自己做饭. 最近发现一家有卖各种蔬菜的超市, 把自己喂的很饱.
姜博士请吃饭
KIST宿舍 这么写来08年幸福的事情好多. 希望09年一切顺利, 开开心心, 还有就是希望写09年总结的时候能出现男友的身影, 呵呵 2008/10/28 非应届生的应届生最近几天一直在网上搜索招聘信息, 在不同的招聘网站上填写简历, 徘徊于各个公司的招聘网站.
简历, 求职信, 简历.... 原来找工作也是门学问, 也是如此的耗时繁琐.
投了不少简历出去, 不过至今还没有什么消息, 在等着人事经理筛选完了的判决.
国内各大公司的校园招聘已接近尾声, 发现我这个悬在空中的"学生"处于非常尴尬的境地:既不是应届生也没有任何的工作经验.
在惊叹于自己08年夏天刚刚毕业却为何不能称为应届生的时候惊醒到应届生指的是那些09年毕业的还在校园里的弟弟妹妹们. 回想自己应该享受应届生待遇的去年,那时的自己还在为project, 每周一次的semianr而夜夜守在电脑前, 大脑里还没有工作的概念.其实到现在工作的概念还没有完全形成.
不过话说回来对在国外求学的我们来说只有毕业了才能开始找工作. 毕业前一连串的实验,项目,论文,答辩, 工作的事情完全被压在了脑后. 不是没想到去找工作, 而是很迷茫, 不知道怎么找工作, 到哪里找工作, 在哪里工作, 甚至思考到最后连自己想要什么都迷茫了.
尽管在韩国呆了两年也已毕业, 在这里生活不方便的地方也还不是一点儿两点儿. 就找工作这件事情来说, 甚至连招聘的消息都听不到. 也难怪, 处的地理位置太局限是一方面原因, 学校不为外国学生提供任何就业信息也是件让人匪夷所思的事情. 只管招收学生却连在韩国工作机会的信息也不提供. 只好自己动手丰衣足食了.
其实每年的9月底到11月初是找工作的黄金时期, 当然是指应届生啦. 各大公司,纷纷到各高校开宣讲会,吸引人才, 当然最后有多少人能够有幸进入这些名企工作就不得而知了. 对于在韩国的中国学生来说这段时间也是个机遇, 各大韩企也会在自己的招聘主页上贴出招聘启示, 期中有专门针对外国人的招聘. 最近我看到的有sk, 三星等.
只是注意到的比较晚, 很多公司的招聘启示已经关闭, sk的是到27号下午5点截至, 我看到的时候已经是27号晚上11点了. 提醒一下有意向进韩企工作的要提前注意啦. 也有一些公司还没有关闭, 三星也开放了一些针对外国人的职位到30号截至.
这里把知道的一些招聘的网址帖一下, 希望给有心人提供一些消息:
SK的招聘网站, 针对中国人的招聘27号已经结束, 不过可以继续关注其他岗位的招聘信息 http://www.dearsamsung.com.cn/front/web/index.aspx
我是应届生2009名企校园招聘 http://www.54yjs.cn/2009/
前程无忧校园招聘 http://xy.51job.com/
暂时提供的消息就这些, 其他的招聘网站如中华英才网, 应届生就业网也不妨常去转转.
好运! 2008/3/20 献血记下午实验室的朋友想去献血, 问我有没用时间一起去。
我对献血持保留态度, 至于在韩国献血更需要考虑一下了。
食堂前面有两辆献血车, 需要先登记一下, 在表格上填写资料。
由于表格是韩语的, 我们两个用英语一边讨论一边填, 写到一半的时候, 一个西装革履的中等身材的50岁左右的韩国人走过来问朋友你从哪个国家来的, 在韩国多长时间了, 朋友回答说是越南人来韩国快2年了, 我说是中国人。 然后他说越南和另外一个国家(没听清是哪儿)不能在韩国献血。 然后说中国人可以, 我直接对他说我没打算献血是陪朋友来的。
作为外国人感觉受到了歧视。
先不说有没有这种规定, 即使有的话还问来韩国多长时间干嘛。
去献血本来是个积极的事情, 结果产生了反感。
2008/3/18 战毒记为期一天半的战毒终于结束, 从昨天中午发现电脑里的杀毒软件开始瘫痪, 有预感中毒了。
接着在网上搜索相关的专杀工具, 凡是关键字中有杀毒软件名字浏览器自动关闭, 从别的电脑上下载的工具也不能运行, 打开任务管理器的话电脑运行速度奇慢并且有lsass.exe 和smss.exe两个以用户名存在的进程。
各项特征显示中了最新的磁碟机病毒了..... 电脑够赶潮流的, 啥厉害中啥......
回想一下, 好像是从网上下载软件的时候下载了setup.exe 而不是要下载的软件, 只要打开这个安装文件就中病毒了。
至于详细特征个大杀毒网站上都有, 提醒没中毒的同志们要小心, 如果发现setup.exe先检查一下大小, 若为100kB左右就要小心了。
隔壁实验室的朋友已经开始重装电脑。
经过各种尝试后发现360安全卫士的专杀工具能奏效,继而用其他杀毒软件进行扫尾工作。
杀病毒的感觉也不错哈....只是太浪费时间了......
虽然编写病毒也需要技术, 可那些人们也有点儿太缺了, 有这技术的话干点儿啥不好。 2008/3/14 终点后的起点第二季也结束了, 期待第三季继续演绎一段幸福的开始.
期待秀丽会有怎样的选择, 完全不同的路不管如何选择都是幸福的.
很累的时候想到一些事情总是幸福的吧, 令人质疑的是如果失去了永远的还会有这般的自信和微笑吗?
不管怎样选择自己内心想要的那种生活吧, 即便要放弃.........
十三姬能够抛却以前的记忆从新开始新的生活吗?
认真而努力的生活吧...
2008/3/13 Do not wait...Don’t wait for a smile to be nice...
不要等到了一个笑容才面露慈善 Don’t wait to be loved, to love. 不要等被爱了以后,才要去爱 Don’t wait to be lonely, to recognize the value of a friend. 不要等到寂寞了,才明白朋友的价值。 Don’t wait for the best job, to begin to work. 不要非要等到一份最好的工作,才要开始工作 Don’t wait to have a lot, to share a bit. 不要等拥有许多后,才开始分享一些 Don’t wait for the fall, to remember the advice. 不要等到失败后,才记得别人的忠告 Don’t wait for pain, to believe in prayer. 不要等到受伤了,才相信愿意祈祷 Don’t wait to have time, to be able to serve. 不一定要等到有时间,才能够去付出服务 Don’t wait for anybody else pain, to ask for apologies... 不要等别人受伤了,才来乞求原谅 ... neither separation to make it up. 不要等到分开了,才想到去挽回 Don’t wait... Because you don’t know how long it will take. 不要等待,因为,你不知道等待需要花费多少的时间 Remember: Friendship is like wine, it gets better as it grows older.
友谊像醇酒,越久越浓。 Love and TimeOnce upon a time, there was an island where all the feelings lived: Happiness, Sadness, Knowledge, and all of the others, including Love. One day it was announced to the feelings that the island would sink, so all constructed boats and left. Except for Love.
Love was the only one who stayed. Love wanted to hold out until the last possible moment. When the island had almost sunk, Love decided to ask for help. Richness was passing by Love in a grand boat. Love said, "Richness, can you take me with you?" Richness answered, "No, I can't. There is a lot of gold and silver in my boat. There is no place here for you." Love decided to ask Vanity who was also passing by in a beautiful vessel. "Vanity, please help me!" "I can't help you, Love. You are all wet and might damage my boat," Vanity answered. Sadness was close by so Love asked, "Sadness, let me go with you." "Oh . . . Love, I am so sad that I need to be by myself!" Happiness passed by Love, too, but she was so happy that she did not even hear when Love called her. Suddenly, there was a voice, "Come, Love, I will take you." It was an elder. So blessed and overjoyed, Love even forgot to ask the elder where they were going. When they arrived at dry land, the elder went her own way. Realizing how much was owed the elder, Love asked Knowledge, another elder, "Who Helped me?" "It was Time," Knowledge answered. "Time?" asked Love. "But why did Time help me?" Knowledge smiled with deep wisdom and answered, "Because only Time is capable of understanding how valuable Love is." 只有时间了解爱的价值 2008/3/12 春天来了, 动起来不知不觉2008年的春天已经来临了。
提到2008总不禁与奥运联系起来。
一年之际在于春, 该活动一下筋骨了, 有时间的话想多在健身房待会儿,当然要是能在春暖花开的日子和朋友在林荫路上慢跑的话最好不过了。又一次想起了玛珈山后临海的小路, 幽幽的好像永远都见不到阳光却总是清凉的即便是在炎炎的夏日。
目前为止还没有探索到适合长期室外跑步的地点, 硬是要挑选一个的话校外的村庄倒是可以去转转, 只是院子里的狗叫声比较破坏心情。 继续发掘中。。。
帖几个跑步的小贴士:
Step 1:
Jog before you run. Every running program, no matter what level, includes some jogging. Step 2:
Run at a conversational pace. While jogging, you should be able to talk comfortably without being winded. Step 3:
Mix running with walking, if necessary. Increase the amount of running and decrease the walking as you progress. Step 4:
Be patient. Understand that although the heart and lungs respond quickly to exercise, muscles and joints take much longer. Step 5:
Build up your running to at least 20 to 30 minutes three times a week, done at a moderate level of intensity. Studies have shown that this is a sufficient amount of exercise for basic cardiovascular fitness. Running more than this amount is for reasons beyond basic fitness. Step 6:
Increase your running by no more than 10 percent a week. This can be measured either in time or distance. Step 7: Use the first month to learn about yourself. Pay close attention to your body; learn to read its signals of fatigue and stress, and find out when you can push beyond them. 2008/1/25 暂别了kut, 2007年韩国的生活明天就要离开韩国, 开始为期三周的假期。
从拿到机票开始无数次盼着回家, 盼着跟父母朋友们见面, 盼着在家里暖暖的睡大觉, 当然还有想到就口水直流的饭菜
随着离开的日子一天天的接近, 竟然很舍不得, 舍不得在这儿的生活, 还有这里熟悉的一切。
2007年的生活很多彩, 感谢在韩国的朋友们, 感谢你们的关心帮助,留下了如此多的幸福的记忆。
不知道yeasoul能不能看到这篇日志, 当然看到的话也不知道我写的什么吧, 虽然我们语言不通很多话不知道怎么交流, 感谢你的理解和帮助, 想起你的话很温暖噢
当然还有严厉温和的教授, 希望您老人家保佑我明年按时毕业。
当然还有另外的特别要感谢的人和事,还有听我唠叨的朋友们, 辛苦你们的耳朵啦, 呵呵。
对大家说一声新年快乐, 来年见。
2008/1/2 新年快乐Time flys away.
日历转入2008年,赶在零点前许下了08年的愿望。
第一个愿望嘛, 秘密
第二个愿望呢, 还是秘密
第三个愿望。。。
愿望说出来就不灵了, 还是保密吧。
不管怎样, 新的一年里希望收获更多, 会一如既往的努力生活。 2007/11/13 100个名著浓缩的一句话1.神要是公然去跟人作对,那是任何人都难以对付的。 (《荷马史诗》) 2007/10/14 (转)在美国的困惑:回国还是不回?在美国呆的时间长了,你会发现,几乎你所遇见的每一个中国人,都是一个祥林嫂。 他们喋喋不休地反反复复地披星戴月地不断追问你追问自己:以后想不想回国?以后想不想回国?以后想不想回国?…… 2007/10/8 生命中的11种感动It hurts to love someone and not be loved in return. But what is more painful is to love someone and never find the courage to let that person know how you feel. 只有付出的爱是痛苦的,但比这更痛苦是爱一个人却没有勇气让那人知道你的感情。 A sad thing in is when you meet someone who means a lot to you,only to find out in the end that it was never meant to be and you just have to let go. 生命中令人悲伤的一件事是你遇到了一个对你来说很重要的人,但却最终发现你们有缘无份,因此你不得不放手。 The best kind of friend is the kind you can sit on a porch swing with,never say a word,and then walk away feeling like it was the best conversation you've ever had. 最好的朋友就是那种能和你促膝而坐,彼此不说只字片语,分别时却感到这是你有过的最好的一次交流! It's true that we don't know what we've got until we lose it, but it's also true that we don't know what we've been missing until it arrives. 的确只有当我们失去时才知道曾拥有的是什么,同样,只有当我们拥有了才知道曾经失去了什么。 It takes only a minute to get a crush on someone,an hour to like someone,and a day to love someone- but it takes a time to forget someone. 迷上某人只需一分钟,喜欢上某人需要一小时,爱上某人则要一天,然而,忘记某人却是一辈子的事情。 Don't go for looks;they can deceive . Don't go for ;even that fades away. Go for someone who makes you smile because it takes only a smile to make a dark day seem bright. 别倾心于容貌,因为它具有欺骗性,也别倾心于财富,它也会消散,倾心于那个能带给你笑容的人吧,因为一个笑容能使漫漫长夜如白昼般明亮。 Dream what you want to dream;go where you want to go;be what you want to be,because you have only one and one chance to do all the things you want to do. 做你想做的梦吧,去你想去的地方吧,成为你想成为的人吧,因为你只有一次生命,一个机会去做所有那些你想做的事。 Always put yourself in the other's shoes. If you feel that it hurts you,it probably hurts the person too. 要设身处地的为别人着想, 如果一双鞋你穿着夹脚, 别人的感觉可能也一样。 A careless word may kindle strife;a cruel word may wreck a life;a timely word may level ;a loving word may heal and bless. 无心快语可能引发争执,无情之词可能折损生命,适时温语可能消弭压力,而关爱之声可能治愈心灵。 The happiest of people don't necessarily have the best of everything they just make the most of everything that comes along their way. 幸福之人并非拥有一切,只是尽力享受生活的赐予。 Love begins with a smile,grows with a kiss,ends with a tear. When you were born,you were crying and everyone around you was smiling. Live your so that when you die,you're the one smiling and everyone around you is crying. 爱情以笑开始,以吻转浓,以泪结束。当你哭着降临人世时,身边的每个人都在为此欢笑,好好生活吧,这样你就能含笑离开人世,而身边的每个人都在为此哭泣。 2007/9/5 济州岛旅行归来 因为8月份有机器人会议在济州岛召开, 在结束暑假的学习和seminar以后我们实验室一行5人去了济州岛--开会加旅行.
星期六早晨7点钟从学校出发, 先开车到了韩国南端的一个港口완도, 然后坐船3个小时到达济州岛.中午在완도吃饭的时候, 特意找了一家中华料理店, 想吃中国菜来着, 可是店里一个汉字也没有, 更别说中国菜了, 最后只吃了一碗豆汤面.
到达的当天晚上去了翰林公园附近的海水浴场. 这个海水浴场虽然不大但是很漂亮, 离海岸不远的地方是个小岛, 游泳过去也不过半个小时的时间. 济州岛的海很清感觉一眼就能看穿海底似的, 天空湛蓝, 石头是黑色的. 在济州岛上随处可见石头制成的돌하르방, 非常可爱的济州守护神.
后来去了好多旅游景点, 电视剧大长今的外景拍摄地, 还有海水浴场, 牛岛, 中文观光园区, 龙头岩等等, 还去爬了汉拿山, 据说是韩国最高的山, 有2000多米, 上山下山整整用了8个小时, 下来以后感觉腿都不是自己的了,没有知觉的在走路, 反正把我给累坏了. 虽然累, 但是到达山顶的心情也不是一般的好.在山外面什么也看不到只是满山郁郁匆匆的树木, 到了1500米左右的时候视野豁然开阔, 陡峭的山崖, 脚下飘着的白云, 远处的城市和大海, 遍布山野的树木和花草, 还有潺潺的溪水. 下次有机会真想再去看一看. 更有意思的是我们登山的时候竟然同行的是韩国的一队篮球运动员, 可是没有机会跟他们说上话, 唉, 这么好的机会给错过了,呵呵. 2007/8/19 今天是七夕情人节今天是中国传统的情人节-- 七夕, 传说中牛郎和织女相会的日子.
祝愿天下有情人终成眷属.
贴一首关于爱情的歌--Nine million bicycles,
Nine Million Bicycles There are nine million bicycles in Beijing 北京城里有900万辆自行车 That's a fact 这是个事实 It's a thing we can't deny 是我们无法否认的事 Like the fact that I will love you till I die 就像我会爱你直到死去 We are twelve billion light years from the edge 我们离时间尽头有120亿光年 That's a guess 那只是个猜测 No-one can ever say it's true 没人能说那是真的 But I know that I will always be with you 但我知道我会一直跟你在一起 I'm warmed by the fire of your love everyday 你的爱天天都温暖着我 So don't call me a liar 别说我是说谎者 Just believe everything that I say 要相信我所说的一切 There are six billion people in the world 在这世上有60亿人 More or less 差不多 And it makes me feel quite small 这让我感觉挺渺小 But you're the one I love the most of all 但你是我最爱的人 We're high on the wire 我们站得很高 With the world in our sight 眼里能看见全世界 And I'll never tire 而且我从不厌倦 Of the love that you give me every night 你所给我的爱 There are nine million bicycles in Beijing 北京城里有900万辆自行车 That's a fact 这是个事实 It's a thing we can't deny 这事我们无法否认 Like the fact that I will love you till I die 就像我会爱你至死 And there are nine million bicycles in Beijing北京城里有900万辆自行车 And you know that I will love you till I die 而你也知道我会爱你直到死去 2007/8/4 Runge-Kutta algorithm (龙格-库塔算法)上一次seminar的任务是用C++实现动态系统的仿真而不是用simulink.
一直都用simulink做仿真, 如果知道系统的动力方程的话只需要拖拽就能建立系统模型. 如果用C++仿真的话, 先不说我那很菜的编程了, 单单是求解二阶微分方程的算法就折腾了一晚上. 虽然看了教授编的程序例子, 还是不能理解, 后来听实验室的俄罗斯同学说教授用的是四阶龙格-库塔算法求解微分方程. 并且对我连这个名字都不知道感到很奇怪, 他说这是个常用的算法, 动态系统仿真经常使用的. 他说他们大学一年级的时候上物理课就用这种算法编程仿真系统的, 精确度很高.
听完我出了一身冷汗. 大学一年级..............., 我们的大学一年级做什么了, 好象山东大学威海分校电子系大学一年级的学生还在学习大学语文, 美术欣赏, 音乐欣赏, 政治思想教育 之类的"文化" 课程吧, 至于专业课, 大概就学了个数学和物理, 还有电路基础. 至于编程是什么, 一年级的学生好多连电脑操作都不熟悉, 更不用说龙格-库塔算法了. 汗啊..................当然这相当程度上也与国情, 学校情况有关, 但是毕竟落后很多啊至少在教育上.
发这些感慨也没有, 我把C++编的算法贴在这里, 希望有人看到的话用到的话做个参考. 至于原理嘛, 我也没有彻底明白, 学习当中
/* Runge Kutta for a set of first order differential equations */
#include <stdio.h>
#include <math.h> #define N 2 /* number of first order equations */
#define dist 0.1 /* stepsize in t*/ #define MAX 50.0 /* max for t */ FILE *output; /* internal filename */ main()
{ double t, y[N]; int j; void runge4(double x, double y[], double step); /* Runge-Kutta function */ double f(double x, double y[], int i); /* function for derivatives */
output=fopen("osc.dat", "w"); /* external filename */
y[0]=1.0; /* initial position */
y[1]=0.0; /* initial velocity */ fprintf(output, "0\t%f\n", y[0]); for (j=1; j*dist<=MAX ;j++) /* time loop */ { t=j*dist; runge4(t, y, dist); fprintf(output, "%f\t%f\n", t, y[0]);
} fclose(output);
} void runge4(double x, double y[], double step)
{ double h=step/2.0, /* the midpoint */ t1[N], t2[N], t3[N], /* temporary storage arrays */ k1[N], k2[N], k3[N],k4[N]; /* for Runge-Kutta */ int i; for (i=0;i<N;i++) t1[i]=y[i]+0.5*(k1[i]=step*f(x, y, i)); for (i=0;i<N;i++) t2[i]=y[i]+0.5*(k2[i]=step*f(x+h, t1, i)); for (i=0;i<N;i++) t3[i]=y[i]+ (k3[i]=step*f(x+h, t2, i)); for (i=0;i<N;i++) k4[i]= step*f(x+step, t3, i); for (i=0;i<N;i++) y[i]+=(k1[i]+2*k2[i]+2*k3[i]+k4[i])/6.0;
} double f(double x, double y[], int i)
{ if (i==0) return(y[1]); /* derivative of first equation */ if (i==1) return(-0.2*y[1]-y[0]); /* derivative of second equation */ } |
|
|