Экспертная система Delphi.int.ru

Сообщество программистов
Общение, помощь, обмен опытом

Логин:
Пароль:
Регистрация | Забыли пароль?

Delphi.int.ru Expert

Другие разделы портала

Переход к вопросу:

#   

Статистика за сегодня:  


Лучшие эксперты

Подробнее »



Вопрос # 6 419

/ вопрос открыт /

Доброго времени суток, эксперты!Ищу площадь фигуры, ограниченной тремя функциями, методом разбиения на прямоугольники.Среда Dev c++ , язык Си.Собcтвенно , сами функции :
f1=2^x+1 f2=x^5 f3=(1-x)/3
Чудо , что вышло у меня , не особо то и работает.В частности выдаёт ошибку "non-lvalue in assignment" в помеченных строках.
В принципе программа должна работать и для других функций , если заменить соответствующие контстанты.Tсли одна из таких функций вдруг окажется делённой на x (5/x например) , то возникает проблема.Подскажите пожалуйста в чём я не прав и где неправильно реализовал метод.
Так же как аргументировать что мы имели право использовать этот метод для решения задачи?

Приложение:
  1. #include <stdio.h>
  2. #include <math.h>
  3. #define a 1
  4. #define b 2
  5. double func1(double X)
  6. {
  7. return pow(x,5)-(1-x)/3; 2&3
  8. }
  9. double func2(double X)
  10. {
  11. return pow(2,x)+1-pow(x,5); 1&2
  12. }
  13. double func3(double X)
  14. {
  15. return pow(2,x)+1-(1-x)/3; 1&3
  16. }
  17. double func4(double X)
  18. {
  19. return pow(2,x)+1;
  20. }
  21. double func5(double X)
  22. {
  23. return pow(x,5);
  24. }
  25. double func6(double X)
  26. {
  27. return (1-x)/3;
  28. }
  29. int main()
  30. {
  31. double M, L, K, shag, EPSILON, n, k=0, S=0, i=0;
  32. printf("EPSILON\n");
  33. scanf("%lf",&EPSILON);
  34. printf("Kol-vo rasbieniy\n");
  35. scanf("%lf",&n);
  36. M=a-((b-a)/(func1(b)-func1(a)))*func1(a);
  37. while(fabs(func1(M))>EPSILON)
  38. {
  39. if(func1(a)*func1(M)<0)
  40. {
  41.  
  42. }
  43. else
  44. {
  45.  
  46. }
  47. M=a-((b-a)/(func1(b)-F1(a)))*func1(a);
  48. }
  49. L=a-((b-a)/(F2(b)-F2(a)))*F2(a);
  50. while(fabs(func2(L))>EPSILON)
  51. {
  52. if(func2(a)*func2(L)<0)
  53. {
  54.  
  55. }
  56. else
  57. {
  58.  
  59. }
  60. L=a-((b-a)/(func2(b)-func2(a)))*func2(a);
  61. }
  62. K=a-((b-a)/(func3(b)-func3(a)))*func3(a);
  63. while(fabs(func3(K))>EPSILON)
  64. {
  65. if(func3(a)*func3(K)<0)
  66. {
  67.  
  68. }
  69. else
  70. {
  71.  
  72. }
  73. K=a-((b-a)/(func3(b)-func3(a)))*func3(a);
  74. }
  75.  
  76. for(i=1; i<=n; i++)
  77. {
  78. S=S+shag*(func4(K+i*shag));
  79. }
  80.  
  81. for(i=1; i<=n; i++)
  82. {
  83. S=S+shag*(func5(M+i*shag));
  84. }
  85. printf ("%lf", S);
  86. scanf("%lf", S);
  87. return 0;
  88. }


ВикторСС. Вопрос ожидает решения (принимаются ответы, доступен мини-форум)

Вопрос задал: ВикторСС. (статус: Посетитель)
Вопрос отправлен: 6 мая 2013, 13:38
Состояние вопроса: открыт, ответов: 0.


Мини-форум вопроса

Мини-форум пуст.

Чтобы оставлять сообщения в мини-форумах, Вы должны авторизироваться на сайте.

Версия движка: 2.6+ (26.01.2011)
Текущее время: 18 июня 2019, 02:11
Выполнено за 0.02 сек.
Рейтинг@Mail.ru