含变量的二维正态分布的概率密度函数

clc
close all
D=5;%
deltaD = 1/100D+0.1；
delta_theta = 3pi/180;
sigma_mz = (deltaD^2+(D^2*delta_theta^2))^(0.5);
theta = 0:pi/180:pi/2;
sigma_x = ((sin(theta).*deltaD).^2+(D.*cos(theta).*delta_theta).^2).^(0.5);
sigma_y = ((cos(theta).deltaD).^2+(-D.sin(theta).delta_theta).^2).^(0.5);
x = Dsin(theta);
y = Dcos(theta);
[X,Y]=meshgrid(x,y);
G = 1/(2pi.*sigma_x.*sigma_y).exp(-1/2(X.^2/sigma_x.^2+Y.^2/sigma_y.^2));
mesh(x,y,G)
title('二维正态分布的概率密度函数')
xlabel('X')
ylabel('Y')
不能运行出图，求助各位大神

你把公式打一下，方便画图，哪些是变量