How to plot this function?

Hello,

Could you please help me plot this function with Matlab ?

<http://img638.imageshack.us/img638/3205/fonction.png >

I'm a beginner. Your help will be valuable to me.

Many thanks

Answer 1

[Matt Fig]

If you want to do it manually, try this:

f = @(z) quadv(@(x) (1/(3*sqrt(2*pi)))*exp(-(x-30).^2/18),z,100);
Z = 0:.25:50;
G = zeros(1,length(Z));

for ii = 1:length(Z)
   G(ii) = f(Z(ii)); 
end

plot(Z,G)
xlabel('z')
ylabel('(1/\surd(2*pi))\inte^{-(x-30)^2/18}')

Answer 2

[Matt Tearle]

Take a look at erfc. A simple change of coordinates should get you there.

Then plotting can be done by

x = linspace(xmin,xmax,npoints);
y = constant*erfc(x);
plot(x,y)

[Michaël]

Thank you for your answer.

But can't I put my function directly (without going through a pre-determined function)?

Thanks

[Matt Tearle]

Why would you want to? This function is defined via an improper integral, which means you'll have to do an infinite-domain integral for each value of z. This will be slow and will introduce all sorts of numerical errors. On the other hand, the error function is a well-known function -- its implementation in MATLAB is stable and tested. Why reinvent the wheel?

But if you really want to, you can use something like quad to do the integral.

Answer 3

[Michaël]

Thank you for your answer.

But can't I put my function directly (without going through a pre-determined function)?

Thanks

Answer 4

[Michaël]

How can't we put : (-1+x.)^2 ?

Answer 5

[Michaël]

Thank you very much !!!