悬赏人民币，求教garch模型参数估计的问题

我想用一个类似pascal的程序语言，编写一个估计ar（1）——garch（1,1）的模型参数的程序，ar（1）没有截距项参数，最终目的就是我输入原始数据序列，输出各个参数的估计值，请问有没有可能，最终会有多少行代码？不用eviews直接去估计的原因是，这个程序语言可以实时更新原始数据，持续自动计算。

但是我看了好多计量经济学的书，一到garch模型的参数估计这块，都是一笔带过，说太复杂，要用数值方法，根本不仔细讲，我完全不知道具体估计这些参数的方法，所以现在想找一个完全了解用 marquardt算法（修正的BHHH，eviews的算法）进行garch参数估计的高手，通过qq指导我具体估计ar（1）——garch（1,1）的步骤，当然编程我自己来，你只要告诉我每一步都要干什么就行。

指导完成后，当然过程可能很麻烦，耗时也比较多，我完全没有上过数值分析的课程，基础只有计量经济学

酬谢500——1000元人民币（有商量），有在校同学愿意的话，请联系qq：285102759，酬谢决不食言

或者有高手用excel 用marquardt估计一个几十行原始数据的ar（1）-garch（1,1）的参数，当然迭代步骤有2、3次就行

随便找个平稳序列就行