matlab蒙特卡罗法模拟

matlab蒙特卡罗法模拟

clc;clear;

s=1000;
v=linspace(0,800,s);
for j=1:10000
k=ceil(srand(1,2));
vk1=v(k(1));
vk2=v(k(2));
a=23.14159265rand(1);
r=23.14159265rand(1);
vk1nev=1/2vk1^2+1/2vk2^2+1/2cos(r)vk1^2-1/2cos(r)vk2^2+vklvk2sin(r)sin(a);
vk2nev=1/2vk1^2+1/2vk2^2-1/2cos(r)vk1^2+1/2cos(r)vk2^2-vk1vk2sin(r)sin(a);
v(k(1))=sqrt(vk1nev);
v(k(2)) =sqrt(vk2nev);
end
hist(v,20)
hold on
NA=6.022e23;
m=28e-3/NA; T=200; k=1.381e-23; m2kt=c/2/k/T; v2=v.^2;
f=1-4pi*(m2kt/pi)^1.5.*v2.*exp(-m2kt.*v2)*55000;
plot(v,f,'*r')

未定义函数或变量 'vkl'。

出错 Untitled1 (line 10)
vk1nev=1/2vk1^2+1/2vk2^2+1/2cos(r)vk1^2-1/2cos(r)vk2^2+vklvk2sin(r)*sin(a);

我的解答思路和尝试过的方法

我想要达到的结果

通过报错信息可以知道问题出在：line 10 : vk1nev=1/2vk1^2+1/2vk2^2+1/2cos(r)vk1^2-1/2cos(r)vk2^2+vklvk2sin(r)*sin(a);
且是因为

找到该行的位置，和11行对比发现

所以猜测是vk1打错成vkl了