Réponse apportée
2D information to 1D colorplot
Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

plus de 6 ans il y a | 0

Réponse apportée
Quiver with a system of ODEs.
For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

plus de 6 ans il y a | 0

Réponse apportée
Solving nonLinear equation and Linearized Eqaution using ODE45
YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

plus de 6 ans il y a | 0

| A accepté

Réponse apportée
How to plot magnitude and directions between 2 points?
try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

plus de 6 ans il y a | 0

| A accepté

Réponse apportée
surface fixed point value Z
I want to know the value of Z from the surfc function Here is it

plus de 6 ans il y a | 0

Réponse apportée
Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?
Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

plus de 6 ans il y a | 1

| A accepté

Réponse apportée
I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?
alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

plus de 6 ans il y a | 0

Réponse apportée
Euler 3D rotation between two vectors
Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

plus de 6 ans il y a | 0

| A accepté

Réponse apportée
Simultaneous differential equations - derivative of y^3 wrt t
Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

plus de 6 ans il y a | 0

Réponse apportée
Numerically derive a continous, non-symbolic function
Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

plus de 6 ans il y a | 0

| A accepté

Réponse apportée
3D Point Cloud - Gaussian Curvature
one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

plus de 6 ans il y a | 0

| A accepté

Réponse apportée
Runge Kutta for system of eqs
I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

plus de 6 ans il y a | 0

Réponse apportée
The problem of finding intersection point of triangle mesh
To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

plus de 6 ans il y a | 1

| A accepté

Réponse apportée
How can i solve a system which contains matrices
Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

plus de 6 ans il y a | 1

| A accepté

Réponse apportée
Using a for-loop to score values in different intervals
Correct form if 0 <= x(i) && x(i) < 1

plus de 6 ans il y a | 0

| A accepté

Réponse apportée
Intersection between line and circle when line ends inside circle (using geom2d)
Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

plus de 6 ans il y a | 1

| A accepté

Réponse apportée
Removing points from meshgrid (mask)
Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

plus de 6 ans il y a | 0

Réponse apportée
How can I find the surface normal of a 8 by 3 matrix?
Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

plus de 6 ans il y a | 0

Réponse apportée
How do I increase the circle size and animate a stick figure?
Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

plus de 6 ans il y a | 0

| A accepté

Réponse apportée
Finding specific data point x for given y
I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

plus de 6 ans il y a | 0

Réponse apportée
Surface between 3D lines
You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

plus de 6 ans il y a | 2

| A accepté

Réponse apportée
two body problem using ode45
- equation of equilibrium Is this condition fullfield?

plus de 6 ans il y a | 0

Réponse apportée
Numerically Integrating the differential equation below?
Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

plus de 6 ans il y a | 0

Réponse apportée
Streamline plot not complete
Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

plus de 6 ans il y a | 1

| A accepté

Réponse apportée
2D Matrix revolution
Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

plus de 6 ans il y a | 1

Réponse apportée
How do I align the vertices and the adjacency graph of a Buckyball?
Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

plus de 6 ans il y a | 0

Réponse apportée
Find the coordinates of a point chosen on a plot
Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

plus de 6 ans il y a | 3

| A accepté

Réponse apportée
How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?
Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

plus de 6 ans il y a | 0

Réponse apportée
Differential Algebraic Equations with vectors
Try f = @(t,u) inv(Z)*(J+p*W*pi*Q)*[u(3) u(4)]' + inv(Z)*(V+p*W*pi*Q)*[u(1) u(2)]'; F = @(t,u) [u(3) u(4) ...

plus de 6 ans il y a | 1

| A accepté

Charger plus