The Merton Jump-Diffusion Model
Merton (1976) [1] jump-diffusion model is an extension to the Geometric Brownian Motion model, with the underlying asset exhibit jumps in addition to have continuous diffusion paths. The asset price evolves as 
dS t =μS t dt + σS t  dW t  + (η−1)dq  
where μ is the drift rate, σ the volatility of S  and dW t   the increment of a Wiener process. The independent Poisson process dq  has a value of 1  with probability λdt  and 0 otherwise. 

The following is a ThetaML implementation of the Merton Jump-Diffusion Model. 
Merton jump diffusion Process simulation
Implementation as ThetaML
model JumpDiffusionProcess
%Jump diffusion model for stock price movement
%Example input parameters:
%S0 = 100, r = 0.03, sigma = 0.2, lambda = 1
......................

#define _CRT_SECURE_NO_WARNINGS

#include<stdio.h> 

void swap(char *a,char *b) 
{ 
char t; 
t = *a; 
*a = *b; 
*b = t; 
} 
void print(char *a,int length,int* t) 
......................

#include"stdio.h"
void main(){
int a,b,c,sum,average,product,smallest,largest;
printf("please input three numbers:\n");
scanf("%d%d%d",&a,&b,&c);//*输入三个整数*//

printf("Input three different integers:%3d%3d%3d",a,b,c);
sum=(a+b+c);//*求和*//

average=(a+b+c)/3;//*求平均值*//

product=a*b*c;//*求三个数的乘积*//
......................

public struct STONES
    {
        public int x;
        public int y;
    };
    public struct Lines
    {
        public STONES[] Point = new STONES[2];//错误
        public int[] items= new int[6];//错误
        public int value;
        public int color;
    };

#define _CRT_SECURE_NO_WARNINGS

#include <stdio.h>
#include <string.h>

void dec2yyy(int source,char* bin)
{
bin[8]='\0';
int j=7;
    while(source>0)
{
bin[j]=source %2+0x30;
......................

// t5.cpp : 定义控制台应用程序的入口点。
//
#define _CRT_SECURE_NO_WARNINGS


#include <stdio.h>
#include <string.h>

int main()
{
FILE* pfile;
int rows[20000][2];
......................

#===============================================================================
# author: yangjie
# mobilePhone:15850792203
# QQ:691633493
# E-mail:691633493@qq.com
# version:2.0
#===============================================================================
from maya.cmds import*

def YJRigZeroGrp(TargetObjList):
    if TargetObjList == 'select':
        TargetObjList = ls(sl=1)
......................

#===============================================================================
# ParentShape
# author: yangjie
# mobilePhone:15850792203
# QQ:691633493
# E-mail:691633493@qq.com
#version:2.0
#===============================================================================

from maya.cmds import*
def YJ_Parent_Shape(ShapeObj,TarObj,Sel,add):
    if Sel==1:
......................

#====================================================================================================
#Create_IKFK_RIG
#
#Designer:YangJie
#Email:691633493@qq.com
#Version:1.0
#============================================================================================
from maya.cmds import*
from YJ_FK_Rig import*

#Create window
if window('YJ_IKFK_UI',q=1,ex=1):
......................

#====================================================================================================
#Create_IKFK_RIG
#
#Designer:YangJie
#Email:691633493@qq.com
#Version:1.0
#============================================================================================
from maya.cmds import*
import maya.mel as mm
from YJ_Ctrl_Colour import*
from YJ_ParentShape import*
#YJ_CreatFkRig('joint1','joint3')
......................