HYDROGEN ENERGY

 // FOR HYDROGEN ENERGY .

Input –

function y1=f(x)

y1=((2/((a0).^(3/2)))*(exp(-x/a0)));

endfunction

function y2=g(x)

y2=(1/((2*sqrt(2))*((a0).^(3/2)))*(exp(-x/(2*a0)))*(2-(x/a0)))

endfunction

function y3=h(x)

y3=(2/((81*sqrt(3))*((a0.^(3/2)))))*(27-

18*(x/a0)+2*((x^2)/(a0.^2)))*(exp(-x/(3*a0)))

endfunction

a0=5

x=0:25*a0;

subplot(2,2,1)

plot(x,f)

xlabel("X-AXIS")

ylabel("Y-AXIS")

subplot(2,2,2)

plot(x,g)

xlabel("X-AXIS")

ylabel("Y-AXIS")

subplot(2,2,3)

plot(x,h)

xlabel("X-AXIS")

ylabel("Y-AXIS")