not able to solve the integration function

Dear Sir/Madam,

I am not able to solve the given integration function. The screenshot of given is attached and the code is given below

Thank you

for j=1:100
x(j)=linspace(0,0.025,101);
R(j)=linspace(0,0.4,101);
f=0.025;
syms x
b(j)= int(R(j)-sqrt(R(j).^2-x(j).^2),0, f);
c(j)=sqrt((1/f)*b(j));

end

[Kundan Prasad]

This type of error is seen

Unable to perform assignment because the indices on the left side are not compatible with the size of

the right side.

Error in diffra (line 46)

x(j)=linspace(0,0.025,101);

[Walter Roberson]

Solve for what value?

Is σ a constant?

Is the task to solve for f for each of a list of R values ?

Your vectors are length 101; why do you only loop to 100 ?

[Kundan Prasad]

σ is not constant and it varies with x.

The value of R can be taken as constant i.e. 0.4 and x should be varied from 0 to 0.025.

the upper limit of f=0.025 and lower limit is 0

I think for loop should be avoided.

[Walter Roberson]

What is being solved for?

Is the question how to calculate the right hand side given R, f, and symbolic x, outputting sigma ?

Normally when the phrase "integral equation" is used, it means that some variable inside an integral or some bound of the integral needs to be modified to find the correct value to match the other side of the equation. For example,

syms x b
solb = vpasolve(int(cos(x)^2 + tan(x/10), x, 0, b) == 1)

solb = 

1.966415013513613065429486177359

[Kundan Prasad]

yes sir

for given R=0.4, f=0 & 0.025, x=0:0.001:0.025,

what will be the value of σ for each x

[Kundan Prasad]

here f=0.025

Answer 1

[Jan]

This cannot work:

x(j) = linspace(0,0.025,101);

On the left is a scalar and on the right a vector with 101 elements.

Seeing only the invalid code does not clarify its intention. What is the purpose of the loop?

[Kundan Prasad]

what need to be change in code. Please clarify??

Thank you

[Jan]

I cannot fix the code such that it does, what you want, as long as you do not clarify what you want to achieve. I ask again: What is the purpose of the loop? I'd expect RT to be a constand and x to be a symbolic variable. Then there is no need for a loop.

Answer 2

[Tesfaye Girma]

indeces should have to match in both side and be compatable as well