Integration and Differentiation in c#

Question

1.00/5 (1 vote)

See more:

Hi ,
I am solving Manning Equation in C# ,but facing problem in implementing Differentiation to the below

f(y)=[ (nq/(c2*s0^(1/2)))3/5)*(c2/c1)]5/8 to f'(y).

i am doing this in console application.can anyone help me with this please.

EDIT: AL
Added code from comments:

sir sorry ,i am not troubling i think ,

this the whole code of my project
and this is the class file
Equation1.cs

C#

public double Find_C1(int Z1,int Z2) 
{
  return (Z1 + Z2) / 2;
}
public double Find_C2(int Z1,int Z2)
{
  return c1=Math.Pow(((1 + Math.Pow(Z1, 2))), 0.5) + Math.Pow(((1 + Math.Pow(Z2, 2))), 0.5);

}
public double Aofy(int y,int b,double c1) 
{

  return a=((b + c1) * y) * y;
}
public double PofY(int y, int b, double c2) 
{
  return p = b +(c2 * y);
}
public double RofY(int y, int b, int c2) 
{
  return a / p;
}
public double find_V(double s0, double n, double R) 
{
  double v= (Math.Pow(s0,0.5) / n) * Math.Pow(R,0.66);
  return v;
}
public double find_q(double a,double v) 
{
  return a * v;
}
public double find_F(double n, double q, double c2, double s0, double c1) 
{
  double f;

  return f= Math.Pow((Math.Pow(((n * q) / (c2 * (Math.Pow(s0, 0.5)))), 0.6)) * (c2 / c1), 0.625);

}

and in main i am calling it as

C#

class Program
{
  static void Main(string[] args)
  { 
    Equation1 obj = new Equation1();
    Console.WriteLine(string.Format("Enter S0"));
    s0 = double.Parse(Console.ReadLine());
    Console.WriteLine(string.Format("Enter Z2"));
    Z2 = int.Parse(Console.ReadLine());
    Console.WriteLine(string.Format("Enter Z1"));
    Z1 = int.Parse(Console.ReadLine());
    Console.WriteLine(string.Format("Enter Y"));
    y = int.Parse(Console.ReadLine());
    Console.WriteLine(string.Format("Enter b"));
    b = int.Parse(Console.ReadLine());


    Console.WriteLine(string.Format("C1:{0}",obj.Find_C1(Z1, Z2)));
    Console.WriteLine(string.Format("C2:{0}", obj.Find_C2(Z1,Z2)));
   Console.WriteLine(string.Format("A(y):{0}",(b+(obj.Find_C1(Z1, Z2)*y))*y));
   Console.WriteLine(string.Format("P(y):{0}",(b+ (obj.Find_C2(Z1, Z2)*y))));
   Console.WriteLine(string.Format("R(y):{0}", ((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))));
   Console.WriteLine(string.Format("V:{0}", obj.find_V(s0,0.1,Math.Pow((((b+(obj.Find_C1(Z1, Z2)*y))*y)/(b+ (obj.Find_C2(Z1, Z2)*y))),2/3))));
   Console.WriteLine(string.Format("Q:{0}",obj.find_q(obj.Aofy(y, b, obj.Find_C1(Z1, Z2)), obj.find_V(s0, 0.1, Math.Pow((((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))), 2 / 3)))));


//F is Wrong
   Console.WriteLine(string.Format("F:{0}", obj.find_F(0.1, obj.find_q(obj.Aofy(y, b, obj.Find_C1(Z1, Z2)), obj.find_V(s0, 0.1, Math.Pow((((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))), 2 / 3))), obj.Find_C2(Z1, Z2), s0, obj.Find_C1(Z1, Z2))));
   Console.ReadLine();

}

public static double s0 { get; set; }

public static int Z2 { get; set; }

public static int Z1 { get; set; }

public static int y { get; set; }

public static int b { get; set; }
}

Posted 27-May-15 3:01am

RAJU-11052090

Updated 27-May-15 4:37am

Andy Lanng

v2

Add a Solution

Comments

Andy Lanng 27-May-15 9:11am

What have you tried? How are you parsing this?

RAJU-11052090 27-May-15 9:19am

sir there are several related expression like

to find Q=A*V
V=((S0)1/2/n)*R2/3
A=(b+c1y)y

C2=(1+(Z1)^2)^1/2 + (1+(Z2)^2)^1/2
C1=(Z1+Z2)/2

i am calculating all these separately and directly substituting the values

and S0,n,b,Z1,Z2 all are constants.

Andy Lanng 27-May-15 9:35am

You have to reply to my comment or I don't get notified that you have replied.

That's great, but show us some code. Click the 'Improve Question' button and add some code to the question ^_^

RAJU-11052090 27-May-15 9:51am

ok sir i'll

RAJU-11052090 27-May-15 9:51am

to find A

public double Aofy(int y,int b,double c1)
{

return a=((b + c1) * y) * y;
}

RAJU-11052090 27-May-15 9:52am

to find V;
public double find_V(double s0, double n, double R)
{
double v= (Math.Pow(s0,0.5) / n) * Math.Pow(R,0.66);
return v;
}

RAJU-11052090 27-May-15 9:54am

to find F(y) :

public double find_F(double n, double q, double c2, double s0, double c1)
{
double f;

return f= Math.Pow((Math.Pow(((n * q) / (c2 * (Math.Pow(s0, 0.5)))), 0.6)) * (c2 / c1), 0.625);

}

CPallini 27-May-15 9:35am

There is no y in your function.

RAJU-11052090 27-May-15 9:51am

i have to find out 'y'

in question i have this line
" Assume y0=2 and do iteration until the tolerance limit is 10^-4

Andy Lanng 27-May-15 10:01am

but f(y) needs to solve for y.
It's saying "function given y"
you define the function as "(nq/(c2*s0^(1/2)))3/5)*(c2/c1)]5/8"
you even start with y0=2
So y1 = f(y0).
so where does y0 fit in the equation?

RAJU-11052090 27-May-15 10:23am

sir i have separately added those sub equations

-->To find A

public double Aofy(int y,int b,double c1)
{

return a=((b + c1) * y) * y;
}

-->To find V;
public double find_V(double s0, double n, double R)
{
double v= (Math.Pow(s0,0.5) / n) * Math.Pow(R,0.66);
return v;
}

which is required to solve the main function.
and Y value will be assumed as 2 for first iteration and then it will run up to tolerance reach 10^ -4

RAJU-11052090 27-May-15 10:25am

Here we are finding A and V for the sake of Q.

Andy Lanng 27-May-15 10:33am

ah. ok - i see now ^_^ - thanks for the explanation.

So you need something like:
'{
''double nextY = 0.0
''double thisY = 2.0;
''double tolerance = thisY-nextY;
''while(tolerance > 0.0001)
''{
'''nextY=FofY(thisY);
'''tolerance = thisY-nextY;
'''if(tolerance <0.0)
'''{
''''tolerance = -tolerance
''''thisY = nextY
'''}
''}
'}
}

RAJU-11052090 27-May-15 10:48am

sir but can please explain this ,
i am not getting that thisy applying once again to f(y) without increment or decrement . please correct me if am wrong.

Andy Lanng 27-May-15 10:50am

oh, yup. You are correct. thisY = nextY should be in the loop. I have updated the comment. It's not too important. I fixed the bug in my solution.

RAJU-11052090 27-May-15 10:53am

i dint receive that updated comment sir :(

Andy Lanng 27-May-15 11:02am

never mind. I have updated the solution.

What values are you using for your tests?
If you tell me S0, Z2, Z1, Y and b then I'll see if I get the correct answer

RAJU-11052090 27-May-15 11:20am

S0=0.00001;
Z1,Z2=4;
Y from 2
And n =0.1

Andy Lanng 27-May-15 11:23am

'n'? where did that 'n' come from? there's a 'b'

RAJU-11052090 27-May-15 11:48am

Sir final f(y) we have to put n value f(y)=(nq/(c2*s0^(1/2)))3/5)*(c2/c1)]5/8

Andy Lanng 27-May-15 11:52am

YOUR code asks for b:
Console.WriteLine(string.Format("Enter b"));
b = int.Parse(Console.ReadLine());

YOUR code does not ask for n
what is the value for b?

RAJU-11052090 27-May-15 12:00pm

I am directly entering it sir

Andy Lanng 27-May-15 12:08pm

Ok i see n. n is there.

WHATS B!!!!!

(btw: I am not at all angry. I'm just teasing. But I still need to know what the value of the interger b is.)

when you run the app, you are asked to enter 5 variable values:
S0 is a double. I enter 0.00001. Then Z2 and Z1 are integers. I enter 4 for both. Y is also an integer. Y is 2. Then it asks for b. b is an integer.
"Enter b" it says. What should I type?
:P

RAJU-11052090 27-May-15 10:33am

sir sorry ,i am not troubling i think ,

this the whole code of my project
and this is the class file
Equation1.cs
------------------------------------
public double Find_C1(int Z1,int Z2)
{
return (Z1 + Z2) / 2;
}
public double Find_C2(int Z1,int Z2)
{
return c1=Math.Pow(((1 + Math.Pow(Z1, 2))), 0.5) + Math.Pow(((1 + Math.Pow(Z2, 2))), 0.5);

}
public double Aofy(int y,int b,double c1)
{

return a=((b + c1) * y) * y;
}
public double PofY(int y, int b, double c2)
{
return p = b +(c2 * y);
}
public double RofY(int y, int b, int c2)
{
return a / p;
}
public double find_V(double s0, double n, double R)
{
double v= (Math.Pow(s0,0.5) / n) * Math.Pow(R,0.66);
return v;
}
public double find_q(double a,double v)
{
return a * v;
}
public double find_F(double n, double q, double c2, double s0, double c1)
{
double f;

return f= Math.Pow((Math.Pow(((n * q) / (c2 * (Math.Pow(s0, 0.5)))), 0.6)) * (c2 / c1), 0.625);

}

-------------------------------------------------------
and in main i am calling it as

class Program
{
static void Main(string[] args)
{ Equation1 obj = new Equation1();
Console.WriteLine(string.Format("Enter S0"));
s0 = double.Parse(Console.ReadLine());
Console.WriteLine(string.Format("Enter Z2"));
Z2 = int.Parse(Console.ReadLine());
Console.WriteLine(string.Format("Enter Z1"));
Z1 = int.Parse(Console.ReadLine());
Console.WriteLine(string.Format("Enter Y"));
y = int.Parse(Console.ReadLine());
Console.WriteLine(string.Format("Enter b"));
b = int.Parse(Console.ReadLine());

Console.WriteLine(string.Format("C1:{0}",obj.Find_C1(Z1, Z2)));
Console.WriteLine(string.Format("C2:{0}", obj.Find_C2(Z1,Z2)));
Console.WriteLine(string.Format("A(y):{0}",(b+(obj.Find_C1(Z1, Z2)*y))*y));
Console.WriteLine(string.Format("P(y):{0}",(b+ (obj.Find_C2(Z1, Z2)*y))));
Console.WriteLine(string.Format("R(y):{0}", ((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))));
Console.WriteLine(string.Format("V:{0}", obj.find_V(s0,0.1,Math.Pow((((b+(obj.Find_C1(Z1, Z2)*y))*y)/(b+ (obj.Find_C2(Z1, Z2)*y))),2/3))));
Console.WriteLine(string.Format("Q:{0}",obj.find_q(obj.Aofy(y, b, obj.Find_C1(Z1, Z2)), obj.find_V(s0, 0.1, Math.Pow((((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))), 2 / 3)))));

//F is Wrong
Console.WriteLine(string.Format("F:{0}", obj.find_F(0.1, obj.find_q(obj.Aofy(y, b, obj.Find_C1(Z1, Z2)), obj.find_V(s0, 0.1, Math.Pow((((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))), 2 / 3))), obj.Find_C2(Z1, Z2), s0, obj.Find_C1(Z1, Z2))));
Console.ReadLine();

}

public static double s0 { get; set; }

public static int Z2 { get; set; }

public static int Z1 { get; set; }

public static int y { get; set; }

public static int b { get; set; }
}

Andy Lanng 27-May-15 10:34am

this should be added to your question. I'll do it

RAJU-11052090 27-May-15 10:35am

oh sorry sir

RAJU-11052090 28-May-15 1:07am

Thank you Very Much Sir. :)

RAJU-11052090 28-May-15 1:08am

Thank you Very Much.

Sergey Alexandrovich Kryukov 27-May-15 10:29am

Numerical (giving you particular figures for each separate data set) or analytic (Computer Algebra System, giving you the expression on output)?
—SA

RAJU-11052090 28-May-15 1:07am

Thank you Very Much

Sergey Alexandrovich Kryukov 28-May-15 1:09am

Hm... for what?
—SA

1 solution

Add a Solution

Add your solution here

Treat my content as plain text, not as HTML

Preview 0

…

Existing Members

Sign in to your account

...or Join us

Download, Vote, Comment, Publish.

Your Email
Password
Forgot your password?

Your Email
This email is in use. Do you need your password?
Optional Password

I have read and agree to the Terms of Service and Privacy Policy
Please subscribe me to the CodeProject newsletters

When answering a question please:

Read the question carefully.
Understand that English isn't everyone's first language so be lenient of bad spelling and grammar.
If a question is poorly phrased then either ask for clarification, ignore it, or edit the question and fix the problem. Insults are not welcome.
Don't tell someone to read the manual. Chances are they have and don't get it. Provide an answer or move on to the next question.

Let's work to help developers, not make them feel stupid.

This content, along with any associated source code and files, is licensed under The Code Project Open License (CPOL)

Andy Lanng · Accepted Answer · 2015-05-27T04:48:00

I got it. Thanks for posting the code. I know the error:

ints are problematic. you would think that (1/2) would equal 0.5? well c# says it's 1.

operations on ints returns an int.

C#

//Every time z1+z2 is odd, you will lose precision.
public double Find_C1(int Z1, int Z2)
{
   return (Z1 + Z2) / 2;
}


//Try this:
public double Find_C1(int Z1, int Z2)
{
   return ((double)(Z1 + Z2)) / 2.0;
}

//also in main:
...
Console.WriteLine(string.Format("V:{0}", obj.find_V(s0, 0.1, Math.Pow((((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))), 2 / 3))));
//2 / 3 = 1, not 0.66.
// try 2.0 / 3.0
...

Console.WriteLine(string.Format("V:{0}", obj.find_V(s0, 0.1, Math.Pow((((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))), 2.0 / 3.0))));
//the 2/3 appears a few time so make sure you catch them all

hope that helps ^_^

also - this is the iterative process to find yn to a tolerance of 0.0001:

C#

double nextY = 0.0;
double thisY = y;
double tolerance = thisY - nextY;
while (tolerance > 0.0001)
{
    nextY = obj.find_F(0.1, obj.find_q(obj.Aofy(y, b, obj.Find_C1(Z1, Z2)), obj.find_V(s0, 0.1, Math.Pow((((b + (obj.Find_C1(Z1, Z2) * y)) * y) / (b + (obj.Find_C2(Z1, Z2) * y))), 2.0 / 3.0))), obj.Find_C2(Z1, Z2), s0, obj.Find_C1(Z1, Z2));
    tolerance = thisY - nextY;
    if (tolerance < 0.0)
    {
        tolerance = -tolerance;
    }
        thisY = nextY;
}




//F is Wrong
Console.WriteLine(string.Format("F:{0}", thisY)); ;
Console.ReadLine();

Integration and Differentiation in c#

1 solution

Solution 1

Add your solution here

Preview 0

Existing Members

...or Join us