clear all;
close all;
clc
sigma1=1;mu1=0;mu2=0;sigma2=1;p=0.5;
y=-4:0.1:4;
syms x
f=1/(2*pi*sigma1*sigma2*sqrt(1-p*p))*exp(-1/(2*(1-p*p))*(((x-mu1).^2)./(sigma1*sigma1)-2*p*((x-mu1)*(y-mu2))./(sigma1*sigma2)+((y-mu2).^2)./(sigma2*sigma2)));
a=int(f,-inf,3);
b=double(a);
plot(y,b),這樣能畫出個圖來.你求積分是對x求的,怎么畫圖里面是x呢?是不是搞錯了,我給你畫（y,a）就可以畫出來了