Matlab code to solve an integral

Hello! I need some help with Matlab.

I would like to solve and plot the following function:

I know everything except for W and x. I want to plot x vs. W, and the values of x should go from 0 to 1.

kw = 91371;

Ca0 = 1.81 * 10 ^ -4;

Fa0 = 0.1233;

K = 2;

Thank you in advance.

[Walter Roberson]

The value of the expression goes to infinity as x approaches 1, and the integral does as well. The integral is

(4116815215161981957*A)/1525375194746818396160000 - (2711395674095616*log(1 - A))/363677786528305625

Answer 1

[KALYAN ACHARJYA]

#Edited

syms x
kw = 91371;
CaO= 1.81 * 10 ^ -4;
FaO= 0.1233;
K = 2;
fun = (FaO*(1+(K*CaO)*(1-x)))/(kw*CaO*(1-x))
int(fun)

[madhan ravi]

+1 fplot()ing the result with the range 0-1 would give the desired plot. I suspect Fa0 is common for both numerator and denominator can’t verify as I don’t have PC to verify. The first paranthesis "(" should be after the first asterisk , I suppose.

[KALYAN ACHARJYA]

syms x
kw = 91371;
CaO= 1.81 * 10 ^ -4;
FaO= 0.1233;
K = 2;
fun = (FaO*(1+(K*CaO)*(1-x)))/(kw*CaO*(1-x))
in=int(fun)
fplot(in,[0 1])

I suspect Fa0 is common for both numerator and denominator, @Madhan I am not sure.

[Walter Roberson]

This is incorrect code. The 1+ K*Ca0 is not being multiplied by (1-x) in the original formula: it is 1 + (K*Ca0)*(1-x) .

[KALYAN ACHARJYA]

@Walter Yes, Edited Thanks sir