calculate x value from normal distribution

I have data x , with nonzero mean. I can plot normal distribution with following

x=randi([1 20],1,60);
a=mean(x);
b=sqrt(var(x));
x= [-130:1:120];
norm = normpdf(x,a,b);
plot(x,norm)

to calculate x value for 97%, I have to see z score for 97% and then use

z=(x-mean)/standard deviation

for calculation of x value but how can I calculate x value directly by code.

Answer 1

[Torsten]

You mean

x = norminv([0.015 0.985],a,b)

?

Or simply

x = norminv(0.97,a,b)

?

Best wishes

Torsten.

[Jay Hanuman]

yes this is thank u.