Coupled differential equation using ODE 45

My program is attached. I am getting dwp0dz and dwp2dz zero. Please help.

[James Tursa]

Please delete the image and instead post your code as text highlighed by the code button. We can't copy & run pictures.

[MAMTA SINGH]

I have removed image and inserted the code.

Answer 1

[Jan]

Yes, of course the values are 0.

The initial value of y(1) is 0. The derivative of y(1) is:

dwp0dz = ((s * alpha0 * y(1) * x0)/w);

This is 0, when y(1) is 0. So the value remains at 0.

The same happens for y(2):

dwp2dz = x1 *((alpha0 * y(1) * xr ) + (alpha0 * y(2) * x0) + (alpha2 * y(1) * x0))
%                       == 0                    == 0                   ==0

If all terms are 0, the sum is 0 also.

This means, that the calculations do exactly, what is expected. Why do you think that you need help?

A hint: c = 2.997 * 10^10 is a multiplication and an expensive power operation, while c = 2.997E10 is a cheap constant.