y=[0.9,0.95,0.975,0.99] % alphas
tvar=[icdf(t10,y);icdf(t50,y);icdf(t100,y)] % VaR of t distribution
format long
x=transpose(1.6:0.0004:4.9) % icdf function doesn’t apply to mix. normal distribution 
nvar10=[x,cdf(gm10,x)] % VaR manually searched   
x=transpose(1.3:0.0004:4.4) 
nvar50=[x,cdf(gm50,x)] 
x=transpose(1.4:0.0004:4.3) 
nvar100=[x,cdf(gm100,x)]
