-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathug.c
5 lines (5 loc) · 2.67 KB
/
ug.c
1
2
3
4
5
#include <math.h>
void diffeqf(double* du, const double RHS1, const double RHS2, const double RHS3, const double RHS4, const double RHS5, const double RHS6) {
du[0] = RHS2 * cos((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4) + (RHS3 * (-1 * RHS2 * 1 + (RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1)) * cos((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4)) / sqrt(((RHS5 * RHS5) * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2) + (RHS4 * RHS4) * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2) + (RHS1 * RHS1) * (RHS5 * RHS5) * pow(sin(RHS6), 2) + -2 * RHS5 * RHS4 * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2)) / pow(RHS5 + -1 * RHS4 * 1, 2)) + (RHS1 * RHS5 * RHS3 * sin(RHS6) * sin((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4)) / ((RHS5 + -1 * RHS4 * 1) * sqrt(((RHS5 * RHS5) * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2) + (RHS4 * RHS4) * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2) + (RHS1 * RHS1) * (RHS5 * RHS5) * pow(sin(RHS6), 2) + -2 * RHS5 * RHS4 * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2)) / pow(RHS5 + -1 * RHS4 * 1, 2))) + -1 * RHS1 * cos(RHS6) * cos((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4) + -1 * RHS1 * sin(RHS6) * sin((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4);
du[1] = (-1 * RHS3 * (-1 * RHS2 * 1 + (RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1)) * sin((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4)) / sqrt(((RHS5 * RHS5) * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2) + (RHS4 * RHS4) * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2) + (RHS1 * RHS1) * (RHS5 * RHS5) * pow(sin(RHS6), 2) + -2 * RHS5 * RHS4 * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2)) / pow(RHS5 + -1 * RHS4 * 1, 2)) + (RHS1 * RHS5 * RHS3 * sin(RHS6) * cos((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4)) / ((RHS5 + -1 * RHS4 * 1) * sqrt(((RHS5 * RHS5) * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2) + (RHS4 * RHS4) * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2) + (RHS1 * RHS1) * (RHS5 * RHS5) * pow(sin(RHS6), 2) + -2 * RHS5 * RHS4 * pow((RHS4 * RHS2 + -1 * RHS5 * RHS2 + RHS1 * RHS5 * cos(RHS6)) / (RHS5 + -1 * RHS4 * 1), 2)) / pow(RHS5 + -1 * RHS4 * 1, 2))) + -1 * RHS2 * sin((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4) + RHS1 * cos(RHS6) * sin((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4) + -1 * RHS1 * sin(RHS6) * cos((RHS6 * (RHS4 + -1 * RHS5 * 1)) / RHS4);
}